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PREFACE. 



The present Treatise is intended a.s an introduction to 
the study of the DiflTerential and Integral Calculus, but 
will be found to contain what it is necessary to know in 
order to pass on to the subjects which presume a know- 
ledge of the Calculus. 

I have endeavoured to make this book suitable not 
only for the mathematical student, but also for men like 
engineers and electricians who require the subject for 
practical applications, to whom even a slight knowledge 
of the notation and methods of the Calculus is becoming 
more and more indispensable. 

Hitherto in this country the influence of Newton, 
although the inventor of Fluxions, has been employed 
to delay the study of this subject and make a know- 
ledge of it the privilege of a select few ; my object in 
writing this treatise has been mainly to present the 
subject in as simple a manner as possible, in order to 
encourage a larger number of students to cultivate it. 

With the object of keeping the size of the book within 
reasonable limits, it is assumed that the reader has 
already acquired a knowledge of the elements of Algebra, 
'Trigonometry, and Coordinate Geometry, as given, for 
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vi PREFACE. . 

instance, in the treatises of Hall and Knight, J. B. Lock 
and C. Smith ; accordingly I have at once proceeded to 
the operation and application of Differentiation with as 
little preliminary explanation as possible. 

I have followed the recent American treatises on this 
subject of Rice and Woolsey Johnson, Byerly and J. H. 
Taylor, in introducing the notion of Time as an inde- 
pendent variable, and the associated ideas of velocity 
and acceleration, in order to afford illustrations of the 
use of the Calculus ; this is after all only a return to the 
Method of Fluxions as invented by Newton, and carried 
out by Maclaurin and other writers in this country, 
until supplanted by the notation of the Differential 
Coefficients of the foreign mathematicians. 

The Doctrine of Fluxions is a useful and rigorous 

method of presenting the elementary ideas of the flow 

of varying quantities, and is employed in the treatises 

of Rice and Woolsey Johnson under the name of the 

Method of Rates; but the notation for a fluxion, for 

instance x the fluxion of x, though easily written is 

difficult to print, and has the inconvenience of not 

indicating the independent variable, so that the notation 

dx 
of Leibnitz, jt- instead of Xy is now used almost uni- 

dut 

versally in printed books ; and to economise space, this 

notation it is now proposed to print in the form dx/dt 

The chief novelties in the present work consist, first, 

in carrying on the subjects of the Differential and of 

the Integral Calculus together, instead of, as is usual, 

completing the Differential before passing on to the 

Integral Calculus ; secondly, in the use of the hyperbolic 

functions in conjunction with the ordinary circular 
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trigonometrical functions, as thereby an exact analogy is 
preserved, which is not apparent when only the exponen- 
tial and logarithmic functions are employed. 

The notation of sinh, cosh, tanh, etc., to denote the 
hyperbolic sine, cosine, tangent, etc., has been employed, 
in accordance with what appears to be now the most 
universal custom. 

I have ventured also, on the grounds of symmetry, 
to introduce the inverse hyperbolic functions, and, fol- 
lowing Ferrers and Byerly, to denote them by sinh"^, 
cosh-^ tanh"\ etc., by analogy with sin-^ cos~\ tan"^ 
etc.; this idea of symmetrical symbolism will be found 
indicated in Bertrand's Integral Calculus, Chapter I., but 
has not been pursued apparently because of the lengthi- 
ness of the notation there employed, namely, sect, sin 
hyp., sect, cos hyp., sect, tang hyp., etc., instead of the 
above. 

By the use of the direct and inverse hyperbolic func- 
tions in conjunction with the direct and inverse circular 
functions, the Calculus is in my opinion considerably 
simplified; and the student is led on more naturally and 
readily to the consideration of the elliptic and other func- 
tions. The consideration of these last functions is how- 
ever beyond the scope of the present treatise. 

To exhibit more clearly the analogy and symmetry 
between the circular and hyperbolic functions, I have 
made a digression in Chapter I. on the formulas of the 
addition equation (as it may be called by analogy with 
elliptic functions) of ordinary trigonometry, showing 
how the formulas may all be deduced from a single 
figure, with the corresponding relations of the hyperbolic 
functions. 
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Numerous collections of examples will be found 
throughout the book, introduced at each point to illus- 
trate what has immediately gone before, and chosen as 
having some bearing on subsequent stages of mathematics. 

The order of arrangement will be found in some 
respects different to what is customary ; for instance, 
the idea of tracing simple curves from their equations 
has been introduced into the first chapter, so far as is 
required for the ordinary applications of the Integral 
Calculus to finding the areas, etc., of these curves; the 
general theory of curves being resumed in Chapter VI. 
Maximums and Minimums also have been investigated 
without the aid of Taylor's Theorem. 

Change of the Independent Variable has only been 
touched upon where necessary ; the general theory of 
Change of the Independent Variable, Lagrange's and 
Laplace's Theorems, and the Elimination of Constants 
and Functions, have been omitted as beyond the scope 
of an elementary practical treatise. 

I have to thank Mr. A. G. Hadcock, Inspector of 
Ordnance Machinery, Royal Artillery, for drawing the 
diagrams, and also for revising the proof sheets and 
preparing the index. 



Woolwich, 
December, 1885. 
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PREFACE TO THE SECOND EDITION. 



In this second edition a rearrangement has been made, 
by which the general theory of Integration has been 
relegated to a supplementary final Chapter; ollly a 
slight preliminary sketch of Integration is now given 
in Chapter II., sufficient to carry on the student to some 
simple applications, and to enable him to follow the 
subsequent articles. 

In this rearrangement, and in many other details at 
the outset, I have received valuable advice and assistance 
from Mr. W. Gallatly, M.A. 

The theories of Change of the Variable, the Elimina- 
tion of Constants and Functions, and the Theorems of 
Lagrange and Laplace have now been included, so that 
this new edition is increased in size. 

By following the custom of giving the examples in 
smaller print, this increase in size might have been 
apparently much reduced; but it was decided to use 
a uniform tjrpe throughout the book, as less trying to 
the eyesight. 

The solidibs notation, for instance dx/dt for -,y, has 

w been used sparingly, in places where no confusion 
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or ambiguity could arise, and thereby some economy of 
space has been obtained. 

A few pages of explanation of the simplest Differen- 
tial Equations which occur in Dynamics and Electricity 
have been added; while the illustrations of Change 
of the Variable have been chosen from among the most 
typical differential equations; and it is hoped that this 
slight sketch will enable the student to solve an ordin- 
ary differential equation which may come in his way. 

The greater part of the diagrams have been redrawn 
by Mr. A. G. Hadcock, who has again helped me in the 
revision of the proof sheets and the preparation of the 
Inde^ 



Woolwich, 

March, 1891. 
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DIFFERENTIAL AND INTEGRAL CALCULUS. 



CHAPTEK I. 

DIFFEEENTIATION. 

1. Introduction, Constant and Variable Quantities. 

The Calculus to be developed in the present treatise is 
the method of reasoning applicable to variable quantities 
in a state of continuous change. 

We call quantities variable when they change gradu- 
ally by increase or decrease ; on the contrary constant, 
when they remain unchanged while others change. 

Thus the abscissa and ordinate of points on a parabola 
vary, while the parameter remains constant. 

Again, the distance of a railway station from the 
terminus, or the latitude and longitude of a rock, of a 
lightship, or of an observatory, are constant quantities. 

But the distance of a railway train from the terminus, 
or the latitude and longitude of a steamer or traveller 
are variable quantities. 

Constant quantities are generally represented algebra- 
ically by the first letters of the alphabet, such as a, 6, c, 
..., or -4, JB, (7, ..., or a, /3, y, ...; while variable quantities 
are represented by the last letters, ... t, u, v, w, x, y, z, 
or ... X, Y,Z, or ... ^, ^, ^. 



2 DIFFERENTIA TION. 

2. Definition of a Function ; and of Independent and 
Dependent Variables. 

One variable quantity, denoted by y or f(a;), is said to 
be a function of another variable quantity, denoted by Xy 
when the value of y or f (a:) depends on the value of x. 

The variable x^ called the independent variable or 
argument, is one to which any value may be arbitrarily 
assigned. 

The variable y or i{x), called the dependent variable, is 
one of which the value is determined, when the arbitrary 
value of the independent variable x has been assigned. 

We may use the notation ix instead of i{x) when the in- 
dependent variable consists of a single term like x. 

Thus x^, x^, x^, ..., x^, j^x, ;^x, sin a:, cos a?, tanaj, cot a;, 
sec a:, coseca:, versa;, sin"^a;, cos"^aj, tan"^a;, cot"^2J, sec"^aj, 
cosec"^a;, vers"^aj, a^, exp x, cosh x, sinh x, tanh aj, ... log x, 
cosh"^a;, sinh"^aj, tanh'^a;, ..., all denoted generically by 
fa;, and most of them presumably familiar to the student, 
are simple functions of x, which we shall require hereafter. 

We proceed to differentiate them, that is to find the 
differential coefficient, which is defined as follows : — 

3. Definition of a Differential Coefficient, 

If fa? denotes any function of a variable quantity x, 
and if f(a;+/i) denotes the same function of a;+/i when x 
receives a small increment h, then the limiting value of 

i(x+h)-{x 
h 

when h is indefinitely diminished, is called the differential 
coefficient of fa; with respect to x, and is denoted by 

dix n, 

—r- or I X. 
ax 
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DIFFERENTIATION. 3 

This definition may be conveniently expressed as 

dix _ , , f (a; + A) — fic 
dx'~ h 

(It) being the abbreviation employed to denote the limit- 
ing value, as A is indefinitely diminished and ultimately 
becomes zero. 

Since i(x+h)^{x is the increment of ix corresponding 
to the increment h of x, therefore 

^* — h — 

is the ultimate ratio of the corresponding increments of 

tx and X, denoted by dix and dx, and called the differ- 

dix 
entiaZs of ix and x ; and thus -r— measures the rate of 

increase or growth of ix ; while -^^ ^ represents 

the average rate of increase of ix from xio x+h. 

The chief object at the outset of our subject is the 
determination of the differential coefficients of functions, 
and the application of them to the discussion of the 
geometrical and analytical properties of the functions. 

The algebraical difficulty in the determination of the 
differential coefficient lies in the reduction of its original 

indeterminate form tt to a determinate limit. 

(Hall and Knight, Higher Algebra, chap, xx.) 
The name derivative or derived function is sometimes 
used instead of differential coefficient. 

The differential coefficient i'x of a function of ix may 

be supposed to derive its name from being the coefficient 

which turns the differential of x into the differential of ix. 

(Wicksteed, The Alphabet of Economic Science, ip, 32.) 



4 DIFFERENTIATION. 

4. Differential Coefficient or Derivative of x^. 

From the definition of a differential coeflScient in § 3, 

dx h 

dx^^^^ jx+hf- x^J) 
dx ho 

.-2xh+h^ 

= ^^— A-" = 0' 

= li(2x+h) = 2x. 

dx~ h ~ 

.,^xVi + ^xh^ + h?_0 
"^* h "6 

= \i{^x^+^xh+h^) = ^x\ 

da^_ J,x+hf-x^ _0 
dx^ h ~ 

_ .Aa?h + 6x%^ +4x^4- h^ _0 

And generally, when n is a positive integer, 

dx''_. (x+hy-x'' _0 
dx~~ h ""0 

(and expanding by the Binomial Theorem) 

nx^-%+^-J^x''-Vv'+ ..,+h^ ^ 
"^* h "0 

By assuming the convergency of the Binomial Theorem 
we can make the same proof hold for showing that the 
differential coefficient of x^ with respect to x is nx^~^^ 
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where n represents any constant number, whether 
integral, fractional, positive or negative. 

But without using the Binomial Theorem or assuming 
its convergency, a proof may be given as follows : — 

(i.) Suppose n a positive integer, and denote x+h 
by 0?! ; then 

dx x^—x 

(and then by Division) 

=lt(x^^-'^+xx^''-^ + ,,.+x''-^x^+x'^''^) = nx'^-\ 

when h = 0, and therefore x^ = x, 

(ii.) Suppose oi a positive fraction p/q, where p and q 
are positive integers; and put x = z^y x^ = z-^\ then 

dx^ = It^i^Z -"" ? — 

(mX fl/j ~" X 

(dividing out z^ — z from numerator and denominator) 

_j^ V+V"^2;+ ... -\-z^z^-^-\-z^-^ 
z^-^-\-z^-H'\- ... +2;i0^-2+2?«-i 

qz^-^ q 
(iii.) Suppose n any negative number — m; then 

dx^ y,x^-'*^—x-'^ ^^x.^—x^ 1 



= lt-l ^ = -lt-1 



'mnt*'n\, 



dx X^ — X X^—X x^^od 

=z —7rix-'^-^ = nx^'^: 



mx'^-^ 



X' 



dx^ 
so that -T- = nx^~^, universally, where n denotes any 

constant quantity. 



DIFFERENTIA TIOS. 



—1. Determine from the definition the d.c. 
(differential coefficient) with respect to x of 

■J"- «'■ ^- '^' "=*■ '''• ^' \ S' i' S" '"+«)" 

(£)",?£i_3^,„+6«+rf+... . 

2. Give a rigorous proof that the differential coefficient 
of x-^ with reapect to a; ia ^2 a;^-^. 

5. Geometrical Interpretation of Differentiation. 




In Figure 1 the coordinates of a point P, referred to 
axes Ox and Oy at right angles, are denoted by x and y, 
where OM=x, MP=y; and in general, in the vertical 
elevation of an object, x is measured positive to the right, 
and negative to the left, y is measured positive upwards, 
and negative downwards; or, on a map or plan, Ox is 
drawn to the east, and Oy to the north. 



DIFFERENTIATION, 7 

Then the equation y = ix represents some curve APB, 
the assemblage of points whose coordinates satisfy this 
equation; so that if OM=x, then MP = ix; the curve 
APB is now called the graph of the function fx 
(Chrystal, Algebra I., chap. xv.). 

If ilfm = /i, then Om = x+h, 

and mp = i(x+h), Rp = i{x+h)'-ix; 

so that i(^ + h)-ix ^ ;fi2 = tan RPp. 

h JrK 

Now if the angle xTP is denoted by ^, then since the 
direction of the tangent TP is the ultimate direction of 
the chord Pp when the point p has approached in- 
definitely neJEir to P, therefore 

tan ^ = It tan RPp = Itp^ 

'..i(x+h) — ix dix „, dy , ,, 

= ^^ h =-dx' °' *^' °' di' °' ^^z*^' 

since 2/ = ^ ^• 

By giving a? a decrement h as well as an increment A, 
we may define a differential coefficient by the relation 



dix 



_^l(x+h)'-i(x-h) 



dx 2h 

which is sometimes more elegant and useful than the 
definition in § 3, to which it is ultimately equivalent. 

Figure 1 shows that dy/dx or tan xfr is positive if y 
increases with x, but that dy/dx is negative if y dim- 
inishes as X increases; and thus dy/dx is zero at the 
turning points, such as £ or C, where y from increasing 
begins to diminish, or from diminishing begins to increase. 
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Ki B^y is said to have a Tnaxmrnvfi value, and at C a 
minimuTn value ; and in each case dy/dx=^0. 

Thus to discover the maximum or minimum value of 
y, a function of x, we must first find the values of x 
which make dy/dx = 0. 

Afterwards we must determine whether, as x increases 
continually, dy/dx changes from positive to negative, in 
which case y has a maximum value ; or whether dy/dx 
changes from negative to positive, in which case y has 
a minimum value ; we shall return to this subject here- 
after, although the present application to simple cases 
is easy. 

Thus ax^x^ has a maximum value Ja^ when x = \a 
(fig. 33); and 2a3^— 9aj2+12a; — 3 has a maximum value 2 
when x=\y and a minimum value 1 when a; = 2 (fig. 29). 

Imagine the curve ABG to be a road, and dx a small 
horizontal forward step, dy the consequent vertical step ; 
then dy and therefore dy/dx is positive if the traveller 
is ascending, and dy and dy/dx is negative if the traveller 
is descending, supposing dx always positive. 

At the turning points, such as the top of a hill or the 
bottom of a valley, where y is a maximum or minimum, 
dy/dx = 0. r 

The road is steepest where dy/dx has its maximum 
numerical value, and then the differential coeflScient of 
dy/dx will be zero. 

We may call tan i/^ or dy/dx the gradient of the road, 
the ascent being dy/dx feet vertical for one foot hori- 
zontal ; or otherwise expressed, as on railway plans, the 
gradient is 1 in cot y\t or dx/dy, meaning 1 foot vertical 
for dx/dy feet horizontal. 
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6. The Tangent and Normal of a Curve. 
If x\ y' are the coordinates of any point P' on the 
tangent TPF to the curve APB at P (fig. 1), then 

x^x dx ^ ^ dx^ ^ 

the equation of the tangent TP. 

TM is called the subtangent, and MQ the subnormal 
at P ; P(? being the normal at P, drawn perpendicular 
to the tangent through P. Therefore 

TM=y cot i/r = ydxjdy ; if (? = y tan i/r = ydy/dx. 
If T is to the right of M, tan yfr is negative, and 

MT= -- ydxjdy, QM= —ydxjdy. 
Also the equation of the normal OP at P is 

, dXf , ' . % 

Let us apply these principles to some simple curves, 
the graphs of (i.) y = x\ (ii.) y = x^y (iii.) y = x'\ 
(iv.) y^ = a?. 

(i.) y = x^y a parabola (fig. 2 i.). 

Here dyjdx = 2x ; so that the equation of the tangent is 
y'—y — 2x{x''-x) ; which, since y = x^, may be written 

Thus Or=iaj, TiJf =|a;, and therefore T bisects Oif. 
The equation of the normal is ^ o — +-^71 ~ ^' 

(ii.) y = x^,'ajioiher parabola (fig. 2 ii.). 
Here dyjdx = |aj - i = Jy/a; ; 

so that the equations of the tangent and normal become 

I2/ a? ' 0; + ^ 22/^+2/ 
thus TO = OM, 0V= F/V, and i/G = i. 



10 
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M G X o\ M Tx 

(ii) (iii) 

Fig.2 

(iii.) 2/ = a;-\ a hyperbola (fig. 2 iii.). 
Here ^^=—03-2= — 0=--; 

dx OCT X 

and the equation of the tangent becomes 

2x^2y ' 
so that OT=2x,OV=2y, 

The equation of the normal is 

icx'-'yy'=x^—yK 
(iv.) y^ — Qi?, or y=aj^, a curve called the aemicubicat 
parabola, and something like figure (i.). 

Here dy/dx=^x\ so that TM=^x, OT=ix = iOM, 
(v.) Prove that in the curve y'^ = cx'^, the equations of 
the tangent and normal are 

^' 2/' 1 1 / . / 2.2 
^- = ", Tnxx +nyy =7nx^+ny\ 

rax ny m n 

* 7. Supposing h again, as at first, to denote a finite 
increment of x, and supposing that neither ix nor i'x 
becomes infinite between P and p, and that they change 
graduaUy between these points ; then in the graph of {x 

* Articles which may be omitted at a first reading are marked 
with an asterisk* 
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the tangent TK at some point K between P and p 
(fig. 3) is parallel to the chord Ff ; and 

!(^+^)^?=taniiP2) = tanajr^=f>+eA^), 

where x + 6h=0H, the abscissa of K; and then 
6 = MB /Mm, a proper fraction, which is some un- 
known function of x and h. 




Therefore i{x +h) = ix+ M'{x + Oh), 
a theorem required subsequently in Chapter IV. 

This leads to another definition of fx, employed by 
Lagrange, according to which i'x is defined as the co- 
efficient of h in the algebraical expansion of i{x+h) in 
ascending positive integral powers of h. 

To illustrate geometrically exceptional cases of tx 
and fx, consider ix = (x — b)/{x — a), or ix = {X'- a)^, taking 
x<a, and x+k > a. 

Given that ix = a+bx+cx^, 

then 6 = ^; and generally, when h is small, we shall find 
0^^ ; the symbol ^ denoting approximate equality. 
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lik denotes a different increment of x, then 

i{x +h)^ix+ ld\x + </>k), 

where <f> is a proper fraction, the same function of x and 
k which is of aj and h ; so that 

i(x+h)-fx _h r(x+eh) 
i(x+k)'-ix ki\x+<l>k)' 

Ultimately, when h and k are suflGlciently small, we 

may put 

i\x+dh) _. 

r(x+<f>k) ' 
without sensible error ; and then 

i(x+h)-'ix_h 

i(x+k)''ix k' 
which is the rule of jproportional parts, employed in 
Trigonometry ; equivalent to neglecting the curvature of 
the arc Pp, and replacing the arc by a straight line. 

Again supposing APp to represent the profile or 
elevation of a road, and i'x to represent the slope or 
gradient at P, i\x+h) at p, then the gradient i\x+Qh) 
at K is the average gradient of the chord Pp, repre- 
sented by {i{x+h)^ix}jh 

8. Instead of the letter h the symbol Aaj is often em- 
ployed to denote the finite increment of x ; and Aj/ is then 
used to denote the corresponding increment of y, where 
2/ is a function of a?, denoted by ix. 
Then, with this notation, 

y+£^y = i{x+^x\ 
and A2/ = f(iB+Aa:) — faj; 

so that ^y = Up, if ^x = PR (fig. 1) ; 

and ^ = lt^. 

dx Ax 
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Li this notation, A (Delta) and d are Hot algebraical 
factors, but symbolic operations ; so that Ace and dx must 
be considered as inseparable quantities, the x being quali- 
fied by the A or d, and not multiplied algebraically by it. 

Suppose for instance it is found that the difference of 
range of a gun over the sea at high and at low water is 
Aic, when Ly represents the rise and fall of the tide ; then 
the tangent of the angle of descent of the projectile into 
the water is Ay/ Ax, approximately. 

9. Again, let the length of the arc AP oi the curve 
APp, measured from any fixed point A to the variable 
point P, be denoted by s (fig. 1). 

Denoting by As the increment of s, corresponding to 
the increment Ax of x, then the arc APp =^8+ As, and the 
arc Pp = A8. 

Now, when the point p approaches to coincidence with 
P along the curve pP, it is assumed as axiomatic that 

uchordPp 1 

It D-- =1' 

arcPp 

arc Pj) 

As 
where Ax^ means (Ax)^ ; 



or, as it is sometimes written, 

it)'- 




2 

+ 
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The assumption, which has been taken as axiomatic, 
that It (chord Pp)/(arc P/)) = l, is proved in Newton s 
Lemma VII. (Principia, Book I, § 1), as follows — 




Fig.4 

As the point p approaches to the point P, let the chord 
Pfy and the tangent Pv be produced to points Q, F, at a 
finite distance, where QF is parallel to pv, and the arc 
PQ is similar to the arc Pp; so that PQF may be con- 
sidered the magnified image of the similar microscopic 
figure Ppv, 

Now let the point p move on the curve pP close up to 
P, and let pv move parallel to itself up to P ; then the 
angle QP V always diminishes and ultimately vanishes ; 
and the chord PQ, the arc PQ and the tangent PV are 
ultimately coincident; therefore also the chord Pp, the 
arc Pp and the tangent Pv, when vanishing together, 
will have a ratio of equality. 

Again (fig. 1) cos xJa = cos xTP = It cos RPp = It- 



= lt^lt 



arc Pp _ dx . 



Pp 



and 




Since 
therefore 



As chord Pp ds ' 
sin yfr = sin xTP = It sin RPp = It 

= li^y It ^^^Pp = ^. 

As chord Pp ds' 
(cos ^)2+ (sin yfrY = 1, 
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The tangent rP = 2/cosec yfr^yds/dy ; and the normal 
PO = y sec \[r = yd8/dx; and the length of the perpen- 
dicular p from the origin upon the tangent at P is 

. , , dy dx 

p = X8inyfr-ycos\lr = x-^-y^. 

Since (sec yf/f = 1 + (tan \l/f, 

and sec yjr = ds/dx, tan xJa = dy/dx ; 

ds du 

therefore ^ = l+5&- 

Similarly, from (coseci/r)2 = (cotV^)^+l, 

T -i W'O Q/X . 1 

we deduce ^— «= -7-0 + 1- 

dy^ dy^ 

Sometimes the gradient is measured by the vertical 
rise for one foot on the incline, that is by dy/ds or the 
sine of the slope yfr; and then the gradient is one in 
ds/dy or cosec xJa, 

We may say that the gradient is measured by the 
tangent of the slope on the plan, and by the sine of the 
slope on the constructed work ; in the inclines of roads 
and railways the slope is sufficiently small for no differ- 
ence to be perceptible in the two methods. 

10. Velocity expressed by a Differential Coefficient. 

If X and y, the coordinates of a point P, moving along 
the curve APp (fig. 1) ; the path of a projectile in the 
air or of a train on a railway, for instance) are given as 
functions of the time t; and if A^ denotes the time 
occupied in moving from P to j9, a distance As ; then the 
average velocity from P to jp is As/ At; and the actual 
velocity of P, in the direction of the tangent TF, is 

i,As _ds 
Ai'dt' 
the rate of increase or growth of s per unit of t 



C' 
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The component velocity of P parallel to Ox is the 
velocity of M\ and treating Jf as a point moving in the 
straight line Oaj, the average velocity from M to m is 
Lx\Lt \ and therefore the actual velocity of M 

_yA.x__dx 

the rate of growth of x per unit of t 

Similarly the component velocity of P parallel to Oj/, 
or the rate of growth of y is dyjdt, the velocity of N, 

Since the resultant velocity of P, dsjdt, has the com- 
ponents dxjdt and dyjdt in two rectangular directions, 

., p cfo^ dx^ . dy^ 

therefore M^W^W' 

Making t = 8, x, or y, we obtain the previous results. 

According to Newton's original definition, the fluxion 

of a variable quantity x is its velocity or rate of growth, 

and was denoted by i; in Leibnitz's notation of the 

dx 
Differential Calculus it is represented by -77 or dxjdt \ 

thus diflerentiation with respect to the time t may be 
considered the typical fundamental conception of this 
Calculus, which investigates the rate at which variable 
quantities change, now not only with respect to the time 
t, but also with respect to any other independent variable. 
11. Differentiation of a Function of a Function. 

CI. Ay Ay Ax 

A^ Ax At 
therefore, proceeding to the limit, 

dy__dy dx 

dt dx dt' 
where a? is a function of t, and y is a, function of x, and 
therefore of t 
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Thus, with the notation of §§ 9, 10, (fig. 1), 

^_^dx ds^ da dy _dy ds __ . .da 

~dt~d^dt-''^^'l'df dt^dadr^^'^'^di' 

showing how the component velocities are obtained 
by resolving the resultant velocity. 

Denoting the velocity da/dt by v, then dv/dt, the rate 
of growth of V, is called the acceleration; and 

dv _dv da _ dv __ d^v^ 

dt da dt da da * 

With different letters, supposing is a function of y, 

and 2/ is a function of the independent variable x ; and 

supposing the increment Ax of x gives y the increment 

Ay, and z the increment Az ; then since, algebraically, 

Az_ Az Ay 
Ax Ay Ax* 

therefore * -^^=^ % 

ax ay ax 

the rule for the differentiation of a function of a function. 
Thus, suppose z = y^, where y is a function of x ; then 

dz ^dy"^ _dy'^ dy _ .^dy , 
dx dx dy dx dx^ 

or,if 2/ = fa;, then ^:!^^=m{ixy^-Wx, 

ax 

Putting z = Fy and y = ix, so that z = F(ix); then 

dz^^ <^ = ry rx = F(ix) i'x, 
dx dy dx 

dz 
Thus if = (a -h hx"")^, -r- =m(a + bx'^)'^-'^nbx''-\ 

With 2/ = faj, then dy/dt = fx,dx/dt; and this, in the 
notation of Differentialay is abbreviated to dy = fx,dxy 
in which the independent variable is left arbitrary ; and 
dx, dy are called the differentials of x and y, as in § 3. 

B 
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12. Differentiation of the Sum, Product, and Quotient 
of Functions, 

Let u and v denote any two functions of x; and let 
Au and Av denote the increments of u and v, correspond- 
ing to the increment Ax of x ; then (i.) for their sum, 

d(u+v) _. Au+Av _../Au Av\_du dv ^ 
die "" Ax " \Ax Ax/^dx dx' 

so that the d.c. of the sum of two functions is the sum of 
the dx.'s of the functions. 

Similarly, for the difference of u and v, 

d{u-'V) _du _dv 
dx dx dx' 
And generally, if a and h denote constant factors, 

d{au + hv) _ du .dv 
dx "~ dx dx 

Thus ^(ax'^+bx'''-cxP+.,.) 
dx 

and if s = ut+\ft^, v=ds/dt=u+ft 

We notice now that the differential coefficient of a 
constant is zero ; as is obvious, since a constant does not 
change with the independent variable. 

(ii.) For the product of u and v, 

duv _. (u+Au)(v+Av)-'UV 
dx " Ax 

ixAu , ,, . , . .Av du . dv 

because It^ = -r-, It^ = -r-, and lt(u + Au) = u. 

Hence the rule for differentiating a product: Differ- 
entiate with respect to each factor and add the results. 
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For instance if u = (ic — a)'^, v = (a; — 6)**, 

= (ic— a)*^-i(aj — 6)"-^{(m+'Ji)ic— m6— net}. 
It is important for subsequent purposes not to alter the 
order of the factors in differentiating ; thus for a product 
of three factors u, v, w, 

duvw du , dvw du , dv , dw 
ax ax dx dx dx ^ dx 

and so on for any number of factors. 

To illustrate geometrically, refer to fig. 1, in which 
P is a point moving along the curve AP, whose co- 
ordinates are x, y at the time t 

In the time A^, P has moved to p, whose coordinates 
are x + Aaj, y+^y \ and the rectangle OMP^ = xy has in- 
creased by the rectangles Prn = Aoj. 2/,and nR =-{x+ Ax) Ay; 

and the formula 

dxy dx , dy 

dt "dt^^ dt 

expresses the rate of growth of the rectangle xy in the 
two directions Ox and Oy. 

(iii.) For u/v, the quotient of n by v, 
,u u+Au u 

V _ ., v+Av V __ 1 . ('M' + Au)v — u(v + Av) 
dx~' Ax "~ Ax{v+Av)v 

All Av dvb dv 
.— -y— u-r— -^v—u^f- 

= lt~^— _^^ — • 

(v+vA)v "" v^ ' 

hence the rule for differentiating a fraction: the dif- 
ferential coefficient is equal to 

d.c. of num\ x denom^. — nvrni'^. x d.c. of denom'^, 

(denominator^ 
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Thusif,=f^3j 

<^y- (3-3a;^)(l--3a;^)-(3a?-g;g)(--6a;) _^/ l+g;g y 
dx (l-3aj2)2 \l-3ic2/' 

If the denominator of the fraction is of the forno 
-y", for instance, if y = u/if^, 

du ^ ^ .dv dub dv 

then ^ — ^^ dx_dx dx 

dx ^ If^ ' 

after removing the common factor v'^-^ from the numer- 
ator and denominator. 

Thusif 2/ = 5^^^#, 

dy _ m(x'-a)'^-'^{x—b)'*—(x'-a)''^(X'-b)'^-^ 
d^ (aj-6)2« 

In such cases it is sometimes preferable to write the 
fraction in the form of a product ; then y = uv-^; 

««.! dy du ^ ^ .dv 

and ^= -y-n_^^^-n-i 

ax dx dx 

du dv 

_dx dx 

as before, but without introducing extraneous factors. 
Thus ii y = ^-, r— , 

writing it in the form 

y^ix — a) -^(x — 6;-^ 

then ^= -m(a;-a)-'^-Xa;-6)-«-(aj-a)-'"'n(a;-6)-'*-^ 
dx 

_ (m+n)aj— ?ia— m6 
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Examples of Differentiation. 
(1) Prove the following differentiations 

(I) !' = VW,|-^^,. 
{Solution,— 'Kqtq y = (2a;)* and therefore 

|=fx(2.)-ix2=(2.)-i=^^.) 



i-i 



(ii.) y=(«a;)» ^=(«a;)- 

,. X 1 civ 1 

(IV.) 2/ = — -^ = -. 

s/{a-x) dx 2(a-a;)^ 

(V.) y-s/(a^+-^i'£-^^y 

(vi) 1,= I dy__x__ 

(vii.) 2/= ? ^-._^' 



V(a32+a2)' da; (x^+a^f 
(Solution. — Employing the rule of § 12, 

d^ x^+a^ (a;2+a2)ff/ 

(viii.) y = ^(2aa; - a;^), -J^- = //f~^ or 
^ / i7 'v \ /* ^^ ^(2ax — x^) 

(2) Differentiate with respect to a;, 
|a;7-faV+a*a;^-a«a;, x^(l-x% x{l-2x^)J(l-x^), 
a — x la^x 2x 3a;+a;^ 4a; — 4a;^ 5a;+10a;'+a;^ 



V^ 



a+x \a+a; l-x2 l + 3a;2'l-6a;2+a;* l + 10ai2+5a;* 
(3) Deduce the rule for the differentiation of a quotient 
from the rule for the differentiation of a product. 
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13. Implicit and Explicit Relations between two 
variables, x and y. 

When the variables x and y are connected together 
by any relation, for instance by 

ax^ + 2hxy + by^ + 2gx+2fy+c = 0, 

the general equation of a conic section, or by 

the equation of a cubic curve, this relation is called an 
implicit relation between x and y, and is denoted gener- 
ically by the notation 

F(x,y) = 0..: (A). 

When however it is possible by solution of the equation 
to obtain y in terms of x, or x in terms of y, then y is 
called an explicit function of x, or x of y, denoted 
generically as before by y=tx, or x = </>y (in the notation 
of inverse functions, § 25, x=i~^y). 

In the above implicit relations between x and y, we 
cannot obtain y explicitly in terms of x, or x in terms 
of y, except by the solution of a quadratic or cubic equa- 
tion ; but, for instance, from the implicit relation 

x^y^'-a\x^—y^) = 0, 
the equation of a quartic curve, we obtain explicitly, by 
mere transposition, 

^ ^a^+x^* ~a^'-y^' 

From an implicit relation (A) between x and y to find 
dyfdx, it is not necessary to express 3/ as an explicit 
function of x ; but we differentiate (A) with respect to x, 
treating 2/ as a function of cc, and obtain what is called 
the^rs^ derived equation and thence determine dy/dx as 
a function of x and y, without solving the equation (A) ; 
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thus from the implicit relation or equation of a conic 
section we obtain the first derived equation. 

2aaj+2%+2/«cg+262/||+2flr+2/J=0. 

so that dy_ ax+hy+g 

dx hx+by+f 

For changing the coordinates x, y ot a. point on the 
conic section to x+Ax, y+Ay, the coordinates of an 
adjacent point on the curve, then 

F(a;+Aa5, y+Ay) = 
a(x + Axf + 2h(x + Ax)(y + Ay) + b(y + Ayf 

+ 25f(a; + Ax) + 2/(2/ + Aj/) + c = ; 
and therefore, by subtracting F{x, y) = 0, 
a{^2xAx+Aai?)+2h{Ax.y+xAy+AixAy) 

+ 6(22/ Ay + Ay^) + 2grAa; + 2fAy = ; 
and dividing by Ax, 

a{2x+Ax)+2My+x^+Ay)+h{2y+Ay)^ 

+2^+2/^=0; 

which in the limit, when Aaj and Ay are indefinitely small, 
becomes the first derived equation. 
Again, from the implicit relation 

a?— 3aaJ2/ + 2/^ = 0, 

we obtain the first derived equation 

3aj2 - Zay - 3aic^ + S^/^^ = 0, 

so that dy_ x^-ay 

dx ax — y^ 
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Examples. — Find dy/dx from the implicit relations 

(1) a?+y^^c\ (4) aj« + 2/T = al 

(2) (x^aY+{y^hf = c\ (5) x'^y'^ = a'^+\ 

a^ b^ ' 

(7) Given the implicit relation 

a;^ — 3axy +2/^ = 0, 
prove that dy/dx = 0, when ic = a ^2, y = a^4 ; 
dx/dy = 0, when x = a^^, y = CL^^- 

(8) Find where dy/dx = and where dx/dy = 0, in the 

curve ic*+2/'' — 5tt2aj22/ = 0. 

14. Tracing Curves. 

The student should at this stage be exercised in tracing 
or plotting simple curves from their equations, thus 
exhibiting to the eye the flow or march of the function 
y = ix, or the graph of the function ix (§ 5). 

Most of these curves should be already familiar to 
the student who has studied the elements of Analytical 
Geometry, as given in Smith's Conic Sections. 

The curves should be drawn to scale as carefully as 
possible; for this purpose the paper which can be bought 
ruled into small squares is useful. 

By drawing these curves and interpreting geometrically 
the differential coefficients obtained from their equations, 
the student will begin to understand the use of the 
Calculus, always difficult to explain to the beginner. 

Examples. — Draw the curves : 

1 1 

(1) y = l, X, x^, x^, x^y fjxy ~, -^,....(on one diagram). 

X X 

1 

(2) y = x — x^, x—x^y x+a?y X'\ — , .... 

X 



DIFFERENTIATION. 25 

1 1 

(3) 2/^=1, Xy x\ a?, Qc^y fjx, -, -2, ....(on one diagram). 

X X 

(4) y^=\—o^, a?— 1, a?+l, x— a;^ aj+a^, o^—x, .... 

l+i=l l + i = l 

_ T^ -^J oT^ 9 — X,.... 

a; 3/ x^ y^ 

m _a:~l (a?-l)(a?-2) (a;-l)te-2) (a;-lXa?-3) ' 
^^y x-i aj-3 '(aj-3Xx-4y(aj-2)(a;-4y" 

(8) 03^— 4a3+2i/ = 0, a:^/ — 205— y = 0, {y — xY = \—x^, 

(9) Prove the properties of the following curves and 

sketch the curves : — 

(i.) In 2/2 = j9ic,aparabola, the subnormal MO is constant, 
(ii.) In 03^+2/^ = a^, a circle, the normal PO is constant, 
(iii.) In aj2+2/^=2a2/, a circle, OT^TP. 

(iv.) In x^+y^=a^y the part of the tangent TF inter- 
cepted by the axes is of constant length. 

(v.) In hy^ = a?y a semicubical parabola, 

TM^ = -i^h.MO'y and generally in 2/"' = ca?", 
MO varies as {TM)^^"'"^)l\ 

111 aj^+^ 

(vi.) In —-H — r = -;r, 02^=--, : and determine p, the 
\ / x** 2/ ^ a^ 

perpendicular from upon the tangent. 

(10) Prove that the equation of the tangent at {xy) of 

(i.) The drcle (fig. 13) x^+y^ = c^ is xx'+yy' = c^'y 
or of the circle (x — a)'^ + (y — by = c^ is 

(a>-a)(aj'-a) + (2/-6)(2/-&) = c^ 
(ii.) The eiZipse (fig. 8) ^'+^2 = 1 is 5-+f-'=^ ^ 
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{Solution : — Forming the first derived equation, 



25+2|?$/ = 0,sothat^=- 



• 

a? * h^ dx "' dx o?y ' 

and then the equation of the tangent TPV is 

ory a/' (r o? ¥ 

(iii) The hyperbola (fig. 14) g-g=l is 5'-^'=l. 

(iv.) The parabola (fig. 2 ii.) y^ = pxis y = 1- 

(v.) The hyperbola (fig. 2 iii.) a3y = c2 is „ — \-^ = l. 

(This is called the isothermal curve, as it exhibits the 
relation between x the volume and y the pressure of a 
given quantity of gas at a constant temperature.) 

(vi.) x'^y'^=a'^^''i& /i , ^/ x + /-. , ^ / x =1- 
^ ^ ^ {l+n/m)x {l+m/n)y 

(This is called the adiabatic or steam cui^e, when we put 
for instance m = 7 and -n. = 5, or -m = 10 and -n. = 9.) 

(vu.) (-) =(f) IS 73 7-T-+Ai — / \ =1- 

^ ^ \a/ \c>/ {l^n/m)x (l'-m/n)y 

Determine also the equation of the normal of these 
curves. 

(11) Prove that the curves aj^y'*=a^+** and - — ^ = 6^ 

cut at right angles. 

(12) Prove that the equation of the tangent to the para- 

x' v' 
bola v^aj+v^2/ = >v/a is —y — h-y- = >v/«, and then 

OT+OV=a; and generally, the equation of the 
tangent of x^ + y^=a'^ is x^-^x'+y^-^y'=a\ 
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(12) Prove that if 

(iii.)2^-2Z., 1 = ^(1+^;). 

(13) Prove that the conic section 

aQi?-\-2hxy+hy^+2gx+'lfy+c=0 
can be described with component velocities 
dx/dt = ?ix+by+f, dy/dt= — aa? — %— gf, 
or dx/dt= J{Qi^-ab)x^+2{fh'-hg)x+P'-hc), 
dy/dt= ^^{(h^-.ab)y^+2(gh^af)y+g^-cu^}. 
Prove also that the equation of the tangent at (xy) is 
(ax+hy+g)x+{hx+by+f)y+gx+fy+c=0. 

15. Algebraical and Transcendental Functions, 
So far we have dealt only with algebraical functions 
of Xy an algebraical function being defined as one which 
is composed of a finite number of the algebraical oper- 
ations of Addition; Subtraction, Multiplication, Division, 
Involution, and Evolution. 

An algebraical function of x may be defined in the 
most general manner as the root of an algebraical equa- 
tion of integral degree, the coefficients of which are 
rational integral algebraical functions of x; and a 
rational integral algebraical function of x is defined as 
the sum of a finite number of terms, multiples of integral 
powers of x, such as 

where n is a positive integer ; and then n, the highest 
power of X, is called the degree of the function. 
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But now we pass on to the functions of Trigonometry, 
such as sin Xy cos cc, ... , which are called transcendental 
functions; every function of x which is not an alge- 
braical function being called a transcendental function. 

Suppose we require to differentiate sin x ; then from 
the definition of a difierential coefficient in § 3, 

d sinic_ ,,sin(a;+A) — sinic__,,2 cos(x+ ^h)sm ^h 
dx k h 

i=]t cos (x+^h) —Y?—=Gosx; 

since It cos(a; + ^7i) = cos x, and lt(sin ^h)l^h = 1, 

when X and therefore h is expressed in circular measure. 

Or, more simply from the definition of § 5, 

d since ^,8in(x+h)—sm(x — h) ., sinfe 

— 7 — =lt — ^^ -Tj — ^^ ^ = ltcosa3— r — = cosaj: 

dx 'Ih h ' 

and so on for the other functions of Trigonometry. 

But if X is given in degrees, or minutes, or seconds, then 
dsina:**__ tt o cZsina;'_ tt , 

~rf^ 180 ^^^^ ' ~d^ 180 X 60 ^^^"^^ 

cZsincc"_ TT^ „ 

~d^ 180 X (50x60^^®^ ' 
so that extraneous factors are avoided by the use of 

circular measure. 

We thus require a knowledge of the circular measure 
of an angle and of the formulas of Trigonometry, a 
subject with which the student is presumably already 
familiar, but on which it will be desirable to make a 
slight digression, so far as is required for our purposes. 

16. Trigonometry. The Circular Measure^ and the 
Circular Functions of an Angle, 

Definitions, — The circular measure {Q) of an angle 
AOP is the number giving the ratio of the circular arc 
AP to the radius OA, 



DIFFERENTIATION, 29 

The unit angle is the angle subtended by an arc equal 
to its radius; this angle is called the radian^ and contains 

180 X 60 X 60-5-7r = 20626r-8, or 57' 17' W'S. 

Here tt denotes as usual the ratio of the circumference 
to the diameter, so that 27ra is the circumference of a 
circle of radius a; thus the cm. (circular measure) of 
a right angle is Jtt (radians), of the angle of an equi- 
lateral triangle is Jtt, and so on; and 7r = 3*14159265...., 
an incommensurable number. 

By Euclid VI. 83, the area of the sector OAP is to that 
of the whole circle as the arc -4 P is to the circumference, 
or as to 27r ; and the area of the circle being xa^, the 
area of the sector is ^a^O, or one half the rectangle of the 
radius and the arc. 




MAT 
Fig. 5 

In fig. 5 draw PM the perpendicular from P on OA, 
and let the tangent at P meet OA produced in T ; then 
according to the definitions of Trigonometry, 

arc^P . OM . MP , . 

QP = 0, yp=cos0, ^=sm0, 

PT . . OT .AM . 

^ = tan0, -^=sece, ^Q=vers0; 

also if TP produced meets OB drawn perpendicular to 

PV OV 

OA in F, Tjp = cot 0, ?^p = cosec 6 ; 

we shall call these functions of 6 the circular functions. 
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Then the triangle O^P = |a^ sin 0, the sector OAP 

= Ja^0, and the triangle OPr=|a^ tan 0; and these are 

obviously in ascending order of magnitude, so that 

sin0, 0, tan0, 

are also in ascending order, when Q is the cm. of an angle 

less than a right angle, that is when Q < Jtt. 

Now make Q indefinitely small; then since, when = 0, 

— ^ Itsin0/tan0 = ltcos0 = l; 

therefore also, lt(sin 0)/0 = 1, and lt(tan 0)/0 = 1 ; 

when = ; as assumed in the last article. 

More generally, when = 0, 

sin -3710 __ sin m0 ?>i_m tanm0_7n.^ 

nQ m0 n n' nQ n ' 

,, 1 7 /v sinfe® tan^° tt 

thus, when A = 0, —fjo—=—f^ = jgQ> 

sin^' tanfe' tt sin A" tan^" tt 



h' K 180x60 K' K' 180x60x60" 

Again, when = 0, 

,.vers0 ,. AM _,^ilf, chord ^P_^ 
lt__ _ It^^^^p-lt-jplt ^^^^p - 0, 

because It -r-rr = 1> and It-rvi = It sin i0 = 0. 

arc^P AP ^ 

As a numerical verification, the student may calculate 
from the Tables the values of (sin 0)/0, (tan 0)/0, and 
(vers 0)/0, when is the cm. of an angle of 3°, 2°, 1°, 30', 

10', 1', 30", 10", r. 

When = 7r/n, the cm. of 180/71 degrees, then 2ifP^ 
and 2Pr are sides of regular polygons of n sides, in- 
scribed and circumscribed to the circle ; so that 2na sin tt/n 
and 2'?i«rtairU/'n. are the perimeters of the inscribed and 
circumscribed polygons of n sides ; and the perimeter of 
the circle 27ra being intermediate in length, we can, by 
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sufficiently increasing n, determine closer and closer 
limits between which tt lies ; in this manner Archimedes 
found that tt lies between 3^ and 3^ J. 

As a numerical exercise the student may calculate limits 
between which ir lies, by taking ti = 180, 180x60, ... in 
n sin Q/n and n tan Qjn, 

17. Differentiation of the Circular Functions, 

Having thus justified the difierentiation of sin x, 

casino; 

— = — = cosaj, 

ax 
where x is the circular measure of a variable angle, we 
proceed similarly with the other functions : — 

d cos X __ , . cos (x+h) — cos a? __ , . — 2 sin(a; + |^)sin ^h 
cte ~ h ~ h 

= --It sm(x + ^h y -^ = — sinaj. 

<i tan aj _ , . tan(aj + A) — tan a? 
dx h 

= ItW— sec(a; + k)sec x = see^x, 

d cot X __ , , cot(aj + A) — cot x 
dx " h 

= — It— T— cosec(a;+A)coseca3= — cosec^a;. 

cJ sec a; _ , sec(a5 + ^) — sec a; _ , , cos x — cos(aj + h) 
dx " h ~" Acos(a;+A)cosaj' 

= sin xlco^^x = sec x tan x. 
Similarly 

dcosecaj . c? vers a? . _ 
= — cosecajcotaj, — = = smaj. 

dx dx 

The formulas of Trigonometry we have assumed in 
these differentiations are, with the usual notation, 
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sin(?-aini)= 2cos J(C'+7))ami((7-7>), 

cosC-co8i>=-2sinJ(C+i)}sini(a-i)), 

tanjl— tanS=3£n(jl — i})8ecjl sec-B. 

But these formulas are proved immediately by means 

of lig- 6, if we denote the angle AQP by C and j40Q 

byBjthen jlOii-lCC+jO), andJiO/' = J!OQ = J(C-Z)); 

™C-sinD.-^-J^-2£^=»£P ^ 

^^ oe oj> jipop 

= 2cosi(C+J)>mi(C-Z)); 

co.,C-eosC-°^-°^--2P^--2?'^ ^ 
co.,(, eosi'-yj, ^^ Zyp- 2^yy 

= -2sin J(C+ J))siu KC-H) j 

and denoting ^Oii by A and iiOP or ROI^ by £, 

tan^-tanit--^''-''e_er_*« «!■ Oe 

tan J ""-s-os os-oB-oglgoB 

=sinC4 — S)sec^ secB, 
18. Geometrical Differentiation of the Circular Fv/nc- 

tions. 




' MA 
Fig.8 

By means of fig. 6 we can illustrate geometrically 
the differentiation of the circular functions. 
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Thus if X is the cm. of A OR, and h of the small angles 
POR, QOR ; then employing the definition of § 5, 
dsina;_,.sin(aj+A) — sin(cc — ^)_,.iP — ifQ 
dx ~ 2h ~~ arcPQ 

arcPQ "PiJ- arcPQ ~*^'"- 
dcosx _, OL—OM 
dx arc PQ 

- I t 2pi? , _.^ chord PQ _ _ . 

With the definition of § 3, 
(Z tan aj_,. tan(a5+^) — tan aj_,.^u_ 
dx '^ h " SiXcrP 

tv arc tw arc rP 0-4 0^ 

tv, Oil i. ^v -1 /o i^\ J arc^ti; 0^ 

smce _=_— ,lt — T- = l (§ 16), and d = ?ti- 

^v OJ. arc^tt; arc rP OJ. 

c?seca3__, 2i;^i6 _,,tt;i6 m6 tv Rrctw 

dx " arc rP ~~ vu tv arc tw arc rP 

= tan X sec cc ; since 

./luu ^ ,,vu .,Au ^ 1. tv ^ axctw 

It — =1, lt7- = lt7c-T = tanaj, It j- = l, 5 = seccc. 

vu tv OA arc tw arc rP 

19. By means of the formula of § 11 for the difierenti- 
ation of a function of a function, we can differentiate any 
power of the circular functions of x, and the circular 
functions of any function of x ; thus 
d(sin x)^/dx = m(sin x)'^ ~ ^ cos x; dsin y/dx = cos ydy/dx ; 
d(sin yy^/dx = m(sin y)^ " ^ cos ydy/dx. 

As an exercise, the student may in the same manner 
difierentiate with respect to x the functions (cos a*)**, 
(tana?)**, (cotaj)^ (seca:;)^ (coseca?)^ (versaj)**, (f aj)^ sin fa;, 
{f (sin x)}^, (cos fa;)^. . . 
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Eosamples. — (I) Prove the following differentiations : 

(i.) 2/ = Jaj + J sin 2x, dy/dx = cos^o;. 

(ii.) y = Ja;— Jsin2aj, dy/dx = 8in^x. 

(iii.) y = Jcos^cc— cosaj, dy/dx = sin^x, 

(iv.) 2/ = sin a; — J sin^a;, dy/dx = cos^oj. 

(v.) y=^x+i sin 2x+^ sin 4a;, dy/dx = cos*a3, 

and write down y when dy/dx = siii^x, 

(vi.) y = tanaj— a;, dy/dx=taxi^x. 

( vii.) y = — cot aj — ic, dy/dx = cot^a;. 

(viii.) 2/ = J tan^oj — tan x+x, dy/dx = tan^oj. 

(ix.) 2/= J tan^aj+tanoj, dy/dx =^sec^x. 

(2) From the definition of § 3, differentiate with respect to a, 

sin 2a;, cos x/a, tan(ma;+7i), x tan x, tan x^, sin a;*" 

(3) Differentiate x^, sinaj, tana;**, with respect to x^. 

20. It will be useful also to draw the graphs of th< 
circular functions ; thus the equations 

y/b = ± sin 2'7rx/a and ± cos 27rx/a 
represent in fig. 7, for two values of 6, the plan of th< 
spirals of a double-threaded screw of pitch a. 
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Exam/plee. — Draw the curves: (1) y=8mx, cosic, tanx, 
cotiT, seciB, coseca;, versa; (on one diagram). 

(2) y = 3in"^a;, C03"^a;, tan "'a;, cot"^a;, sec"^, cosec"^, 
\ vers'^a; (on one diagram) ; and superpose them on 

\ the preceding diagram. 

(3) y = a;sina:, (3ina;)/3;, (aina;)^ sina^, ain^a;, arcot Jtrb. 

(4) 8ina;=siny, coaa3 = c08y, taQa;=tan^, 

(ain a;)*+ (sin y)^ = 1, (tan a;)^+ (tan y)^=l. 

(5) Prove that in the ellipse {xjaf + (y/bY = 1, we may put 

3;=acos9, i/ = fe8in9 (9 is then called the excentrio 
angle); and A0p = 6, and the equations of the 
tangent and normal are (fig. 8) 
(a;/ar)coa 6 + {y/b)am $—1; axsec6 — by cosec d=a^~ 6'. 




(Fig. 8 shows the dlipti^ trammel or comjpaeses, the 
theory of which is obvioua from this example.) 
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If z — Cy a constant, then V is the volume of a right 
cylinder on the base -4, so that 

F= cjjdocdyy and A ^ffdacdy \ 

and the area A may be considered as built up of the 
infinitesimal elements of area djxdyy the limits of integra- 
tion being taken so as to include all the elements of area 
in A, 

135. If the surface %—i{x, y) becomes a plane, then 

z=lx+'my+n, 
and Y^ff{lx + Tny + n)dxdy = (Ix + my + n)A = zA, 
where i is the ordinate of the plane standing on the 
centroid of the base A \ this is the formula for the volume 
cut off a cylinder by a slant plane and a cross section, 
or by two slant planes. 

By means of this principle we can calculate the volume 
of an earthwork dam across a valley, of which the contour 
lines are assumed to be parallel straight lines ; also the 
volume of a groin formed by the intersection of two equal 
barrel vaults, crossing at right angles. 

136. With polar (cylindrical) coordinates r, Q, and z, 
and double integration, the infinitesimal element of area 
enclosing a point (r, Q) in the plane oi z = may be sup- 
posed bounded by circles of radii r — \dr and r+ \dr, and 
vectors from making angles Q — ^dQ and Q+\dQ with 
Ox ; so that its area is 

\{{r+ldrf-{T'-^dvf}d6=dr,rd6\ 

and A—ffrdrdQ, V^ff%TdrdQ\ 

the limits of r and Q being taken so as to include all the 
elements in the area A. 
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In the general change from the mdependent yariabl< 
x.yio any two new variables u, v, the element of area 
bounded by the curves v,, v, u+du, v+dv ia ultimately a 
parallelogram, the coordinates of whose angular points 
are, relatively to the comer (cc, y), 

^^J? ^J ^J . ^Jl ^.Vj .^7 '^J ^J 

dw du ' du dv du dv ' dv dv 
and therefore the area of this elementary parallelogram is 

\du dv dv du/ d{u, v) 

so that now, with z expressed as a function of u, v, 

JJ 3(u, V) JJ 3(u, V) 

Thus, for instance, in Thermodynamics, the pressure 'p 
and volume v of imit mass of a perfect gas are connected 
with 0, the absolute temperature, and ^, the entropy, by 
the relations (Maxwell, Heaf) 

and, taking 0, ^ as independent variables, 

~pd0 vdO & p'd4> vd4> ' 

ndp^mdv^^^ ndp.7ridv_^^ 

pdd vdd~ ' pd^ vd^'' ' 

sothat, !?^::!^^= m m^3p^_l 

p 30 p d4> 

m—n dv n m—n dv 



V de & V 30 



= 1 



and ?tej?^ = _l_2?!=l. 

3(0, 0) m-ti c0 

Therefore the area on the p, v diagram of the Camot 

cycle bounded by the two isothermals O^, 62 and by the 

two adiahatics tp^, ^g is (^i""02)(0i""^2)- 

T 
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137. The Sign to be attributed to an Area. 

Suppose a point to travel once round the closed ov 
area A, an indicator diagram for instance, so as always 
have the interior of the curve on the left hand (fig. 44). 




Fig. 44 



Then A = //docdy = /xdy -jx -^U 

taken round the perimeter of the curve. 
From Bio C along BPG, dy/dt is positive, 



'x^^dt= QxefiMBPGN) 






and ix 

and from G to B along CQB, dy/dt is negative, 

and l^f^^ ^ "" ^^®^ MBQCN ; 

so that, taken round the curve, 

/x-^dty 01 /xdy=^ Ay the area of the closed curve. 

But JJ^y^^Jy^'yy'dt^^' 
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and from EtoD along EPD, dx/dt is negative, so that 



fy 



y^.dt= -area LEPDK; 



A 



dt 
and from D to E along DQE, dx/dt is positive, so that 

5§d^= ei.ve&LEQDK; 

and therefore, taken round the curve, 

Jy-^^> o^Jy^ = - ^ ; 

Therefore, taken round the curve, 

J Vd^'^^r^* oxj{ydx+xdy) = 0\ 

and ydx+ocdy = d(xy) 

is called a perfect differential; its integral between two 
limits is independent of the intermediate values of x and 
y and of the path described between the limits; so that, 
taken round any closed path, the integral is zero. 

I 

^d xdy—ydx is not a perfect differential, so that its 
Vtegrated value depends on the path taken between two 
.mits ; and in a closed path it represents twice the area 
inclosed by the path. 
Similar theorems hold for the integrals 

JJzdxdy dindi JJzdydx, 

138. Changing to polar coordinates by putting 

a5=rcos0, 2/ = ^sin0, 

then «,J_2,_-=^_, 

80 that A =f¥'^^^ =/Jr«c?d, 

taken round the curve. 
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For an origin inside a closed oval curve, the limits of 
Q may be taken as and 2x. 

But if the origin is outside the area, draw the 
tangents OF^ OH to the curve from ; then along FPHy 
dO/dt is positive, and 



A 



i7^%It=: Bxes^OFPGHO; 



dt 
but along ERF, dO/dt is negative, and 



A 



ir^^dt= -area OFQHO ; 



so (hat, taken round the curve, 

A =J\i^^'t =Ji't^de =JixHv, if 2//a; = w = tan 0. 

Although not a perfect differential, xdy—ydx can be 
made one by dividing by a?+y\ and then 

so that f^^^^iy^=fde^2^. or 0, 

according as the origin is inside or outside the curve. 

Consequently if each point of the contour is displaced 
through small distances c^xjr^, c^yl^^, parallel to the axes, 
the change in the area will be 2x0^ or zero, according as 
the origin is inside or outside the contour. 

139. A convenient independent variable to take is 8, 
the length of the arc of the perimeter measured from a fixed 
point; so that the point P may be supposed to move 
with unit velocity ; and now, with p denoting the length 
of the perpendicular from on the tangent at P (§ 9), 

dy dx 
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so that A =/\jpd8=^/^p^dr 



~/hP s®c 0^^ =/jl^^^/V(^ ""?*)• 



When fig. 44 represents an indicator diagram, and KL 
the reduced stroke of the piston, while the ordinate y 
represents the pressure of the steam, the pencil will de- 
scribe the contour with the area to the left, when the 
steam pressure is urging the piston from £ to i^. 

The diagram taken on the return stroke from the other 
end of the cylinder will be described in the opposite sense, 
with the area on the right hand of the describing pencil. 

Sometimes a loop is found on the diagram, described 
in the reverse sense ; this loop shows that a cushioning 
efiect takes place, which requires attention, as the nega- 
tive area of this loop represents so much lost work. 

In the general case, where the perimeter cuts itself a 
number of times, then the area obtained by integrating 
once round a loop will be positive or negative with the 
above formulas, according as the area is on the left or 
right of the describing point, as it travels round the curve. 

(Clifibrd, GoTnmon Sense of the Exact ScicTices ; 
Cremona, Oraphical Statics.) 

Familiar instances of such looped contours are seen in 
Lissajous's figures (Ganot, Physics, § 281), whose general 
equation may be written as (§ 103) 

msin"^/a— 7isin"^/6=a constant, 
or x/a = sin (71^ + e), y/b = sin(m^ + eO- 

Also with the polar curves 

r=acos(m0/'M), 
where m and n are integers, curves seen on the back of 
an engine-turned watch. 
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140. The Planimeter. 

This instrument in the most usual form, that invented 
by Amsler of Schaffhausen, consists of two bars OA^ AP, 
pivoted at and jointed at A, and carrying in PA 
produced a small graduated roller R, with axis fixed 
parallel to PA (fig. 45). 




Fig- 45 



The instrument is used to measure areas ; the pointer 
P is carried round the perimeter of the curve whose area 
is required, an indicator dii^am, or the section of a ship 
for instance ; the roller M, which rolls and slides on the 
plane of the paper, then registers the area. 

If precision is required, the point P may be carried 
say ten times round the perimeter, and the reading of the 
roller be divided by ten. 
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To explain the theory of the instrument, we shall 
suppose the pointer P to travel round a finite area 
PP^P^P^ (something like a Carnot cycle in Thermo- 
dynamics), in which PPi, P^z ^^ circular arcs described 
round as centre, and PxP^^y ^^Pt are arcs round A^, A 
as centres ; in this way we analyse t^e motion due to 
the joints and A separately. 

Let OA = a, AP=b, AR=c; and let the direction of 
a positive motion of the roller, as marked by the gradua- 
tions, be that which on a right-handed screw would give 
a motion in the direction AR, 

Drop the perpendicular 01 from on AR. 

(i.) Fix the joint -4, and move P to P^ by rotation 
round on the circle PP^ through an angle 0; the 
angular velocity dQjdt will give to the roller R the 
component velocities 01, dO/dt in the direction IR and 
I R. dQjdt perpendicular to IR\ the first component 
drags the roller over the paper, and the second com- 
ponent makes the roller turn with circumferential velocity 
IR . dQjdt ; and therefore the whole travel of the roller, 
or its graduations, will be IR . ft 

(ii.) Fix the joint 0, and move Pj to Pg by rotation 
round the joint A^ through an angle <f>\ the angular 
velocity d(f>jdt of AP will communicate a circumferential 
velocity — cd(pjdt to the roller; and the travel of the 
roller will therefore be c^, backwards. 

(iii.) Fix the joint A, and move P^ to P^ by rotation 
round through an angle Q ; the roller will now move 
backwards, and its travel will be I^R . Q. 

(iv.) Fix the joint 0, and move Pg to P by rotation 
round A through an angle ^; the travel of the roller 
will be C(f>, forwards ; this cancels the travel of (ii.). 
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In completing the finite circuit PP^P^^, the total 
forward travel of the roller will then be (IR'-I^R)6, 
But the area PP^P^^ = area PP^QjQ 

= sector OPPi- sector OQQi 

=.\{OP^^OP^)e 

= l{OA^ + AP^-¥2AI.AP'-OA^'-AP^''2AI^.AP)e 

= b{AI^AI^)d-=b{IR'-I^E)e 

= 6 times the travel of the roller. 

Thus by altering the length of b by an adjustment of 
the instrument, which allows the arm AP to slide in the 
sleeve -422 and be clamped, the area can be read off 
in any required unit, say the square inch or square 
centimetre. 

Any irregular area, such as for instance an indicator 
diagram or the cross section of a ship, must be supposed 
built up of infinitesimal elements formed in the same 
manner as PP^PJP^; and will be read off when the pointer 
P completes a circuit of the perimeter, both joints being 
now free to turn simultaneously. 

When / coincides with i2, the roller will not turn, and 
then P describes a circle called the zero circle, represented 
-by the middle dotted circular line (fig. 45) of radius 
J(0B'+RF')=^^{a^^c^+{b+cy} = ^(a^+b^+2bc). 

When / lies on the same side of -B as -4, the travel of 
the roller is reversed, but in a complete circuit the reading 
is unaltered. 

If however the origin is taken inside the area to be 
measured, the area of the zero circle must be added to the 
reading of the roller. 

Prof. G. B. Mathews explains the theory of the Plani- 
meter in a slightly different manner, by 
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(i.) moving OA to 0-4p and P to P^, keeping AP 
parallel to itself; and then the curvilinear area APP-^A^ = h 
times the travel of the roller ; 

(ii.) moving AJP^ to A^P^ by rotation round A^ ; 
(iii.) moving OA^ to OA and Pg to Pg, keeping -^.^Pg 
parallel to itself; and then the curvilinear 2k£e2i,AA^P^^^h 
times the travel of the roller, backwards ; 

(iv.) completing the circuit by moving AP^ to AP by 
rotation round A, when the sector -4PP3= sector A^PJP^;, 
and the travel of the roller cancels the travel in step (ii.). 
Therefore the 

area PPiP2P8=area -APP^ili— area A-^P^^A 

= 6 times the travel of the roller. 
The end A may be guided in a slot of any form and 
the area will be read off as before ; a straight slot is often 
employed, with the advantage that the pencil P can then 
cover a greater area ; and with appropriate mechanism 
can be made to register the moment of the area, and its 
moment of inertia about the straight line of the slot. 

(J. F. Bramwell, British Aaaodation, 1872 ; 
H. S. Hele Shaw, Proc. I.CE,, 1885.) 

Exam/pies, 

(1) Prove that the area is 7r{Bc—bC) of the curve, an 

ellipse (Ex. 4, p. 203), given by 

a3=a + 6cosd + csin 6, y=A +J5cos0 + (7sin ft 

(2) Trace and find the area of the curve 

(r—a cos Of = a^cos 20. 

(3) Prove that the volume cut off from the surface 

2^=Aa^+2Hxy+By^ 
by the plane z=c ib n/(n+l) of the cylinder on 
the same base. 
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141. Functions of three or more Independent Variablea. 

A function of three independent variables, x, y, z, de- 
noted by {{x, y, z), may be supposed to represent some 
function of the position of a point in space whose co- 
ordinates are x,y,z; for instance, the density or tempera- 
ture or pressure at the point. 

Then {{x, y, z) = G, bl constant, would imply a relation 
connecting x, y, z, and would be the equation of a surface; 
for instance, a surface of equal density, or pressure. 

If F denotes the volume contained in a closed sur&ce 8, 

then Y —fffd/xdydZy 

the integration including all the infinitesimal brick shaped 
elements of volume dxdydz (fig. 43) contained in S. 

When the density p within the surface 8 is variable 
and a given function of x, y, z, then the mass ilf contained 
by the surface 8 is given by the triple integration 

M=//fpdxdydz. so that ^3^=^. 

and the mass is the space integral of the density p 
throughout the volume F"; while x, y, i, the coordinates of 
the centre of mass, and kx, kyy k^^ the radii of gyration 
about Ox, Oy, Oz, are given by 

^^ =///oopdocdydZy yM ^^fffypdocdydz, 

zM=^fffzpdxdydz ; 

Mkx^=///{y^+z^)pdxdydz, Mky'=///{z^+x')pdxdydz, 

Mk^ -=-///{x^ + y^)pdxdydz. 

Functions of more than three independent variables 
cannot be interpreted geometrically without the introduc- 
tion of the fiction of space of more than three dimensions, 
a thing which is inconceivable. 
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A function F(^, x, y, z) of four independent variables 
t, X, y, z, may however be interpreted, as in Hydrodynamics, 
as representing the velocity, or density, or pressure, at the 
time ^ at a point in space whose coordinates are x, y, z, 

142. We have used the words mass and density; the 
Tnasa of a body is the quantity measured by the balance 
against certain standard lumps of metal, called weights in 
the Acts of Parliament (French, poids, German, Gevdchte), 
the standard in this country being the Pound Weight, and 
in the Metric System the Kilogramme of 1000 grammes. 

The density of a body is defined as the number of 
units of mass in the unit of volume ; with British imits, 
the density is the number of lbs. in a cubic foot of the 
substance, and with Metric units is the number of 
grammes in a cubic centimetre, or of tonnes of 1000 kilo- 
grammes per cubic metre. 

The units of length are thus the foot in the British 
System, and the Tnetre or centimetre in the Ifetric 
System; while the unit of time in universal use for 
theoretical investigations is the secondj the mean solar 
sexagesimal second. 

When therefore we speak of a time t, we mean t seconds; 
and coordinates x, y, z are measured in the unit of length, 
which is either the foot or else the tnetre or centimetre. 

For practical purposes there are only three systems of 
fundamental units which need be considered, 

(i.) the British foot-pound-second (f.p.s.) system ; 

(ii.) the C.G.s. (centimetre-gramme-second) system ; 

(iii.) the metre-kilogramme-second (M.K.S.) system; 
and from these fundamental imits of length, mass, and time, 
all other units, of area, volvmie, density, velocity, accelera- 
tion, momentum, energy, force, etc., may be derived. 
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*143. Spherical Polar Coordmatea. 

In this system of coordinates the position of a point P 
on a sphere with centre at is defined by 6, its angular 
distance from a fixed pole JN on the sphere, and by <f> the 
angle which the plane ONP makes with a fixed initial 
plane ; so that on the terrestrial sphere, ^ will be the 
longitude and Jtt— the north latitude, if JV denotes the 
north pole. 

By taking r the radius of the sphere as variable, we 
can define the position of any point in space by means 
of the three quantities r, d, ^, called the spherical polar 
coordinates in space. 

With ON coincident with Oz, and with the plane xOz 
as the prime meridian, these coordinates are connected 
with the orthogonal coordinates oj, y, z, by means of the 
relations x=rsin6 cos <p,y=r sin 6 sin <p, z=r cos 6. 

If we replace rsind by tJJ, then ts, 0, z, are the co- 
ordinates in the cylindrical system (§ 132). 

The element of volume cut out by the spheres r-^^dr 
and r + ^dr, by the cones 6—^d6 and 6 + ^dd, and by the 
planes ^--^dip and <f> + id(f>, will be ultimately 

dr . rdO . r sin 6d(f> ; 

so that by triple integration 

r=///rhm edrded(t>y and ^^ = r^sin 6 ; 

while the mass M and its centre of mass, for variable 
density />, are given by 

M =///'pr^sm edrdOd^, 
xM =:J^pr^ain^6 cos ^rd6d(f>y 
yM=:y//'p7^am^e sin ^drded<f>, 
zM ^//jfp^^sin e cos OdrdOdip. 
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*144. Space, Surface, and Lime Integrals. 

Consider a fixed closed surface S, and a function X of 
the coordinates x, ^, of a point in space. 

Then in the triple integration extending throughout the 
volume enclosed by the surface S, called a apace integral, 

where Z^ Xg, Xg, . . . denote the values of X where a point 
moving from ~ 00 to 00 parallel to the axis Ox succes- 
sively enters and leaves the interior of the surface S. 

Denoting by l^, l^ {3, ... the cosines of the angles which 
the outward drawn normals of the surface S at these 
points make with Ox, then 

dydz^ —IjdSi^l^dS^^ ^l^dS^==... , 
supposing the infinitesimal prism on the base dydz parallel 
to Ox to cut out the elements of surface dS^, (i/Sfg, dS^,..., 
on entering and leaving the surface S. 

Therefore, denoting the element of volume by dV, 

the double integration extending over the surface S, and 
this is called a surface integral; so that a volume integral 
can always be expressed as a surface integral 

Similarly, with T, Z other given functions of x, y, z, 

where m, n denote the cosines of the angles the outward 
drawn normal of the surface S makes with Oy, Oz. 
Therefore, by addition 

^(l+f+S)'^^=#(^^+-^+-^ (^-^ 

Thus V=/J/dxdyz = J //'(Ix +my + nz)dS. 
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For instance, integrating over the surface of the ellip- 
soid, p denoting the perpendicular from the centre upon 
the tangent plane, 

fjpdS = ^irahc, SMdydS/p = ^Ttibc/a + ca/b + ab/c) ; 

and, as an exercise, the student may calculate JclS/p and 
/d8j^ for the hyperboloids, as well as for the ellipsoid. 

Again, suppose X, F, Z are the component forces per 
unit of volume acting throughout a fluid at rest, in which 
the pressure at any point is represented by p ; then the 
equation of equilibrium of the fluid within the closed 
surface 8, obtained by resolving parallel to Ooj, is 

///XdV^f/lpdS. 

But changing the surface integral into a space integral, 



ff^'jift^y. 



so that Tr- = ^\ and similarly ;rS= F, -^^Z\ 

or dp=Xdx+Ydy+Zdz, 

so that the space va7%ation of the pressure of a fluid at 
rest in any direction is equal to the component force per 
unit of volume in that direction ; and surfaces of equal 
pressure are cut Orthogonally by the lines of force. 

As another illustration, we may suppose X, Y, Z to 
represent the components of flux (estimated with British 
units in lb. per square foot per second) of a fluid in 
motion ; then ff(^X + m F+ n Z)dS represents the number 
of lb. which is flowing out across the surface 8 per second; 
while if p denotes the density, in lb. per cubic foot, of the 
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fluid at any point of the interior of S, so that the mass M 
within 8 \& fffpdiVVo,^ then 'dMj'dt represents the rate of 
increase, in lb. per second, of the quantity of fluid inside S. 
Equating this gain and loss, 

ZMjU +//{IX+ mY+ nZ)dS = ; 
and, replacing by space integrals, 

leading to the equation of continuity in Hydrodynamics, 
when we replace X, F, Z by pu, pv, pw, where u, v, w 
denote the components of velocity of the fluid. 

In a plane, dA denoting an element of the ai'ea A, of 
which ds represents an element of length of the closed 
contour s, equation (ii.) becomes 

#(S+f>^=>^+-^>^ • (-)' 

thus expressing a surface integral by a Zme integral; 

and I = dyjds, m = — dxida ; 

or, if F(aj, 2/) = is the equation of the contour, 

A similar theorem connects the surface integral on a 
curved surface /S, which is a portion of a closed surface, 
with a line integral round the edge of 8 ; 

#{C-f)^»'(i-D-''(f-f)h 

-A^^^^O <"■) 

(Maxwell, Electricity, i., p. 25), this surface integral 
vanishing, by equation (ii.), for a closed surface. 
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For a plane surface, we may take Z=0, m=0, «=l; 
and now, as in (iii.), 

*145. QreerCs Theorem, 

Now suppose U and U' ^e given functions of Xyy,z\ 
then from equation (i.), integrating by parts, 

and therefore, denoting by — y^ the operator 



#■( 



3E7 3U' . -dUW . 3Cr3 



_^3D'3?7^ , 3tr3?A^y^ 



3aj 3fl5 dy dy dz dz / 
and therefore, by symmetry, 

^^^f'^^+j^^^^^''^^' (-^•)' 

where &+mP^+n^^=^^, 

dx dy dz dv 

representing the rate of growth of U in the direction of 
the outward drawn normal of S. 

Equations (vi.) and (vii.) constitute Oreen's Theorem, a 
theorem of great use in the mathematical theories of 
Electricity and Magnetism. 

(An Essay on Electricity and Magnetism, by G. Green ; 
edited by N. M. Ferrers.) 
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*146* Change of the Variables in Space Integrals, 
Generally in changing from x, y, z to any new inde- 
pendent variables u, v, w, we may consider that space is 
divided up into elements of volume bounded by the 
surfaces for which u, v, w are constant; and now the 
element of volume will be changed' from doodydz to 



d{u, V, w) 



J\ ' y* — I denoting the determinant 

^(t6, V, w) 



dx dy dz 
du du dw 
dx dy dz 
dv' dv' dv 
dx dy dz 



dw' div dw 

called the Jacobian of x, y, z with respect to u, v, w. 

Denoting by ds, the element of length, then df? becomes 
changed from dx^ + dy^ + dz^ to 

fdXj dXj dx. V (dy. \^,{'dZj V 

= AHu"^ + Wdi^ + G^dw^ + 2Ddvdw + 2Edwdu + 2Fdudv, 

J2_^^ 32/2 dz^ • 7)_'9^ ^ ,'^y ^2/ .'^^ '^^ 

du^ dv?' dv?^"'* dv dw dv dw dvdw*"" 

It is convenient to choose the new variables u, Vy w, so 
that the corresponding surfaces cut at right angles, and 
then D, E, F vanish ; while 

1 dx 1 dy 1 dz 
A du A du A du 

represent the cosines of the angles the normal to the 

surface u makes with the axes Oar, Oy, Oz. 

u 
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But with Xf y, z as independeDt variables, and denoting 

, 1 du 1 du 1 du 

h^ dx* h^ 'dy' h^ dz 
are also tba cosines of the angles which the normal to the 
surface u makes with Ox, Oy, Oz; so that 

1 du__l 'dx 1 du__ 1 3t/ 1 dii^^l dz 
F^^^A du* F^^ " A du' h^^'~'Adu* 
and similarly, 

1^=-^ 1—=:L^ A5?^=^?^- 
K^dx B dv' h^ dy B dv' h^ dz G dv' 

Idw^ldx j^5^^1% 1^ = 2??- 
Ti^dx Udw' F^dy C du/ h^ dz Bdw' 

A 1 • 7 ^^ 7 . ^^ 7 . ^W^ 7 

Also, Since a «t = ;r- ace + ;^ at/ + ^r- a^;, 

dx dy ^ dz 

. dw ^ , dw J , dw J 
div=^7^dx+-:r-dy+^r-dz; 

dx dy ^ dz 

therefore, multiplying these equations by 

l du 1 dv 1 dw 

h}dx'h}^' h^dx' 

, ,,. , 1 3u , . 1 dv J . 1 dWj 
and adding, <fo;=^, ^dn+^i ^"+ V ^'^^''5 

so that 7r- = i-o ;r", ... ; or -4A,=fi/t2 = CAq=l. 

du h^ dx ^ ^ ^ 

Now denoting bj^ c^j, cfeg* ^3. ^^^ elements of the 
normals, intersections of the surfaces u, v, w, 

ds^ = Adu, da^ = Bdv, ds^ = Cdw ; 

30 the.t^^== ABO, and ^^=h,h^=l/ABC 
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Now if V denotes any function of x, y, z, or u, v, w, then 

also, 

which are readily proved by taking the axes parallel to the 
normals to u, v, w at the point. 

With spherical polar coordinates r, 6, <p (§ 144), 
c?82 = dx^+ dy^ + dz^ = d^^ + rH0^ + r^sin ed<f>\ 
so that A = l, B = r, (7=rsin 6; 

^ ^"r^diK dr J ' r^sin 6 ddV^^ dO Jr^am^O ^02• 

■*147. Confocal Quadrics. 

A familiar instance occurs with confocal quadric sur- 
faces, where, with X, /x, v for new variables, 

1^70 , ^ 1^ 2iA ' 2i i^'" — ■■•> 9~; — I • • • — J- ) 



a?+\ 6HX c^+X aHjtt a^^,, 

whence, by solution, 

and A" (X-M)(X-y) ^^ ^2_ 

^"'^ ^ ~4(a2+X)(6''+X)(c2+X)' ■••'• 

Denoting by p the perpendicular from the centre on 
the tangent plane of the surface X, then we can show that 

p" {a? +Xf {h^~+ \f {c'+Xf ^x a* + X' • • ■ ' ^ ~ *^ ' 

, s^x . a^x , 3^ „ 2/ 1 ^ 1 , IN 

and thence that F=fX will satisfy the condition v*r=0, 
provided f X = C/{ {a? + X) (6* + X) (c^ + X) } - ^d\. 
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Examples, 

(1) Interpret geometrically the equations 

fr=O,f0 = O,f0 = O,f(0,0) = O,F(r,0) = O,F(r,0) = O. 

(2) Prove that the quadric surfaces of revolution, 

'^'■'^%.-i=^ (i.). 



sec^i6 tan^u 
se'ch^ tanh^i; 

X^ /.tZ 



= t^ (ii.), 



2 ^2- =0 (iii), 

cos-^tc; sin^ty ^ ^ 

are (i.) oblate spheroids, (ii.) hyperboloids of one 

sheet, (iii.) planes, forming a system of confocal 

orthogonal surfaces ; and sketch the figure. 

Prove that they intersect in the point 

x = c sec u sech v cos w, y = c sec u sech v sin w, 

and that the generating lines of (ii.) are given by 
cc-csecusechvcosty __ 2/'"^secttsecht;sint(; __ 2;-ctanutanhv 
sech'i;sin(u±t(;) ""-sech'ycos(u±t(;) ~" Itanh-y ' 
the angle between them being 2cos"^(cost6sechi;). 
Write down the corresponding equations for a 
system of confocal prolate spheroids and hyper- 
boloids of two sheets. 
(3) Prove that the equations 

y^ z^ 
— ?^i— + - 1.21 =8a(acoshu — a?) (i.), 

— ^1 ■ 91 =8a(acosh'y — a;) (ii.), 

cos^Jv sinr^v ^ / \ /» 

y^ z^ 

. fo| — I 1:01— = ^cl(<^ cosh ty +aj) (iii.) 

sinh^^ty cosh^^ty ' ^ ^ ^ 

represents a system of confocal paraboloids inter- 
secting at right angles in the point 
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X = a(cosh u + cos v — cosh w)y 
y = ^a cosh |u cos \v sinh Jic;, = 4a sinh |t6 sin ^v cosh \w. 
Determine A^ B, Cfor this system, and the equa- 
tions of the generating lines of (ii). 

Prove also that v^('^> t;, or ty) = ; as also in Ex. 2. 

(4) Verify that 1/r, 0, log tan |0, and their products are 

annihilated by the operator v^. 
Also r'^cos 0, r"%in Q cos(0 + a). 

(5) Prove that, if A and B are the ends of a diameter of a 

sphere of radius a whose centre is at 0, the function 

y_ 1 _J 1_, AN+AP 

AP BP AB ^^ BN+BP' 
where N is the foot of the perpendicular from P 
on AB, is annihilated by v^J ^^^ *lso that 
dV/dr = 0, when r=a. 

Give the physical interpretation of this result. 

(6) Prove that 

^1^ ^2^ ^3^ "* ^n^ 

" lfl3r\ drJ^W^Ku^ 30i/"^"*"*"30n-iW-i B0n-i/J 
where 

ojj = r sin Q^sin Og- • •^^'^ ^n-i. ^i = ^'^ 

ajg = r sin O^sin 6^. . .cos On-i, ^^2, = ^^sio^^^, 

ajg = r sin ©jsin 02- • -cos 0n_2, 1^3 = r^in^^^sin^^g, 

a;n-i = ''8in O^cos 0^, Un.\==r%m^Q^9\n^Q^,,.^\v?0ri-ii 
Xn = rcos6i, Tr=r'*-^sin~-20iSin~-302"*sin0n-8. 

Examine the case of 71 = 3. 
Prove also that 
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148. Quantws, 

A rational integral homogeneoua algebraical function 
of the 71*^ degree in m variables x, y, z, ... is defined to 
be ** a function in which the sum of the indices of the 
variables in each term is constant and equal to n, the 
indices being positive integers'' ; such a function is called 
a quantiCf and is denoted by (x, y, 0, ...)^ 

Thus (x, yY, or (a, 6, c, ,,.){x, yf 

represents the binary quantic in oj, y of the ti*^ degree 

ax''+nbx^-'^y + \n{n — l)cx^-^y^+.,.. 

Denoting the general quantic by u, then 

This is proved for each single term of the quantic, say 

for, denoting this term by v^ then 

'dv dv dv 

, 'dv ^ 'dv 

and x^ — \ry^^+Z;r-+ ...=nv, 

ox ^oy oz 

More generally 

and so on; in which it is important to notice that the 
expression in brackets on the left hand side must be ex- 

panded by the Multinomial Theorem as if a?, ;r-, y, — .,,, 

ox 'dy 

were independent algebraical quantities, and then u sup- 
plied in the partial derivatives. 



\ 
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For if \X' — hy- — \',..ju meant that the operation 
x—+y—+,., was repeated k times, we should have 

The operator is therefore a particular ease of the operator 

cc':^+ 2/:^+..., denoted by A, and called by Klein the 

polarizing operator. 

These theorems are proved by expanding the quantic 
u^ = {x+lix, y+hy, z+hz, ...)^=(1+A)**u 
in the form (§ 1 26) 

^ \ bx ^dy dz ) 

and equating coefficients of like powera of A. 

These theorems are called Euler^s Theorems of Homo- 
geneous Functions, or Quantics. 

A non-homogeneous function can be always made to 
appear homogeneous by the introduction of an appropriate 
factor to each term, consisting of the requisite power of 
some new variable, which may afterwards be interpreted 
as unity. 

Thus the general equation of a conic section 
ax^+2hxy+by^+2gx+2fy+c=0, 
can be made homogeneous, and a ternary quadratic in 
^yVj ^> hy writing it 

a^+by^+cz^+2fyz+2gzx+2hxy = ; 
and now x, y, z may be considered as trilinear coordinates; 
and replacing z by unity gives the ordinary Cartesian 
coordinates, with oblique axes. 
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*14j9. Arbogast'a Method of Derivation. 
Denoting f(x+x% y+y% z+z't, ...) by f, then 

so that the polarizations of f(a;, y, z,...) by the successive 
operations of A are the coefficients of t, t^/2\, ^/3!,... in the 
expansion of f in powers of t 
Again, if we denote 

aQ+a^x+a2-^^+ ... +ar^+ ...byjf, 

and thence generally 

Therefore, expanded in powers of a? by Maclaurin's Theorem, 

80 that the coefficients of the powers of x are derived by 
the repeated operation of 

this process is called Arbogast's Method of Derivation. 
In Arbogast's second method of Derivation 

fy = {(a+bx+cx^+dx^+ex^+...), 
is expanded in powers of x, in the form 

^f^^=A+Bx+Gx^+Da^+Ex^+ ..., 
>^^^^kand now we find 
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B={'a.b 
0=fa.c+r'a.^ 

Z)=ra.d+f"a,6c+f"'a . ^, 

E=f'a . e+ra(ic?+ic«)+f"'a . i6«c+f""a. -^ 

and 80 on, by a simpler mode of Derivation. 

For example, applying this method to Ex. 3, p. 234, 

here a = l, 6= — |, c = J, rf= — ^, 6 = |, ... , 

and fa=ra =f"a= = e ; 

thence ^ = 1, 5= -^ G=ih ^= - i^, ^=fli*, ... • 

*150. Theorems of Lagrange, Laplace, and Burmann, 

Put y = a+bx+cx^ + doc?+...y 

and (iiy = b + cx+dx^+ex^ + . .. ; 

so that y = a+ x<f>y, 

Tien |.,,+^-,|,|.l+^V|. 

and thus ^=z(hy-^. 

dx ^^da 

Then if u = fi/, 



and since 

3 

therefore 

3a;2 






■I(*^^l)4(*^'l)=^{(*^)'£}^ 
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and generally, by Induction, 

Now make aj = 0, therefore y = a, and u = fa; and 
u = f2/ = fa+5^:^,{(^a)«f'a}; • 

this is called Laaranffe's Theorem , 

'hus, for instance, putting iy = <py, we find 

Again, suppose y = ix + \My^ - 1 ) ; 

then ^='^+2^!^V)- 

But, on the condition that y = fx when /i = 0, 

so that (ex. 6, p. 276) 

" 7i! d/xA 2 / • 
Laplace has given a slightly more general form to the 
theorem of Lagrange by putting 

z = ¥y^Y{a+x<j>y)\ 

and now ;r- = d>y7c-y as before : and 

iz = f (Fa) + 2 J ^J{0(Fa) }»f '(Fa)]. 
Next put X or (y — a)/(f>y = \[/y ; then 

called BurmanrCs Theorem, giving the expansion of fy in 
powers of any other function \fry, y = a being a root of 



CHAPTER VI. 

CURVES IN GENERAL. 

151. Curve Tracing in Cartesian Coordinates, 
A number of such curves have already been introduced 
previously, which presented no difficulty in tracing from 
their equations; a slight sketch will now be given of a 
systematic method of treatment, but for a more complete 
account the student is referred to the treatises on Curve 
Tracing by Frost or Woolsey Johnson. 

Given the equation of a curve in the rational integral 
algebraical form of the implicit relation f(aj, ^) = 0, first 
group the terms in binary qnantics (§ 148), arranged in 
descending degree 71, ii — 1, ... , 3, 2, 1, ; and denote 
them by Um Un-u ... W3, itg, Up Uq] so that the equation 
of the curve becomes 

Un + Un-1+ ... +n^ + U2 + u^ + UQ=0 ; (A) 

then n is called the degree of the curve. 

If Uq does not vanish, the curve does not pass through 
the origin; but by changing the origin to a point on the 
curve we can make Uq vanish, and now the new u^, 
equated to zero, gives the tangent at this point. 

Also ui2+u^ = is the equation of a conic section, 

osculating or having a contact of the second order (§ 118) 

at the point. 

315 
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With rectangular axes, (7(aj2+t/^)+Ui = will be the 
equation of a circle touching the curve ; and if 
u^=Ax+By, u^^aa^+2hxy+by\ 
then this circle will be the circle of curvature, if 

Similarly u^+U2+Uj^ = will represent a cvinc curve, 
having a contact of the third order (§ 118), and so on. 

If i/j also disappears, then ^2 = represents two straight 
lines, real or imaginary, through the origin, which are 
tangents at this point ; if they are real, the origin is a 
double point, and the curve crosses itself; if imaginary, 
the origin is a conjugate point, that is, an isolated point, 
the coordinates of which satisfy the equation of the curve. 

If 1^2 is a perfect square, the tangents are coincident 
and the origin is in general a cusp (§ 104). 

If u^ also vanishes, then u^ = denotes the tangents at 
the origin ; so that if Ug has three real linear factors, the 
origin is a triple point ; and so on. 

Generally to find the multiple points of a curve, that 
is the points where the curve crosses itself, consider the 
first derived equation of (A) (§ 125) 

dx dy dx 
this gives in general a determinate value of dy/dx ; 
except when d{/dx = 0, and dt/dy = 0, when the value of 
dy/dx becomes indeterminate, and the second, or third, 
. . . derived equation must be employed to determine dy/dx; 
so that to determine the multiple points of a curve, we 
must find the values of x and y from the equations 
d!/dx = 0, d{ldy = 0, which also satisfy the equation of 
the curve {{x, y) = 0; and then a change of origin to one 
of these points will indicate its nature by inspection. 
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Newton's Analytical Parallelogram. 

In tracing a curve in the neighbourhood of one of these 
singular points, taken as origin, it is important to know 
which terms in the equation can be neglected in com- 
parison with a pair of others. 

In Newton's Parallelogram a term of the equation 
Ax^y^ is represented graphically by a letter A at the 
point whose coordinates are m, ti on a diagram; and now if 
another term Bx^y^ is placed on the diagram, the line AB 
will divide the diagram into two parts, such that a third 
term Gx^y* can be rejected if its representative point G lies 
on the side of AB remote from the origin, as being an 
infinitesimal of a higher order than the terms A and B, 
which have been retained. 

When all the other terms of the equation can be 
rejected in comparison with Ax^y^ and Bx^y^, then 

Ax'^y''+B(cPy^ = 0, or Ax^^-P+By^-^'^O 

will give a close approximation to the curve. 

Examples. 
(1) Determine the tangents at the origin of 
(i.) xHx+y)-a\x-y)=^0; 
(ii.) aj^3/2 — a^(a;+2/)=0; 
(iii.) x^y^ — a\x^ — 2/2j = ; 
(iv.) x^+y^+eax^y-Say^^O; 
(v.) x^ + Sax^y + 2axy^ — ay^ = ; 
(vi.) a;^+2/^ = a^ 2aa;, 2ay; 
(vii.) a?+y^ = a^, a^x, cj?y, ax^, Saxy, ay^\ 
(viii.) 0?*+^* = ^*, o^x, ..., a?'x^, o?xy, aa?^ ax^y; 
(ix.) a^+y^ = a^ a*aj, a^x\ a^xy, a^a?, ax^, cia^y, a^y\ . . . ; 
and sketch the curves near 0, bv means of 
Newton's Parallelogram. 
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(2) Determine the equation of the tangent of a circle, 

parabola, ellipse, hyperbola, and generally of a 
conic at any point, by changing the origin to the 
point and back again. 

(3) Prove that chords of a conic, which subtend a right 

angle at a fixed point in the conic, pass through 
a fixed point on the normal at 0. 

152. Asymptotes. 

To determine the nature of the curve at an infinite 
distance from the origin, consider the geometrical inter- 
pretation of the equation Un = 0, which represents n 
straight lines, real or imaginary, through the origin. 

The real straight lines will approximate to the nature 
of the curve at infinity, and will be parallel to the 
rectilinear asymptotes, if asymptotes exist. 

A rectilinear asymptote is defined to be **a straight 
line, at a finite distance from the origin, to which a branch 
of the curve continually approaches, and ultimately at an 
infinite distance becomes indefinitely near." 

The equation of an asymptote will therefore be of the 
form y = m^ + Cy when y = mx is the equation of one of 
the straight lines represented by Un = ; the asymptotes 
of a hyperbola are familiar instances. 

The problem of finding an asymptote is then, from the 
implicit relation f(cc, 2/) = 0, to expand y by reversion of 
series (§ 85) in descending powers of oj, in the form 

y = mx + c+px-'^ + qx~^+.,. 

Substituting this expansion of y in terms of x in the 
equation f(aj, y) = 0, and treating the resulting equation 
as an identity, then by equating to zero the coefficients 



V 
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of x^, a;""\ aj'^~-, ..., sufficient equations are obtained to 
determine m, c, p, ?, .... 

Then m determines the direction and c the position of 
the asymptote, while p, or in its absence g, ... determines 
the side of the asymptote on which the curve lies ; for 
this reason it is generally useful to expand y in descend- 
ing powers of x as far as three terms. 

We notice that if u„_i is absent, and y — mx is not 
a repeated factor of u^, then c = 0, and the corresponding 
asymptote passes through the origin. 

Thus, the general equation of a conic section, of ex- 
centricity e, and semi-latus-rectum Z, is 

y^=2lx+{e^-l)x\ 
when the origin is at a vertex ; and putting 

y = mx+c+px''^+,,,, 

we lind m^ = e^ — 1, c = l/m, p= '-n'^/'Z'ni, .... 

The asymptotes are therefore real only for hyperbolas, 
in which e > 1 ; and their equation will then be 

When e = l, the conic is a parabola, and the asymptote 
is given by m = 0, c = x; so that it lies at an infinite 
distance, and does not therefore satisfy the definition of 
an asymptote. 

If Un has a factor cc, then to determine the corre- 
sponding asymptote, we must expand x in descending 
powers of y, in the form 

x = c'+p'y-^+q'y-^+ ,.. ; 
or we may put, in general, 

x=^m'y+c'+p'y-'^+...y 
and determine m', c\ p\ ... as before. 
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Sometimes the expansion of y in descending powers of 
a; or aj in powers of y must be written 

and then y= Ax^+mx+c, or x = A'y^+7n'y + c' 
is called a parabolic asymptote ; and so on. 

When it is possible to obtain y explicitly in terms of 
Xy or X in terms of y (§ 13) from the implicit relation 
{(x, y) = 0, the asymptotes are then readily determined 
by expanding by the Binomial Theorem and other 
algebraical operations in descending powers. 

Thus if ic^ ^ 2/3 -_ fjfZ^ 

then 2/ = 4/(ct^ — a?) = — aj(l — a^aj"^;*= '-x+^a^x-^+...; 

or a; = 4/(a^-2/^)= -2/(1-^32/"^)*= --2/+J«.V"^+--- 
Also if x = a makes y = oo , or y = b makes aj=oo , then 

a; — a = is an asymptote; and so also is y — 6 = 0. 
For instance, if the equation of the curve is 

then y^ = bV/{x^ — a^), and x^ = a^y^/{y^ — b^) ; 
so that aj ± a = 0, and y±b = are asymptotes. 

The preceding considerations are in general suflScient 
for tracing a curve whose Cartesian equation is given, 
but considerations of symmetry are also useful ; thus if 
even powers only of x appear in the equation, the curve 
is symmetrical right and left of the axis of y : if even 
powers only of y appear, the curve is symmetrical above 
and below the axis of x, as if reflected in Ox. 

If X and y are involved symmetrically, then x — y = 
is an axis of symmetry, in which the curve is reflected, 
and on being doubled along this line the two halves of 
the curve will come into coincidence. 
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Eicamples, 
(1) Determine the asymptotes of the following curves, 
and draw them : — 



(•• • ' 
111. 

(< 

(vii/ 
(viii. 

{< 

(xi/ 



.) a? — y^=^a?\ (ii.) {x/a)^ — (y/b)^ = l; 

.) y^ = 2ax+x^; (iv.) x(x^+y^) — ay^ = ; 

.) x^y+xy^ = a^] (vi.) a^-'xy^+ay^ = 0; 

.) a^y^ = a\x^ + y^)y or a\x^ — y^)\ 
) 03^ — t/* — a^aJ2/ = 0; (ix.) xy{x^ — y^)=^a^\ 
) xy(x^+y^) = a(o(?+y^)] 

.) a/3y = c^, where a, j8, y are the perpendiculars on 
the sides of an equilateral triangle. 

(2) Determine the asymptotes of the curves in the pre- 

ceding set of examples (p. 317), and draw them. 

(3) Determine the equation of the cubic curve which has 

the asymptotes 

2x^y+Sa=0y x+y—Sa=0, x+y+a = 0, 
and cuts the axis of cc at at an angle tan"^( — 2). 

(4) Draw, for different values of c, the curves 
(i.) x^+2cxy + y^ = l; 

(ii.) x^ + 2cxy+x^=l^c^, 

equivalent to 
c = xy + J{l-'X^)J{l-y% or xy+^{x^-^l)^(y^^l); 

(iii.) x^ + 2cxY + y^=^; 
(iv.) x^+2cxY+y^ = l'-c^, 

(5) Draw the curve (x^-a^f+iy^-b^f^c^ 

(i.) c'>a^ + b\ (ii.) a^+¥>c^>a\ (iii.) c^ = a\ 
(i V.) a^>c^> b\ (v.) c* = b\ (vi.) 6* > c\ 

(6) Draw the curves 

(i.) tan - \c/a + tan " ^y/b = tan - ^c ; 
(ii.) tanh " ^x/a + tanh " ^y/b = tanh " ^c ; 
(iii.) (tanflj)2+(tan2/)2 = J, or 3. 
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153. Polar Coordinates. 

Id discussing properties of a curve coDnected with 
straight lines radiating from any origin 0, it is convenient 
to change to polar coordinates (r, 6) by putting (§ 22) 

aj=rcos0, 2/ = ^8inft 

Substituting in the rational integral equation (A), p. 315, 

u/-'*+(^ cos 0+5 sin 0)r-'*+i 

+ (acos20+2A sine cos 0+6 sin2%-'»+2+...=0; 

an equation of the n^ degree in r~^; so that a straight line 
cuts a curve of the n^ degree in n points, some pairs of 
which however may be imaginary. 

Denoting by r^, r^, rg, ... r„ the roots in r of this equa- 
tion and their harmonic mean by r, then 
r-^=^(r^-^+r^-^+..,+rn'^)/n=-'(Acose+Bsine)/nvQ. 
or ilrcosfl+firsin 0+nuQ = 0, 

Ax+By + nuQ = 0, or u^ + nuQ==Oj 
the equation of a straight line, the locus of a point P, such 
that OP is the harmonic mean of OP^, OP^y ... OPn» where 
Pj, Pg, ... Pn are the n points in which the straight line 
cuts the curve ; this straight line is called the polar line 
of 0> by analogy with the polar line of a conic section, 
with which it coincides when n = 2. 

The straight line Ui+Uq = is the locus of P when 

r-i = ri-i+r2-^+...+rn-^. 
In the same way the polar conic of 0, 

Ug + (n - l)ui + -y "2 X = ^» 
is the locus of P when 



and so on. 
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Writing, as usual, u for 1/r, then since for the points of 
intersection with the curve, 

2u = — (4 cos fl + -B sin 0)/nuQ, 

therefore licPu/de^ + u) = 0. 

But cPu/d6^ +U — C cosec^0, 

by § 93, where c denotes the curvature, and <f> the radial 
angle at which the vector OP cuts the curve ; so that 

2c cosec^0 = 0, 

at the points of intersection of OP with the curve. 
(Dr. Routh, Quarterly J. of Math,, xxiv., p. 257.) 

We have found that in a central field of force (§ 84), 
the acceleration to the centre 

P = h^u\cPulde^+u\ 

so that the orbit is a straight line, and d^u/d6^+u = 0, 
when P = 0; but the orbit is concave to the origin when 
P and therefore d^u/dO^+u is positive ; and convex when 
i they are negative; and at a point of inflexion, where the 
curve changes from concavity to convexity, or vice versa, 
d^u/dO^+u vanishes and changes sign. 

i Definition, — A curve is said to be concave with respect 
I to a point or line when it lies on the same side of its 
tangent as the point or line in the neighbourhood of the 
point of contact ; and convex when it lies on the opposite 
side of its tangent ; and at a point of inflexion the curve 
crosses its tangent. 

In interpreting the above, some of the p's must be 

^ negative, and we shall take p as positive or negative 

according as the curve is concave or convex to the origin. 

Incidentally we deduce that the three points of inflexion 

on a cubic lie in a straight line. 
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154. Equation of the Chords Tangent, Asymptote, and 
Normal of a Curve in Polar Coordinates, 

It is convenient to employ u, the reciprocal of r (§ 23), 

and now the equation of a straight line can be written 

u = A cos 6+ B sin 6; 
or, more generally, 

u=: A cos {6 — a) + B sin (0 — a\ 

equivalent in Cartesian coordinates to 

l=:A(x cos a + 2/ ^^^ a)+-B(t/cos a— ajsin a), 
an equation of the first degree is x and y, and therefore 
the equation of a straight line. 

To find the equation of the chord of the curve, whose 
equation is u = fd, which passes through the two points 
whose vectorial angles are a±^, we must determine A 
and B from the equations 
f(a+^j = ^cos/3+5sin^, f(a-i8) = ^cos^-£8in)8; 

and therefore 

,_ f(a+/8)+f(a-|8) p f(a + /3)-f(a-^) . 

2coS|8 ' 2sin;8 

so that the equation of the chord is 

2 cos p ^ ^ 2 sin p 

To determine the tangent where = a, put /3=0 ; then 

^ ^f(a+/3)+f(a-^) ^ , ^ /(a+)8)-f(a-i8) ^ j,^ . 
2cos/3 ' 2sin/3 ' 

and therefore the equation of the tangent at = a is 

u = fa. cos(0 — a)+fa.sin(0 — a) (ii.); 

so that, if 0P= 1/fa, then 0^= - 1/f a (fig. 10) ; 

and looking along OP from 0, then 0^ must be drawn to 

the right if fa is negative ; to the left if positive. 

Suppose fa = 0, but fa is finite ; the point of contact is 
then at an infinite distance from 0, but the tangent 
remains at a finite distance, and the tangent is therefore 
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an asymptote, its equation being 

u = fa.8in (0 — a) (iii.)- 

If fa = 00, the asymptote passes through 0, and its 
equation is given by = a. 

The equation of the normal will be 

ti' = fa.cos(fl — a) — ^^ sin(0 — a) (iv.) ; 

1 a 

for it is the equation of a straight line through the point 
of contact of the tangent, where Q = a\ and this straight 
line is perpendicular to the tangent, as is readily seen by 
changing to Cartesians. 

When the equation of the curve is given in the form 
r = F0, the equation of the normal will be found to be 

l_ cos(0~a) sin(0~a) 

7"~ Fa "^ F'a ^ ^^ ' 

so that Og = Fa, if OP = Fa (fig. 10). 

For instance, the equation of the chord, tangent, and 

normal, of the conic lu = l+e cos 6, are 

lu = e cos 6 + sec /3 cos(0 — a), lu = e cos 6 + cos{6 •— a). 

Zu = l±^^5^{esin0+sin(0-a)}. 
sma 

Examples. 
(1) Find the asymptotes of 

(i.) r = asec6, 6cosec0, asecO+bcosecO, asec(0— a), 

b cosec (6 — )8), a sec (0 — a) + 6 cosec (6 — /3). 
(ii.) r = a tan 0, a cot 6. 

(iii.) r = a(sec 0— cos 0), a(8ec fl+tan 0) ; and deter- 
mine their equations in x and y, 
(iv.) r = sec 20, cosec 20, sec 30,... secnO, cosec 7i0. 
(v.) r"2= a^sec 20, a^cosec 20. 

(vi.) r = a + b cosec (the conchoid of Nicomedes). 
(vii.) r = a/e, ae/{6-a), (10^/(6^- a% asmmO/smO. 
(viii.) l/r= ±l+sec/3cos (the hyperbola). 
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Fig. 46. 



155. Pedal Curves, 

The locus of Y, the foot of the perpendicular on the 
tangent of a curve drawn from the origin 0, is called 
the pedal of the curve with respect to 0, and is called 
the pole of the pedal. 

Thus with respect to a focus, the pedal of an ellipse or 
hyperbola is the auxiliary circle, and the pedal of a para- 
bola is the tangent at the vertex of the parabola. 

If OYP is a rigid right angle, of which OY passes 
through 0, and YP touches the curve, then 7, the foot 
of the perpendicular from on the normal at P, is the 
centre of instantaneous rotation of the right angle (§ 21) ; 
so that /F is the normal of the locus of Y. 

Since lY is a diameter of the circle described on OP 
as diameter, it follows that the envelope of circles 
described on the variable vector OP as diameter is the 
pedal with respect to 0, the locus of F. 
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This can easily be generalized for the case of a rigid 
angle PYP' touching two fixed curves, at P and P' ; the 
centre of instantaneous rotation will be at I, the point of 
intersection of the normals at P and P', and therefore 
/Fis the normal of the locus of Y; and since the circle 
circumscribing the triangle PYP* has the same normal 
at r, it follows that the locus of Y is the envelope of 
these circles, described round the triangle P YP\ 

Returning to the pedal of the curve (fig. 46), then 
since the angle YI is equal to the angle OPI in the same 
segment, it follows that the pedal curve cuts F at the 
same angle as the curve cuts OPy or the pedal and the 
curve have the same radial angle <j> (§ 22) at correspond- 
ing points. 

Therefore, denoting F by p, and the angle -40 F by w, 

cot <A = 7YTy=— {— , so that YP= f-. 
^ OY pcUo dw 

We denote 01 by g, so that the polar coordinates of 
1 are (g, ylr) and q^dp/dcoj ^ = a)+^7r; and since Q, the 
centre of curvature at P, is the point of contact of the 
normal PI of the curve AP, or of the tangent PI of the 
evolute BQ at Q, therefore the locus of I is the pedal of 
the evolute with respect to 0, and IQ = dq/d\fr = d^p/d(a\ 
Therefore the radius of curvature at P, 

ds/dw = p = 2^ + dq/d\[r =p + d^p/dw^ =p+ qdq/dp. 
The equation of the tangent YP being 

xco8co+y sin (iD=p; 
and of the normal PI being 

X cos\fr + yam\fr = q or — ccsin w + y cos co = c^p/dco ; 
it is therefore derived from the tangent by difierentiating 
with respect to co. 
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The equation of the normal to the evolute at Q will 

th erefore be —x cos w •— y sin w = cPp/dco^ ; 

so that the coordinates of Q are 

. dp cPp dp . d^p 

— sin ft)^ — cos ojt-^, cos ft) j- — sin ft^^-^- 

Ctft) . Otft)'^ (Xft) cUo^ 

The relation p = fft), connecting p and o) for the curve 
AP, is called its tangential polar equation ; and j9 = fft) is 
the polar equation of the pedal curve A F, with p and 
ft) as polar coordinates; and then g' = r(a) + Jtt) will be the 
polar equation of the pedal of the evolute BQ. 

The original curve AP is the envelope (§ 105) of 

X cos (o + y sin ft) = fft) ; 
and with reference to the curve A F, the curve AP is 
called the first negative pedal of AY with respect to 0. 
Thus the first negative pedal of a circle is an ellipse or 
hyperbola, according as the pole is inside or outside the 
circle; the first negative pedal of a straightlineisaparabola. 

For instance the directions of motion of the parts of a 
rope in contact with a moving pulley are at any insiant 
tangents to a conic, of which a focus is at the instantaneous 
centre of rotation of the pulley. 

We may take the pedal of the pedal, positive or 
negative, with respect to the pole~ 0, any number of 
times n, and then we obtain a curve called the n*^ pedal, 
positive or negative. 

Suppose rays of light proceeding from are incident 
on a reflecting curve AP ; the reflected ray will pass 
through ^ on OY, produced so that 0Z=20Y, and will 
be parallel to YI, and therefore normal to the locus of Z; 
the envelope of the reflected rays, called the katakaustic 
(§ 105), is therefore the evolute of the locus of Z, a curve 
double the size of the pedal curve A Y, 
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(i) 



(•••') 




Fig. 47 

156. LimagonSj the Pedals of a Circle. 
Consider the pedal of a circle, centre C and radius a, with 
respect to any point 0, where 00 =b. 

Then p = OY=IT+OI=a+b cosca, 

the equation connecting the polar coordinates p and co of 
the locus of F, and this curve is called a limagon. 

The limaQon is called a conchoid of a circle, because it is 
described by producing the vector of the circle r = b cos 6 
a constant distance a. 

The equation of the lima^on may be written 

p= ±a + b cos ft), 
corresponding to parallel tangents of the circle ; so that the 
chord of the limaQon through is of constant length 2a. 

li b < a, is inside the circle, and the pedal consists 
of a single oval curve (fig. 47, i.). This oval has points of 
inflexion, if 6 > Ja. 

If 6 > a, is outside the circle, and the pedal is looped, 
having a double point (fig. 47, ii.). 

If 6 = a, is on the circumference, and p = a( 1 + cos ft)), 
the equation of a cardioid (fig. 47, iii.). 

If in fig. 8 we fix the bar PQR, and move the cross 
OQRy then any point fixed in OQ or OR will describe 
a lima9on. 
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157. Orthoptic and laoptic Curves. 

If two tangents PE, QR to the curve APQ intersect at 
a constant angle a, in the point R, the locus of R is called 
am isoptic curve of the curve APQ; and if the angle a is a 
right angle, the locus of R is called the orthoptic curve. 

Thus the isoptic curve of a circle is a concentric circle 
the orthoptic locus of an ellipse or hyperbola is a circle 
sometimes called the directoi^ circle ; the orthoptic locus 
of a parabola is its directrix, while the isoptic locus of a 
parabola is a confocal hyperbola. 

The equation of the isoptic locus for tangents inclined 
at an angle a is obtained by eliminating ta between 
X cos(ft) +\a) + y sin(ft) + ^a) = f(ft) + \a\ 
X cos(a) — \a) + y sin(a) — \a) = f(a) — \a) ; 
and a = |7r gives the orthoptic locus. 

The normal RI of the isoptic locus at R will pass 
through 1 the point of intersection of the normals at the 
points P, Q; since 1 is the centre of instantaneous rotation 
of the constant angle PRQ. 

Examples, 

(1) Prove that the equation of the pedal of an ellipse 

with respect to the centre is 
p^ = a^cos>^(it) + bhin^w, or {x^+y^)^ = d^x^+b'^y^. 
Show that the pedal has points of inflexion when 
b^/a^ < ^, or 6^ > I ; and that, when 6 = 0, the pedal 
reduces to two circles. 

(2) Prove that the pedal of a parabola with respect to 

the vertex is the cissoid 

p = a (sec CO — cos to), or y^ = x^/(a — x). 

(3) Prove that the isoptic locus of a parabola is a hyper- 

bola, of excentricity sec a ; and explain how the 
two branches of the hyperbola are formed. 
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(4) The isoptic locus of an ellipse is given by 

(aj2 + 2^2 _ ^2 _ 52)2 = 4 coi^a(pt?y^h^ +x^- a^b^). 

(5) The isoptic locus of a cycloid is a trochoid. 

(6) The orthoptic locus of x^+y^ = a^ is ^^2 a cos 2ft 

(7) The orthoptic locus of the cardioid is composed of a 

circle and a lima^on (Wolstenholme). 

(8) Prove that if p^, p^, p^, ... are the radii of curvature 

of the envelopes of the sides of a polygon, whose 
sides are a, 6, c, ... then ap^ + hp^+cp^+ ... is 
equal to twice the area of the polygon. 

(9) Prove that the polar equation of the curve OP of 

§ 119, in the neighbourhood of 0, is 
r=^2p6+ippe^+.,.. 

(10) Prove that the pedal of an involute of a circle, with 

respect to the centre of the circle, is a spiral of 
Archimedes (fig. 36). 
Apply this to the theory of a weighing machine, show- 
ing the weight on a dial provided with equal graduations, 
when the body is weighed against a fixed weight, sus- 
pended by a rope which wraps on the involute of a circle. 

158. Roulettes, 

When a curve, carrying a point P fixed to it, rolls on 
a straight line (or any given curve), the path traced out 
by the point P is called the roulette of P with respect 
to the straight line (or given curve). 

Thus, when a circle rolls on a straight line, the roulette 
of a point on the circumference is a cycloid^ and the 
roulette of any other point fixed in the plane of the 
circle is a trochoid (§ 21). 

An involute of a curve (§ 95) is thus the roulette of 
a point on a straight line which rolls on the curve. 
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Fig. 48 

A remarkable analogy, pointed out by Steiner, exists 
between the roulette of a point with respect to a straight 
line and the pedal of the rolling curve with respect to 
the point as pole. Steiner's Theorems assert that 

(i.) The length of the arc of the roulette is equal to 
the length of the corresponding arc of the pedal ; 

(ii.) The area bounded by an arc of the roulette, the 
ordinates at the ends of the arc, and the straight line on 
which the curve rolls is twice the area bounded by the 
corresponding arc of the pedal and the vectors from the 
origin to the ends of the arc. 

For, if AP is the roulette of the point P when the 
curve is rolled on the straight line Ox (fig. 48), and if 
PM is the perpendicular from P on Ox, the tangent at 1 
to the rolling curve, then relatively to P the locus of M 
is the pedal of the rolling curve with respect to P ; and 
therefore relatively to M the locus of P is the same 
curve, and the subnormal IM of the roulette is the q or 
dp/dco of the rolling curve, or the polar subnormal of its 
pedal. 
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We may suppose the pedal A'P rolled on the roulette 
AP, so that M is always vertically over P if Ox is hori- 
zontal ; and the pedal, if loaded so that the centre of 
gravity is at M^ will rest in neutral equilibrium on the 
roulette, provided the friction is sufficient, or else that 
teeth are cut, to prevent slipping. 

The arc -4P of the roulette will then be equal to the 
corresponding arc A'P of the pedal, which is Steiner s 
first theorem. 

The arc a in the roulette is therefore the same function 
of y as in the pedal of p. 

Also if the pedal is rolled into a consecutive position 
so that M comes to M\ and the point p of the pedal 
comes into contact with the point P' of the roulette, then 
the element MM'P'P, which is the increment of area of 
the roulette, is ultimately double the element MPp, 
which is the increment of area of the pedal, or 

lt(area Jl/JlfTT)/area MPp = 2 ; 

and therefore, by integration, the area OMPA of the 
roulette is double the area A' MP of the pedal, which is 
Steiner s second theorem. 

Exam'ples, 
(1) Draw the figures and compare the arcs and areas of 
the following pairs of curves, the first curve of a 
pair being fixed, and the second rolling on it, so 
that its pole describes a straight line Ox (or Oy) ; 
thus illustrating Steiner's Theorems. 

(i.) The parabola 2/^ = 2te, and the spiral r = W. 
(ii.) The circle aj^ + ^2 __ ^,2^ ^nd the circle r = a cos ft 

(This principle is employed in the parallel motion of 
DeleuiFs air pump. Deschanel, Physics.) 
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(iii.) The elUpse (x/ay + (ylby=l, 

and r=b cos{b6/a) or a cos(a6/b), 
,) The hyperbola (aj/a)2-(y/6)2=l, 

and r = a cosh (aO/b) or b sinh (bO/a). 
.) The exponential curve y = 6e*/S 

and the reciprocal spiral r=c/6. 
.) The cycloid and the cardioid. 
.) The trochoid and the lima9on. 
.) The straight line y = x tan a, 

and the equiangular spiral r=a6^**^*. 
.) The catenary y = a cosh x/a, 

and the straight line r=a sec ft 
.) The modified catenary y = 6 cosh x/a, 
and the Cotes's spiral r cos (bO/a) = b. 
.) The sinusoid y = b cos aj/a, 

and the Cotes's spiral r cosh (6fl/a) = b. 
.) The sinusoid y = a'-ae cos (03/6), 
and the ellipse ?/r = l+«cosft 
.) The hyperbola xy = c^, and the spiral r^0 = Jc^. 
.) The semicubical parabola ay^=Qi?, 

and the spiral r0^=8a. 
.) The curves {xjaf^^iylbY and r'^-n^rrrc'^-w 
.) The curves aj*~y^ = a'"+~ and r'^^rnQm.^fm-\-n 
Rectify these curves, when possible ; and prove that 8 
is the same function of y in the first curve as of r in the 
second ; also that ydy/dx and ydx/dy, the subnormal and 
subtangent, are the same functions of y, as dr/dO or 
r^dO/dTf the polar subnormal and subtangent are func- 
tions of r. 
(2) Prove that the roulette of the pole of 

l/r=l + sec a cos (6 sin a), or 1 + sech a cosh(0 sinh a), 
with respect to a straight line, is a circle. 



(iv.; 
(v.; 
(vi.: 

(vii.] 

(•♦• > 
viu.^ 

(ix/ 



{< 



(xi/ 

(xii/ 

/ • • • 
(xin. 

(xiv. 



(xv/ 

(xvi.] 
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(3) Prove that the curvature of the roulette with respect 

to a straight line is d sin <f>/dp, where tp is the 
radial angle and p the perpendicular on the tangent 
from the pole on the rolling curve. 

(4) Prove that if the centre of an ellipse is fixed at a 

distance b from a plane, and the ellipse is rolled 
on the plane, the point of contact will describe 
the Cotes's spiral r cosh (aed/b) = ae. 

159. Centrodes, 

When a moving plane figure slides or turns on another 
plane, which may be considered fixed, then a point / in 
the moving plane can always be found which has no 
velocity; this point / is called the centre of instantaneous 
rotation (§ 21), and the relative motion of the two planes 
is assigned by the position of / and by the angular 
velocity n of the moving plane round /. 

The point / will in the general case describe a curve in 
the moving plane and a curve in the fixed plane, and the 
motion of the moving plane will be given by rolling the 
first curve on the second ; so that any point carried by 
the rolling plane will describe a roulette of the first curve 
with respect to the second curve ; these curves described 
by / are now called the centrodes of the relative motion 
of the two planes. 

Another point J can always be found of which the 
acceleration is zero ; and now the acceleration of any 
other carried point P at a distance r from J will be 
composed of component accelerations tiV towards J and 
nr perpendicular to JP, n denoting the angular acceler- 
ation dn/dt 
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Fig.49 



Therefore the resultant acceleration of any point, such 
as P, makes the same angle a = coi~\nj'n'') with JP\ and 
the lines of resultant acceleration at any instant are 
equiangular spirals of radial angle a, the magnitude of 
the acceleration at distance r from J being ii^-coseca. 

Denoting by R the radius of curvature PP" of the 
trajectory o£ P, by v the velocity of P, and by the 
angle IPJ, then (g 92) the normal acceleration of P, 
\^jR=n^r C03 6 — nr sin Q = n^r cosec a sin(a — 0). 
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The points whose trajectories have zero curvature, and 
which are therefore at this instant describing straight 
paths (as M in fig. 48), or rather are passing through 
points of inflexion on their trajectories, are obtained by- 
putting = a, and therefore lie on a circle /QJ' passing 
through I and J, called the circle of inflexions (the circle 
PMI in fig. 48) ; and this circle will touch the centrodes 
at /, since the acceleration at M is in the normal. 

Now if IP meets this circle in Q, then the curvature 
of the trajectory of P is 

1/jB = tiV sin(a — 0)/v%in a 

= n\PQIv' = PQ/PI^. 

Thus PF^PPjPQ, 

IP'=PP/PQ-PI=^PI.1Q/PQ, 

IP' pi.iQ iQ pr 

or L + Jl =Jl. = !5^ 

PriF IQ ID ' 

where ID is the diameter of the circle of inflexions, and 
where <f> now denotes the angle between IP and the 
common normal of the centrodes at /. 

A point P inside the circle of inflexions will thus 
describe a trajectory convex with respect to /, but concave 
if the point is outside this circle; so that the rolling 
centrode, if loaded so that its centre of gravity is at P, 
will be in stable or unstable equilibrium with IP vertical 
according as P is inside or outside this circle. 

The diameter ID of the circle of inflexions is inferred 
by placing the carried point P for a moment at C, the 
centre of curvature of the rolling centrode at /, when it 
is easily seen that the centre of curvature of the roulette 
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through O-will be at (7, the corresponding centre of 
curvature at 1 of the fixed centrode ; and now, with ^ = 0, 

ID Cria p^p'' 

or ID — ppj{p + p), where />, p denote the radii of curva- 
ture Gly IC of the rolling and fixed centrodes, reckoned 
positive when the centrodes are convex to each other. 



Now ^+^ = sec<^(i + ^^,), 



P P 

and the symmetry of this relation shows that P is the 
centre of curvature of the roulette of P' with respect to 
the former moviug centrode ; also that OP, CP' intersect 
in a point T on IT the perpendicular to IP, 

When a curve rolls symmetrically on an equal curve, 
the roulette of any point will be similar to the pedal of 
the curve with respect to P, but enlarged to twice the 
scale ; the reason being that the reflexion of the carried 
point P in the tangent at / is a fixed point. 

Suppose for instance, as drawn in fig. 49, that the 
fixed and rolling centrodes are equal ellipses, and that the 
carried point P is at one of the foci; the roulette of P 
will be a circle, and the centre of curvature P' will lie at 
a focus of the fixed ellipse. 

We can pivot these ellipses at the other foci and 
0', and now revolve them in contact with each other ; 
teeth may be cut to prevent slipping, and the revolving 
foci P and P' may be connected by a link to prevent 
separation. This mechanism is sometimes employed, and 
the ellipses roll on each other as if connected with a 
crossed parallelogram of bars OP, PP\ PV\ OV. 
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160. Tlve Area of a Roulette. 

nAs / moves along the rolling and fixed eentrodes with 
the same velocity dsjdt^ the moving plane wiD turn in 
the time dt through an angle <U^ which is the sum of the 
curvatures dsjp and dsp of the equal ares ds of the 
rolling and fixed eentrodes ; so that 

dt \p p/dt 

The normal PI of the roulette of P now sweeps out an 
area, denoted by (P/), at a rate which, by Guldin s Theorem 
generalized (§ 62), will be measured by the product of PI 
and of the component velociiy of the middle point of PI 
perpendicular to PI. 

This component velocity, being the arithmetic mean of 
the velocities of / and P in the same direction, is equal to 

J cos (f>(dsldt)'\-^n . PI ; 
so that 

Of these two terms, the first is the rate at which PI 
sweeps out polar area on the rolling centrode with respect 
to the pole P (§ 56) ; and the second represents the 
growth of the M.i. round P of the perimeter of the rolling 
centrode, supposing the perimeter replaced by a wire of 
variable density ^(c+c'), the arithmetic mean of the 
curvatures c = l/p and c' = l/p' of the rolling and fixed 
eentrodes. 

Thus if the rolling centrode returns to its original 
position after one or more complete revolutions, the area 
of the roulette of P exceeds the sum of the areas of the 
rolling and fixed eentrodes by the moment of inertia of 
this perimeter of variable density of the rolling centrode. 
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By taking the carried point P at G, the centre of 
gravity of this perimeter, we obtain the minimum 
moment of inertia (§ 64), and therefore the roulette of 
minimum area; and the area (P) of the roulette of P 
will exceed the area ((?) of the roulette of Ohy M . GP^, 
where M is the mass distributed on the perimeter. 

But day/dt denoting the angular velocity of the rolling 

centrode, M =f\{G + c')cZs =f\dw = titt, 

for n complete revolutions of the rolling centrode; so that 

(P)-((?) = n7r.(?P2. 

161. Theorems of Holditch, Elliott, Leudesdoif, and 
Kempe. 

A simple relation connects the areas (P), (A), (R) swept 
out by three points P, J., ii in a straight line on a bar, 
which we may take to be the bar of the planimeter of 
fig. 45, carried round by the rolling centrode, and making 
n complete revolutions. 

For if OK is the perpendicular from G on ^P produced, 

then (P)^(A)==n7r(0P^-GA^) = mr{KP^-KA^) 

= nirh{KP+KA\ . 
and {A)- (R) = n7rc(KA + KR\ 
since AP = 6, RA = c ; so that 

or (6 + c)( JL) - h{R) - c{P) = nirhc(b + c) ; 

this is called HolditcKa Theorem. 

If the closed contours described by A and P lie entirely 
outside each other, the bar AP can only oscillate between 
two extreme positions, and n=0. 
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In using the planimeter, the pivot is generally fixed 
outside the contour described by P, and the joint A 
oscillates on the arc of a circle of radius a, so that 

(A) = and ti = 0; and thus h{R)+c{P) = 0, 

attending to the sign of the area. 

Thus if P describes a circle of radius r not enclosing 0, 
the motion of the bars OA, AP in similar to that of the 
beam and connecting rod of a steam engine; and now 
(jR)= — 7rr^c/6; this is independent of the length of the 
beam, and is therefore the same for the ordinary direct 
action steam engine, as applied in the locomotive engine, 
where A oscillates in a straight line. 

A new form of Planimeter has been lately brought out, 
in which the end A of the bar AP is constrained to move 
in a straight or curved slot, while the pencil P is carried 
round the area to be measured; but now the area is 
registered by the motion of a wheel W which can slide 
and turn on a graduated round bar CD^ projecting at 
right angles to AP from a point 0, which we may take 
to be the middle point of ^P. 

The rate at which AP sweeps out area is the product 
of AP into the component velocity of C perpendicular to 
AP ; and this component velocity is equal to the velocity 
with which the wheel W slides along GD\ the wheel 
being supposed to roll, but not to slide on the paper ; in 
this manner the sliding motion of the roller R on the 
paper of Amsler's Planimeter is obviated. 

Mr. Elliott has generalized Holditch's theorem by 
showing that the bar PAR may be replaced by an elastic 
thread, which stretches, but keeps the ratio hjc constant, 
when a similar Theorem holds. 
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Mr. Leudesdorf has extended the theorem to the relation 
connecting the areas of the roulettes (J.), (5j, (C) of any 
three points A, By C with the area (P) of the roulette of 
any fourth point P ; and he finds that 

{P) — x(A) — y(B) — z(G) = mr times the rectangle 
on the segments of a chord through P of the circle circum- 
scribing the triangle ABC, when Xy y, z denote the ratios 
of the triangles PBC, PGA, PAB to the triangle ABC. 

For {P)-{A) = nir{GF^'-GA^), ... ; 

so that (P) -x{A)'-y{B)-z{C) 

= nir(GP'^-x . GA^^-y . GB^-^z.GC^) 
=^n7r(x . BC^+y . CA^+z , AE^) 
which can be shown to lead to Mr. Leudesdorf s result; and 
reduces to Holditch's Theorem when P lies in a side BG. 
Mr. Kempe points out that the locus of P, for which 
(P) is zero, is a circle ; and the locus of P, for which (P) 
is constant, is a concentric circle, exceeding in area the 
first circle by an amount proportional to (P) ; but that if 
71 = 0, this system of concentric circles must be replaced 
by a system of parallel straight lines. 

(Messenger of Mathematics, vol. vii., 1878.) 

162. The Envelope of a Garried Line, 

Similar theorems hold for the envelope of a carried line 

NXy fixed in the moving centrode and carried round by it. 

The point of contact N on the envelope is the foot of 

the perpendicular IN drawn from / on the carried line ; 

and now if we denote by y the ordinate IN of the point 

/ of the rolling centrode with respect to the line Nx, by 

X the corresponding abscissa on the line Nx, and by S the 

length of the arc described by iV ; then 

dS __dx . dco 
dt'~di'^ydi' 
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so that S^x =:Jydw ^f\(c + (!)yd8 ; 

or 8 the length of the arc of the envelope exceeds x the 
projection of the corresponding arc of the rolling centrode 
on the carried line by the moment about the carried line 
of the same distribution of density \{c + d) as before. 

Carried lines which have envelopes of the same length 
will therefore, touch a circle, with the centre of gravity 
O of the perimeter as centre. 

Again, if 0' denotes the angle between /iV and the 
common normal of the centrodes at /, the radius of curva- 
ture R of the envelope of the carried line at 'R is given by 

n, dS dx , .fds , 

supposing Nl meets in Q' the circle 1Q'D\ which is the 
reflexion of the circle of inflexions in the common tangent 
of the centrodes at /. 

Therefore, for all carried lines passing through D\ 
R^ = 0, and the envelope has a cusp lying on this circle, 
called for this reason the circle of cusps ; while in addition 
the circle is the locus of the centres of curvature of all 
carried lines, the centre of curvature Q' being obtained by 
dropping the perpendicular UQ on IN. 

The normal NI of the envelope will now sweep out 
area {NI) at a rate 

dt ~^'Adr dtJ'^ydt^^y'di' 

so that the area swept out by the normal NI exceeds the 
area swept out on the rolling centrode by the ordinate 
NI by the M.i. of the perimeter, of variable density 
^(c+c'), round the carried line. 
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Thus from Theorem II., p. 128, it follows that if two 
carried lines intersect at right angles in P, the sum of the 
areas of their envelopes exceeds the area of the roulette 
of P by the area of the rolling centrode. 

Carried lines which sweep out equal areas are such 
that the M.i. of the perimeter round them is the same ; 
they are therefore equimomental lines, and therefore 
tangents of an equimomental ellipse; and by varying 
the area we obtain a system of equimomental confocal 
ellipses, by well known theorems. 
(Besant, Roulettes and Olissettes ; 

Kempe, Messenger of Mathematics ; Nature, 1878 ; 

Darboux, Bulletin des Sciences Mathematiques, 1878; 

Larmor, Nature, 1881 ; Proc. Gam. Phil, Soc, 1890; 

Minchin, Uniplanar Kinematics, 1882.) 

163. When the fixed centrode becomes a straight line, 
c' = ; and the area swept out by the normal PI 
(fig. 48) in a complete revolution of the rolling curve 
is the same as that swept out by the ordinates MP, and 
is therefore double the area of the pedal of the rolling 
curve with respect to P. 

Now, if the rolling curve is coated with matter of 
linear density |c per unit of length, the mass will be 



l\cds = /\d(a = TT ; 



and therefore the area of the roulette of P will exceed 
the area of the roulette of G, the centre of gravity of this 
wire, by tt . GP^ ; and consequently the area of the pedal 
of the rolling curve with respect to P will exceed the area 
of the pedal with respect to (? by ^tt . GP^, 



CURVES IN GENERAL, 345 

Thus, for instance, the pedal of a circle with respect to 
its centre being the circle itself, will have an area iro? ; 
and consequently the area of the lima9on p = a+6cosco 
will be Tra^+^irh^, while the area of the corresponding 
trochoid will be 2ira?+irh^\ and with a = 6, the area of 
the cardioid p = a(l+cosft)) will be f Tra^, and of the 
cycloid ^Tra\ 

In a roulette with respect to a straight line, 

1 1 _sec0 

or NR^J^^ = pEcos<l>, 

or 1_1 pcos^ 

WN iV'2 * 

If the rolling curve is an ellipse and the carried point 
at a focus, then (fig. 48) 

p co8<f> = N(2a-N)/a = 2N"Nya, 

^1, 4 111 

so that -Fi+^r=-- 

it N a 

If this roulette sweeps out a surface by revolution 
round the fixed straight line, the curvature of this surface 
will be everywhere 1/a, a constant; a soap-bubble film 
of revolution will assume the shape of this surface, or a 
surface similarly generated by the roulette of the focus 
of a hyperbola or parabola. 

For instance, the roulette of the focus of a parabola is 
a catenary and we obtain the catenoid of fig. 16. 

Again, the roulette of the pole of the involute of a 
circle with respect to a straight line is a parabola ; and 
the roulette of the centre of the rectangular hyperbola 
7^cos26 = a^ is a curve in which R — ^N, or Ry = ^a\ an 
elastica or lintearia. 
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164. Ejncycloids mtd Hypocycloids. 

These curves are the roulettes of a point P on the cir- 
cumference of a circle which rolls on the outside or inside 
of a tixed circle. 

Let denote the centre and a the radius of the fixed 
circle, C the centre and c the radius of the rolling circle, 
and / the point of contact of the circles ; then IP is the 
normal of the roulette of P, because / is the centre of 
instantaneous rotation of the rolling circle (fig. 50). 




Fig-BO. 
Draw the diameter PCQ of the rolling circle, and sup- 
pose Q originally in contact with the fixed circle at B, 
and that P is then at ^ ; ^ is then an apse and B a cusp 
of an epicycloid (§ 104). 
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Denoting the angle xOI by Q (radians), then the arc 
]B = a6=^eircIQ, so that the angle ICQ = a6/c; and the 
coordinates of P in terms of 0, for the epicycloid, are 

x = (a+ c)cos 6+c cos(l + a/c)6, 
2/ =:(a+c)sin 6+C8m(l+a/c)6 ; 

and for a hypocycloid, change c into — c. 
Thence we find, for the epicycloid, 

^=2(a+c)cos2^, 

and integrating, the arc AP 

=8=4-(a+c) sm ^ ; 

while the arc 5Q = 4- (a + c) vers ^ . 

With the notation of 'p and co of § 154, 

p = (a + 2c)cos(a0/2c), and « = (1 + a/2c)0 ; 
so that p = (a+2c)cos{aft)/(a+2c)}, 

the polar equation of the pedal of an epicycloid, of the 
form r = 6 cos 7710. 

Ai«^ T%\r ^V .Oft) . aO , 

Also PY= — :/- = asin — — r-=asin-^; 

dco a+'2c zc 

so that r2 = 0r2+PF2 

= (a + 2c)^cos2^ + a^in^^ 

r^-a-4c(a+c)cos^g=M^p., 

of the form 7^ — a^ = ( 1 — 7n?)p\ m = a/(a + 2c), 
the relation connecting p and r in the epicycloid; proving 
incidentally that the roulette of the centre of an epicycloid 
with respect to a straight line is an ellipse. 

If m = 0, then c = x, and the curve becomes the in- 
volute of a circle, in which p^ = r^^a^. 
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Again s = 4-(a + c)sin — ^^, 

of the form 8 = Zsinm^, with ^ for co, i for 4c(a + c)/a, 
and m for a/(a + 2c), as in § 104. 

And p= =_A_r_j^cos— -^ = -T-=(l-wi2)p. 

165. The Teeth of Wheels, 

By cutting teeth on wheels we can make them engage 
and transmit power without slipping ; we thus secure the 
condition called perfect roughness by the theoretical 
mathematician ; perfect smoothness between two bodies 
on the other hand is sought practically by the inter- 
position of wheels. 

Suppose 0' to be the fixed centre of a wheel of radius 
0'I=a\ which is to be made to revolve in contact with 
the wheel of centre and radius a, without slipping. 

If the circle, centre C and radius c, rolls inside this 
circle, centre 0' and radius a\ and describes the hypo- 
cycloid A'P, then if the epi- and hypo-cycloids AP and 
A'P start with the vertices A and A' in contact, the two 
curves will roll and slide on each other (tig. 50) so that 
the common normal at P passes through /, and therefore 
the constant velocity ratio of the wheels is maintained ; 
for this reason the wheels of teeth are shaped by epi- and 
hypo-cycloids. 

Only a small portion of each curve in the neighbour- 
hood of a cusp is made use of to form a tooth ; and the 
tooth is completed on the circle 0' by a portion of an 
epicycloid, and on the circle by a portion of a hypo- 
cycloid, each described by the rolling of a circle of the 
same radius c\ 
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For instance, if c'=^^a, then the hypocycloid in the circle 
is given by aj = 0, y = a sin ; so that the hypocycloid 
degenerates into a straight line, a diameter; and the 
inside portion of the tooth is straight and radial. 

Any number of change wheels of a lathe may be made 
to work together indiscriminately, provided the teeth are 
shaped by rolling circles of the same radius c. 

In iig. 50 we have taken a/a' = ^y and c = c' = ^a; and 
four teeth on the circle engage with six teeth on 0\ 

When we make c and c' infinite, the teeth are shaped 
by involutes of the circle and 0' ; involute teeth have 
the advantage of preserving the velocity ratio of the 
wheels unchanged for variable distances between the 
centres of the wheels, and are employed in rolling mills. 

When the radius a' is made infinite, the corresponding 
wheel becomes converted into a rack, and the teeth on 
the rack are shaped by cycloids. 

Sometimes to ensure greater regularity of working, 
lidical teeth are employed, and now the tooth of one 
helix on the cylinder of radius a works with the tooth 
of the helix of equal pitch on the cylinder of radius a' ; 
so that when one helix rolls on another of equal pitch on 
a parallel axis, any point of the helix describes an epi- 
cycloid. (MacCord, Kinematics ; G. B. Grant, Odontics). 

166. The Double Generation ofEpi- and Hypo-cycloids, 

Produce PI both ways to meet the circles and 0' 
again in H and H\ and draw EPE' parallel to 00' to 
meet OH and O'R in E and E\ 

Then E and E' are the centres of circles, of radii a+c 
and a'—c, which touch each other at P and the circles 
and 0' at H and H' ; so that the epicycloid AP and the 
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hypocycloid A'P can be described by the rolling of these 
circles on the circles and 0'; we thus perceive that 
there is a double tnode of generation of the epicycloid and 
hypocycloid. 

Relatively to the circle 0\ any point on the circum- 
ference of the circle describes an epicycloid; and a 
hypocycloid if it is made to roll inside the circle 0\ 

The envelope of the diameter QCP is another epi- 
cycloid, described by a rolling circle of radius |c, as 
is readily perceived when we notice that the point of 
contact of the envelope is at the foot of the perpendicular 
drawn from / on PQ, 

The envelope of any other carried straight line, say 
parallel to PQ, will be a parallel curve to the epicycloid, 
the envelope of PQ^ and will therefore have an epi- 
cycloidal evolute. 

167. Epi- and Hypo-Trochoidsy or Bicircloids, 

A point fixed in CP at a distance kc from G will 
describe a curve, given by 

x=:(a+ c)cos 6+Jcc cos(l -h a/c)6, 

y = (a+c)sin 6+Jccsm{l + a/c)6 ; 

these curves are called epi- or hypo-trochoids, and some- 
times bicircloids. 

The relative orbit of two planets is a bicircloid, if the 
planets describe circles round the Sun; figures of the 
relative orbits of the Earth and the different planets, 
drawn mechanically by Mr. Perigal, are given in Proctor's 
Geometry of Cycloids. 

For if a, ^ denote the (mean) distances from the Sun, 
and 71, m the mean motions, the relative orbit is given by 
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x = a cos nt+^ cos mt, y^asinnt + P sin mt ; 

where nW = 7n?P\ by Kepler's Third Law (§ 179); so that 

a+c^ttj lcc = p, and l+a/c=m/ti=(a//8)^; 

or a = a — ^/a^ c = ^/a^. 

When c = a, or m = 2n, the bicircloids are limagons 
(§ 156), as is otherwise geometrically obvious. Thus if 
we assume that the period of Mars is two years, the 
relative orbit of the Earth and Mars will be a lima9on. 

Examples, 

(1) Prove that if an epicycloid rolls on a straight line, 

the centre will describe an ellipse. 

(2) Prove that, if a helix rolls on a straight line, any 

point on the helix will describe a cycloid. 

(3) The shadow of a helix on a plane perpendicular to 

its axis, thrown by parallel rays of light, is a 
trochoid. 

Find when the trochoid will be a cycloid. 

(4) The path of a steamer turning uniformly in a current 

will be a trochoid relatively to the land. 

(5) Show that a variable velocity ratio of two wheels 

can be attained by cutting teeth in equiangular 
spirals of the same radial angle. 

(6) If an equiangular spiral rolls on a circle, the pole 

will describe an involute of a circle. 

(7) Show that the path of the Moon relatively to the 

Sun has no points of inflexion. 
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(8) The isoptic locus of an epicycloid is an epitrochoid 

(Wolstenholme, Proc. London Math. Soc, vol. iv.). 
If the isoptic locus passes through the centre, 
its equation is of the form r = 6 sin {a0/{a + 2c)}. 

(9) Prove that the epicycloid has (i.) one cusp when 

c = a, 2a, 3a, 4a, ...; and then m = J, ^, f, ...; 
(ii.) two cusps when 2c = a, 3a, 5a, ...; and then 
m = ^, J, j^, ... ; and draw them. 

(10) Prove that the equation of the Tricuap hypocycloid 

can be written 

(i.) p = Jacos3a), (ii.) 8 = fa sin 3^, 

(iii.) (x^ + 2/2)2 ^ sax(x^ - Sy^) + 2a\x^ + y^) - Ja* = 0, 
(i V.) r* + f ar^cos 30 + 2a2r2 - Ja* = 0. 

Prove that if the tangent at P meets the Tricusp again 
in Q and iJ, 

(i.) the length Q/J = |a, 

(ii.) the tangents at Q and R intersect at right angles 
in a point T on the inscribed circle of the Tricusp, 

(iii.) the normals at P, Q, R intersect in a point N on 
the circumscribed circle, 

(iv.) NT passes through the centre. 

168. Inversion and Inverse Curves. 

When the vector OP of a curve is produced to Q, so 
that OQ is inversely proportional to OP, or OQ^^c^/OP, 
then the locus of Q is called an inverse curve of the locus 
of P with respect to the origin 0, or with respect to the 
circle of the centre and radius c ; and c is then called 
the constant of inversion. 

Thus if u = iO is the polar equation of the locus of P, 
then r = cH9 is the polar equation of the inverse curve, 
the locus of Qj with respect to the origin 0. 
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The inverse of a circle (or sphere) is another circle (or 
sphere); except when the circle (or sphere) passes through 
the centre of inversion, when the inverse is a straight 
line (or plane). 

For if OPQ meets a circle (or sphere) in P and Q, then 
OP . OQ is constant and equal to 02^ for an external 
origin 0, where OT is a tangent to the circle (or sphere) ; 
so that a circle (or sphere) is its own inverse with respect 
to any origin ; and varying the constant of inversion 
gives a similar curve (or surface) ; in this case another 
circle (or sphere). 

Maxwell calls the rectangle on the segments of a chord 
through made by the circle (or sphere) the power of 
the circle (or sphere) with respect to \ and now if the 
distances from to the centre of the circle or sphere and 
of its inverse are denoted by a and a\ and the radii by 
b and h\ then a^-h^, a'^-6'^ are the powers of the circles 

or spheres, and _==-. = -- — _ = _ — 

^ ah a^ — b^ c^ 

Any curve Pp and its inverse Qq with respect to will 
cut the vector OPQ at supplementary radial angles ; for 
since OP. OQ = c^=Op . Og, a circle can be drawn through 
PQqp, and therefore the angles pPQ, pqQ are supple- 
mentary, and Pp, Qq are ultimately the tangents at P 
and Q, when p and therefore q are brought up close to P 
and Q. 

Or, otherwise, denoting the radial angles by (f> and <f>, 

cot = d log r/dd = - d log u/d0 = - rO/id = - cot 0', 
cot0+cot0' = O, or + 0' = 7r; 
so that the radial angles are supplementary. 

Any two curves and their inverses will therefore cut 
at the same angle. 
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A small element of area is consequently unaltered in 
shape by inversion, but the alteration of linear scale is 

So also in space, an element of volume will be undis- 
torted by inversion, but will be changed in volume in 
the scale c®/^® or r'^c®. 

If accented letters refer generally to the inverse curve, 

then Ot' = c^lOg, Og'=(?IOt) 

p'=^cyPg, p = cyQg'i 

Pt.Qt'=chec?<p; 

r' . r dp , dp' dr sin d> , dr'smd> « . ^ 
p p dr ar dr dr 

and the circle of curvature obviously inverts into the 
circle of curvature of the new curve. 

Lines of curvature (§ 130) will correspond on a surface 
and its inverse ; and if a circle of principal curvature is 
drawn, the normal plane of this circle and the sphere 
through the circle and origin of inversion will invert 
into the sphere and plane through the corresponding 
circle of principal curvature. 

When the Cartesian equation of a plane curve is given 
in rectangular coordinates x and j/, the equation of the 
inverse curve with respect to the origin is obtained by 
writing c^x/(x^+y^) for x, and c^y/{x^ + y^) for y. 

Thus the inverse of the parabola y^=px with respect 
to the vertex is, writing a for c^/p, the cissoid 
ic(x'^+y^) = ay^j or y^==x^/(a—x); 
or in polar coordinates, 

r=a(sec0--cos 0), 
the locus of ^ in a circle, for an origin on the circumfer- 
ence (ex. 11, p. 42). 
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169. Mechanical Invertors and exact Parallel Motion, 
When P describes a given curve, the point Q can be 
made to describe an inverse curve by means of the 
mechanical invertors, invented by Peaucellier and Hart. 

Peaucellier's motion consists of a rhombus iP, PM, 
^Qi Q^f formed of four links of equal length, jointed at 
X, M, P, Q, and of two other equal links OL, OM, pivoted 
at a fixed point (fig. 51). 




Then, in whatever way the link motion is displaced by 
the motion of P, 

OP .OQ = OE^-EP^ = OU-LP^=:c\ 
a constant, so that P and Q describe inverse curves. 

For instance, if P is compelled to describe an arc of a 
circle by a link GP, pivoted at a fixed point C, and of 
length GP—OC^ then Q will move in a straight line 
perpendicular to OG. 

A ^parallel motion is the name given in Mechanism to 
an arrangement of bars intended to guide a point Q, 
the head of a piston rod, in a straight line ; and this is 
eflFected with exactness by Peaucellier's arrangement, 
but only approximately by Watt's parallel motion. 

Hart's parallel motion accomplishes the same purpose 
with four bars, while Peaucellier's requires six. 
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A jointed parallelogram of rods, FOHK, is taken, and 
the longer rods are crossed ; any fixed point in one of 
the rods FG is taken, and OPQR is drawn parallel to FH 
or GK, to meet FK in P, GH in Q, and HK in R (fig. 51). 

Then P, Q, R are also fixed points in the bars, such 
that when again opened into a parallelogram, OPRQ will 
form another parallelogram, the sides of which are parallel 
to the diagonals of FGHK, and the area of which bears 
a constant ratio to the area of FGHK ; so that OP . OQ 
is constant; and thus the parallelogram FGHK, when 
crossed, will act as an invertor, P and Q describing in- 
verse curves when is fixed. 

We can show Hart's parallel motion working in con- 
junction with Peaucellier's by drawing FG parallel to 
PM or LQ, GQH parallel to Oi, and FPK parallel to 
OM; and now 0, P, R, Q are the middle points of the 
bars of the crossed parallelogram FGHK. 

We may join LR and MR by bars, and now the two 
rhombuses LPMQ, OLRM are said to make a complete 
Peaucellier cell. 

When P and Q are inside the cell, the cell is called 
'positive (fig. 51), where it is shown acting as the parallel 
motion of a beam engine. 

But when P and Q are outside, it is called a negative 
cell (fig. 52) ; and now it forms a compact mechanism for 
drawing not only straight lines, but circles of very large 
radius, as required in Architecture. 

(Kempe, How to Draw a Straight Line.) 

While Q describes a straight line and P a circular arc, 
the point R will describe a curve whose polar equation is 
of the form r = a sec ± 6 cos 0, 

the inverse of a conic section with respect to a vertex. 
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Fig.52. 



If equal bars QJ5, BE are pivoted at Q, B, and a fixed 
point E in the straight line described by Q, then, as in 
fig. 8, any point D in the bar QB will describe an arc of 
an ellipse ; and the normal DI will pass through 7, fixed 
in EB produced, so that EI^IEB, 

Any complicated system of linkwork in mechanism 
can always be analysed into three bar motion, ^\x(^ as OA, 
AP, PC in fig. 45, when a bar PC is pivoted about a fixed 
point at G\ thus representing the beam, connecting rod, 
and crank in a beam engine ; the beam being practically 
of infinite length in the direct action engine. 

If OA, OP produced intersect in 7, then 7 is the 
instantaneous centre of rotation of AP ; and for uniform 
velocity of the crank P, the velocity of A will be a 
maximum when 00 and AP intersect in a point J, the 
foot of the perpendicular from 7 on AP, 

A point Q in AP will describe a. figure of 8 curve ; 
Watt's parallel motion consists essentially of a three bar 
motion, such as AO, AP, PO, by means of which the head 
of the piston rod Q is guided for a short distance on each 
side of the point of inflexion on this figure of 8. 
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The curves described by any point Q on the connecting 
rod AP in three bar motion have been studied by Roberts, 
Cayley, and Clifford ; any such curve being capable of a 
triple generation by means of bars (Clifford, Kinematic). 

When a crossed parallelogram is employed, as in figs. 
51, 52, and in fig. 49, the relative motion of the short 
bars FG and HK is imitated by rolling an ellipse with 
foci F, 0, on an equal ellipse with foci H, K\ and the 
relative motion of the long bars GH, FK is imitated by 
rolling a hyperbola on an equal hyperbola ; and relatively 
to one bar any point carried by the opposite bar will 
describe the pedal of an ellipse or hyperbola. 

170. Polar Reciprocals, 

The inverse of the pedal of a curve with respect to 
the same pole is called a polar reciprocal of the curve. 

For instance, the pedal of a circle with respect to any 
point is the limaqon (§ 1 56) 

r = a+6cos Q\ 
and therefore the polar reciprocal of a circle is the inverse 
of a lima^on, and its equation is 

r=c7(a+6cos 0), 
the polar equation of a conic section, with a focus at the 
pole, excentricity e = b/a, and semi-latus-rectum l = €^/a. 

Again, the polar reciprocal of an epicycloid with respect 
to the centre is r cos m6 = b, a Cotes's spiral (§ 164). 

With the usual notation of 0, r, 0, p, p, ,,, for the 
original curve, and ffy /, 0', . . . for the polar reciprocal 
with respect to the origin, then from §§ 155, 168, 

+ 0' = 7r, / = c2/^,p' = c7r; 
*°^ ff^'^^dp''^ dp' ^^^^ cosecs^. 
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171. Orthogonal and Oblique Trajectories. 

A curve cutting at right angles a system of curves is 
called an orthogonal trajectory of the system ; and a 
curve cutting the system at any constant angle y other 
than a right angle is called an oblique trajectoi^y. 

In § 101, it was shown that the orthogonal trajectories 
of the cycloids described by all the points on the circum- 
ference of a wheel are equal cycloids. 

Familiar instances of orthogonal trajectories are seen 
with (i.) horizontal and vertical straight lines, (ii.) straight 
lines radiating from a centre, and the concentric circles, 
(iii.) confocal and coaxial parabolas, (iv.) rectangular 
hyperbolas, the asymptotes of one system being the axes 
of the other system, (v.) confocal ellipses and hyperbolas ; 
while the oblique trajectories of system (i.) are inclined 
straight lines, of (ii.) are equiangular spirals, of (iii.) are 
confocal but not coaxial parabolas, and of (iv.) are con- 
centric, but not coaxial rectangular hyperbolas. 

Since any two curves and the corresponding inverse 
curves cut at the same angle, it follows that the inverse 
of a system of orthogonal or oblique trajectories forms a 
new system of orthogonal or oblique trajectories. 

We have seen (§ 133) that the graphs of the conjugate 
functions u and v, given by u+iv = {(x+iy), when u or 
V is equated to a constant, form an orthogonal system. 

Also, since in an oblique trajectorj'^ cutting the curves 
u and V (§ 133) at constant angles y and iTr — y 

cfojCOS y + cZsgSin y = 0, 
and dsjdu = dsjdv = J-i; 

therefore du cos y+dv sin y = 0, 

or u cos y+ V sin y = a constant, 

is the general equation of an oblique trajectory. 
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Mr. Larmor has shown (Proc, London Math, Sac, vol. 
XV.) that in a field of force (§ 83), of which the potential 
is J, a particle will describe an oblique trajectory with 
velocity F, if started so that JF^ = ./; or that a ray of 
light will describe this oblique trajectory if the refractive 
index varies as t/"^. 

Then, if It denotes the radius of curvature of the 
trajectory, resolving normally, 
V^ '2J dJ .3 J . ji/dJ . 3 J . \ 

1 3 J4 . 3 J4 . 

or -n=:^— cosy+-,^— smy. 

Any small contour in the (cc, y) diagram will be 
transformed into a similar contour in the (u, v) diagram, 
since angles are unaltered in this transformation, J^ 
denoting the scale of the transformation, or the ratio of 
the linear dimensions of corresponding elements ; this 
is the condition to be satisfied in maps and charts. 

Thus, for instance, if 

xi-\-iv={x-\- iyY = 9**cos nQ + i7'*^sin nQy 
by De Moivre's Theorem (§ 111), then, on putting 

u = c'*^-W = r^cos nO, v = c^~^y' = r'^sin nO ; 
and a new system of curves is obtained in which r is 
changed into r^/c^"^ and 6 into n0; so that we may put 

and now the radial angle <f> is the same in the old and the 

new system ; and 

r r dr'sind) dr sin d) , -w - , 

n—,—-=n — j-r-^ ^ -^ = (?i — l)sin0. 

p p dr dr ^ 

The case of n = — 1 is equivalent to the operation of 

inversion {§ 168) combined with reflexion in Ox. 
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Fig.53 

172. Critical Orbits, 

The orthogonal trajectories given by 

u = r'^cos nQy v = r^sin nO, 
are of important interest ; the oblique trajectories being 
u cos y+v sin y = r^cos(n0 -- y) = a constant. 

These curves are called CHticcd Orbits by Clifford 
(KinematiCy p. 113) because they are described as orbits 
under a central force varying as some power, 2n — S, 
of r (§ 84), while the velocity varies as some other power, 
71 — 1, of r. They are also catenaries, the curves assumed 
by a uniform chain under a central force varying as the 
power n — 2 of r. 

Writing the equations, 

r^cos nO = a**, r**sin nO = &**, r^eos(n6 — y) = c** ; 
then (i.) n= ±1 gives a system of straight lines and 
circles, (ii.) n= ±^ a system of confocal parabolas and 
cardioids, (iii.) 'Mss ±2 a system of rectangular hyperbolas 
and lemniscates. 
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Changing n into -n gives the inverse sTstem of curves 
r"=a"cos7i6, »" = 6"sinn^, r''=c'\»s(n0 — y); 
and now, in r"= a"cos nQ, 

cot = d log rjdQ = — tan nQ = cot(^7r + n0), 
or ^ = Jx+Jnfl; ao that (fig. 53) 

POY=ne,&xiA AOY=w={n+\)e. 
Then •p = r cos w6 = 7"+ '/a", 

or p=a(cosn0)i+>'-=a(coa ^^)^'. 

p"=a''VMsmUi», where m = n/(n+l); 
so that the pedal curve AY oi the critical orbit ^P is 
also a critical orbit ; and the equation of the polar recip- 
rocal of AP may be written 

r*cosm0=a'", 
- another critical orbit 

Consider the path of a ray of %ht in the Earth's 
atmosphere, on Simpson's assumption that the refraction 
.index n varies inversely as the n+l"" power of the 
Tfloce r from the centre of the Earth. 

u'^2>.^sin0 or fip is constant along the ray, or 
by De Moivre's Theo^i the path of the ray is either the 
u,+iv==c"' »" oblique trajectory, 
w = c" - ix' = r"co8 710, w = c^ of curvature through 
and a new system of curves is obtl)> 
i into !■"/(;"" ^ and Q into nB; sJ.*- 
/=.r"/c»-i, e'=TO0; 
e radial angle is the same in ip pole P of the 
; and he straight 

n= — 1 is equivalent to the operation of 
) combined with -refierion in Ox. 
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Differentiating l<^rithmi«allT nilb respect to a, 

m dy , d^ 

— -ji!=-ianiW-; 
y as '^ tU 

, , du . . (i« 

and since -^=-~smvr, -r-r=0, 

as "^ djfr "^ 

therefore /> = y sec ^. m = PJ m, 

so that p the radios of carvature of the roulette U 1 m 
or 1 + 1/n times the length of the normal PI, and there- 
fore m/(m+ 1) times the radius of curvature of the rolling 
pedal A'P ; and the evolute of the roulette will possess 
a similar property. 

For instance, if i( = l, the rolling curve is a circle, and 
the roulette is a cycloid ; and then p=2PI (§ 101), and 
the evolute of AP is an equal cycloid. 

As applications of these roulettes AP, we may instance 

(L) that AP is a catenary curve with Ox horizontal, 
when the line density is made proportional to the 
2 + 1/71"' power of sec^, or of the tension; reducing, 
tor n= —^, to the ordinary catenary, which is therefore 
the roulette of the focus of a parabola ; 

(ii.) that, with m negative and=— ^, AP can be 
described as a trajectory by a projectile in a resisting 
medium, in which the retaliation due to the resistance is 

where v denotes the velocity at P, and v^a,t A ; reducing, 
for p = ^,to the ordinary parabolic trajectory, the roulette 
of the pole of r*co8i0 = o*, the 6rBt n^atiye pedal 0' 
parabola r^cos Jd=a* ; 

(iii.) that AP is also the path of a ray of li^t, 1 
the refractive index varies aay^; a caienaiy foi; f 



reducing, 
e roulette^ 
dalDfl^^ 
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174. Vectors and Vectorial EqvMiona. 

The quantity z=x+iy i^ now called the vector or step 
from to the point (oj, y) ; being compased of a step of 
length X parallel to Ox and then a step of length y parallel 
to Oyy according to Argand's representation. 

Then if %'^x'+iy' represents another vector, 

z—z'=x—x'+i. y — y' 
represents the vector or step from (a?', y') to (ar, y\ 

It is convenient in physical application to consider the 
logarithm of the vector, and now 

u»=log(0— s?') 
or it+i'y=log(ic— a;'+i. y — y') 

= logV{(^-^T + (^-2/T}+itan{(2/-y')/(^-^0} 
gives u the velocity function, and v the stream function 

at {x, y) of a source or electrode at {x, y') ; and we may 
call w the vector potential. 

Suppose for instance we place n equal positive elec- 
trodes at points A^, A^ ... ^n-i» equally spaced on a 
circle of radius a ; the vector of the point -4, being 

a cos 287r/n + ia sin 2s7r/n = a exp(2is7r/ri ), 
the vector potential 

8=n-l 

IV = log n{0 — a exp(2is7r//i) } = logf^" — a**) ; 

« = 

and now, with polar coordinates, 
z =r cos 0+ir sin = re^^, 
z^ = r V"^ = r^cos nQ + ir'^sin nQ ; 
w = log(r'*cos 710 — a** + tr^sin n&), 
u + iv = logy/(r^^ — 2a^r~cos nO + a^) 

+ i tan " ^{r^in n6/(r^cos nO — a**)} ; 
or r^** — 2a^r^cos 7i6 + a^ = e^, 

?^^sin('y -- 7? 0) = a**sin v, 
representing an orthogonal system of curves. 
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Denoting by r„ 0, the polar coordinates of a point P 
relatively to the origin -4,, then 

u=log r^r^r^ ... rn_i, t;=0o+^i + 02+ ••• +^n-i. 
giving De Moivre's and Cotes's properties of the circle. 
Suppose for instance n = 2, then 

u = log r^r^y or ?•* — 2a^r^cos 20 + a* = e^" ; 
and V = 00 + Oj, or r2sin(20 — v) = — a%in v, 

representing Cassinians, and the orthogonal rectangular 
hyperbolas, passing through -4^ and A^ 

When n equal negative electrodes are placed on the 
circle, at 5^, B^, ... 5n» midway between the positive 
electrodes, the vector of 5, being 

aexp{(28— l)7r/ii}, 
then their vector potential is 

log n {0-a exp(28- l)7r/'yi} =log(2?^ + a**) ; 

so that for the system of positive and negative electrodes 

ty = log(0^ — a«)/(2;*» + a") 

— 11 y^^ — 2a^r'^cos nO + g^^ 
""2 ^S rfin ^ 2a~r"cos 7i0 + a^" 

.. ., r^sin-n^ .. , r**sinti0 
+ ttan ^— T. r— ttan-^ 



r^cos 710 — a" r**cos 7i0 + a** 

. , ^ia'^r^cos 710 .^ -2a^r^in'}70 

giving another orthogonal system of curves, 

9'2n + 2a^7'~cos 710 coth u + a^" = 0, 

r2» + 2a"r'*sin 7i0 cot t; — a^** = 0. 

When 71 = 1, these curves are circles, the dipolar circles 
of the next article. 
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175. Dipolar Circles, 

When we put x+iy=^ciei,Ti^{u+iv), and write it 

. _ sin ^{u+ iv)QOs\ {n'- iv) _ sin ii+i sinh v 
^ "" cos \{u + iv)cos J(u — iv) ~" cos u + cosh v 

^, c sin u c sinh v 

then a; = — r — ^ > 2/ = — r — : > 

cosh V + cos u ^ cosh v + cos u 

thus dividing taji^{u+iv) into a real and imaginary part. 
Also tan u = tan { J(u +iv)+^(u-^ iv) } 

__ x+iy+x — iy ^ 2cx 
"" c^'-x^—y^ " c^—x^—y^' 

and tan iv = tan{ J(i6 + iv) — J(u — iv)} 

_ x+iy—x+iy __ 2icy 
■"^ c2+aj2-f2/^ "c^+x^+y^' 

so that a;2+y2+2cflj cot it — 0^=0 (i.), 

x^+y^ — 2cycoihv+c? = (ii.), 

the equations of a system of orthogonal dipolar circles 
(fig. 54. i.) ; called dij9oZar because the system (i.) represents 
circles passing through two poles N and S, while system 
(ii.) represents orthogonal circles having the same radical 
axis Ox ; thus representing the meridians and parallels in 
the stereographic projection of the terrestrial globe. 

Putting u or v = nt, we obtain circles described under 
excentric centres of force ; and writing them 

(a;+ccotu)^+y^ =c^cosec% (i.) 

aj^+(2/ — ccoth v)^ = c2cosech^i; (ii.), 

then (i.) represents a circle, centre ( — c cot u, 0), and radius 
c cosec u ; and (ii.) a circle, centre (0, c coth v) and radius 
c cosech V ; and with respect to 0, the power of (i.) is — c*, 
and of (ii.) is c\ 

At any point P, the angle NP8=Tr — n, 
while the ratio SPjNP^e^ or v=\og(SP/NP). 
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If X denotes the latitude of any parallel of latitude 
circle (ii.), then since on the stereographic projection the 
radius of this circle is c cot X, therefore (§ 34) 

cot X = cosech v, tan X = sinh v, sin X = tanh v, 

cos X = sech V, tan JX = tanh Ji;, or X = gd v, 

V = log (sec X + tan X) = log tan(| x + JX). 

while X = gd(ucoty) represents a loxodroTne or rhuTnh 
linCj cutting the meridians at a constant angle y. 

The inverse system of curves with respect to a pole N 
or 8 will be radiating straight lines and concentric circles, 
representing meridians and parallels in the neighbourhood 
of a pole, when projected stereographically from the other 
pole, as in Godfray's Great Circle Chart. 

A great circle on this chart will be represented by a 
straight line ; and a rhumb line, that is an oblique 
trajectory of the meridians and parallels, will be an equi- 
angular spiral ; so that on the stereographic projection 
it will be an inverse of an equiangular spiral. 

The general operation of inversion with respect to any 
pole, combined with a change of origin, is represented by 
the linear substitution (Schwarz) 

z' = (az+b)/(Az+B), 

by which a polygon in the z plane, bounded by straight 
lines or circles, is transformed into a polygon bounded 
by circles, cutting at the same angles, in the z' plane. 
In Godfray's Great Circle Chart we transform to 

z' = x'+iy'=:cexp{ia — /3) = ce-^(coaa+i8ina), 
so that a is the longitude, and the latitude X is given by 

r' = ce-^ = c tan(i7r — ^X), 
or ^ = log(sec X + tan X), X = gd ^8. 
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176. Mercatcyt^'a Chart, 

Draw any contour on the (u, v) diagram (fig. 54, ii.) 
corresponding to a contour, the outline of a country, in 
the (x, y) stereographic diagram ; we thus obtain a new 
map, called Mercator's projection, in which the meridians 
and parallels are orthogonal systems of parallel straight 
lines ; and while a length u represents the longitude, a 
length v = log(secX+tanX) will represent the latitude X; 
also cZ'y/cZX = secX, so that the minute of latitude in Mer- 
cator's chart increases as the secant of the latitude, being 
equal to the minute of longitude only at the equator. 

Twenty-four standard meridians at equal intervals of 
15° in longitude from Greenwich mark the standard 
time at a place ; the standard time being the mean solar 
time at the nearest standard meridian. 

Taking the Godfray Chart as representing any system 
of parallels and meridians with respect to an axis through 
any zenith Z, say in longitude and latitude gd 6, then 
the representation on Mercator's Chart will be given by 

taxi^{u+i , v~'b) = exip (ia— /3) = x' +iy' =z' ; 

or, dropping b as representing a mere change of origin, 

ia -- /8 = log tan J(u + i'y) 

, sinu+isinhv ., .sinhv ., ,coshv 

= loe: ; r — =^tan"•^-^-T coth"-^ ; 

° cos u + cosh V sm u cos u 

so that sinh v = tan a sin u, cosh v = coth /3 cos u ; 

the representation on Mercator's Chart of a system of 

parallel small circles and their meridians ; these lines are 

called Sumner lines (p. 203), being the lines on Mercator's 

Chart on which a celestial body in the zenith at Z will at 

any instant have the same altitude. 

2a 
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Fig. 66 

177. Confocal Ellipses and Hyperbolas, 

As another practical illustration of conjugate functions, 
put x+iy^c cosh (u + iv) 

= c(cosh u cos v+i sinh n sin v) ; 
then a; = c cosh u cos Vy 2/ = c sinh u sin t; ; 

and alternately eliminating u and 'y, 



X' 



r 



c^cosh^ie, c^sinh^u 



,(i.) 



x^ 



y 



c^cos^-y c^sin^v 



= 1. 



.(ii.) 



representing a system of confocal ellipses and hyperbolas 
for constant values of u and v the polar reciprocals with 
respect io N ov 8 of the circles i; of § 175; and sech u will 
be the excentricity of the ellipse, sec v of the hyperbola 
(fig. 55); while (§ 35) v will be the excentric angle or 
excentric anomaly of a point P on the ellipse AP, and u 
will be hyperbolic excentric anomaly of the point P on 
the hyperbola KP, 
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Put v^nty then 

x+iy = ccos\i {u+int) = ccoB (nt—iu) 

represents confocal ellipses described under a central 

attraction to 0, varying as the distance, the period being 

2'7r/n; while hyperbolas are described with u = nt, and 

x+iy = c cosh (nt + iv). 

In an oblique trajectory 

u cos y+v sin y = constant, 

we may put v = 7i{t—T)cosy, u= --'n(^ — T)siny; 

so that x+iy==ccosn{(t'-T)co8y+i(t'-T)smy}; 

and the equation is obtained by eliminating t between 

x= c cos {7i(^ — t)cos y } cosh {n(t -- T')sin y } , 

y= — c sin{'w(^ — t)cos y}sinh{7i(^— T)siny}. 

The equations of the tangent and normal of the ellipse 

at P, which are also the normal and tangent of the 

hyperbola at P, are 

cosv . sinv ..... 

X — r — h y-r-T — = c ( m. ), 

coshu "^smhu 

coshtt sinhu .. . 

X y—. =c (iv.). 

cosv ^ smv 

Denoting by p and q the lengths of the perpendiculars 

from upon the tangents at P of the ellipse and hyperbola, 

p = c cosh u sinh u/(cosh% — cos^)^, 

q = c sin v cos v /(coshhi — cos^i;)^ ; 

The pole of (iii.) with respect to the hyperbola (ii.) will 

be (c cos^v/cosh t6, — csin^v/sinhu), and this pole will 

therefore be Q, the centre of curvature of the ellipse at 

P (§ 98) ; and PQ = p = c(cosh% -- cos^i;)*/(cosh u sinh u). 

Similarly (c cosh%/cos v, — c sinh%/sin v), the pole of 

(iv.) with respect to (i.), is the centre of curvature of the 

hyperbola at P ; and the radius of curvature is 

c(cosh^ — cos^)*/(cos V sin v). 
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We may consider w=u+iv b& representing a vector of 

curvilinear steps n and v\ u along a hyperbola and v 

along a confocal ellipse, meeting at the point 

z=x+iy = ccosh(u+iv). 

If P' is any other point given by 

z=x+ iy' = c cosh(t6' + iv'), 

then the vector P'P is given by 

0— 0'=c{cosh(t6+iv) — cosh(u'+iv')} 

= 2csinh i^{u+u'+i.v+v')sinh J(t6— u'+i. v— v'); 

so that PP^ 

= 4c%inh J(u+u'+i . i;+t/)sinh J(u--u'+i . v— v') 

sinh J(w+u'— i.i;+r')sinh|(u— u'— i.t;— v') 

= c2{cosh(u+u')— cos('y+v')}{cosh(u— u')— cos(t;— v')}. 

This expression is unaltered by an interchange of t«, < 

or V, v' ; so that if the points p, p\ called corresponding 

points, are given by (u, v'\ {u\ v), then pp' = PP\ 

Further, if denotes the angle between PP" and pp\ 

,, . , cosh(t6+iv) — cosh^u'+ivO 

then ^0 = log — tt — r-^-4^ ,;/../ 

^ ° cosh(i6 + tv ) ■- cosh(u + ^v) 

__ , tanh ^(i^ — u') + i tan |(v — v') 

"" ° tanh J(u — u') — i tan ^(v — i;') 

tanh J(u— u') 

Captain Weir's Azimuth Diagram consists of a series 
of confocals, the hyperbolas marking the longitude or 
hour angle v, while the ellipses are marked with the 
degrees in gd u, to represent latitude or declination. 

Now, if the hour angle is the complement of v, and if 
the latitude and declination are gd u and gd cZ, a straight 
line joining the points (t6, v) and (d, 0) will give the 
azimuth or bearing of the object from the meridian Oy, 

(Godfray, Astronomy^ § 222.) 
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Changing to the focus £f or /S' as origin 

a; — c + iy = 2c sinh^J(u + iv), x+c+iy = 2c cosh^ J(u + iv), 

so that if r, r' denote the focal distances /SP, S'P, 

r = 2c sinh \{u + i'y)sinh \(u — iv) 

= c(cosh u -- cos v) = a(l — c cos v) ; 

r'= c(cosht6+cos'y) = a(l+ccost;). 

Then r'+r = 2ccoshu = -4J.' (v.), 

7''— r=2c cos v = KK' (vi.); 

and hence 8P=^AK, S'P = KA'; also OP=KB, 

Denoting the angles ASP, AS'P by 6, ff, then in 

Astronomy Q is called the ^ruc anomaly of P, reckoned 

from peHhelion at J., the Sun being supposed to be at 

the focus S ; but ff is the true anomaly from aphelion A, 

with the Sun now at the other focus S\ 

Then 

^ x — c cosh u cos v—1 ^, cosh t6 cos v+1 

COS0= = r , cosy = r ; ; 

r coshu— cost; coshu+cosr 

. ^ sinh u sin V . ^ sinh it sin v 

or sin0= — i , sm0^= 



cosh u — cos V "" cosh u + cos v ' 

tan ^6 = coth ^u tan ^v, tan ^ff = tanh |u tan |v ; 

and r = >Sil cos^v sec^O, r' = /Sf J-'cos^Jt; sec^^Q'. 

. . . l+cos0cosh2>6 e+cos6 

Again, since cos v = r — ; tt = -r-. 7\> 

° coshtt+cos0 1+CCOS0 

or (1 + c cos 0)( 1 — e cos v) = l — e^; 

therefore r = c sinh^i6/(cosh u + cos 6), 

of the form r = l/{l + e cos 6), with l = c sinh%/cosh u, 

the polar equation of the ellipse, with origin at 8. 

Similarly / = c sinh^u/(cosh u — cos 6") 

is the polar equation of the ellipse, with origin at S' ; 

while for the hyperbola the corresponding polar equations 

c sin^t; , c sin^t; 

are r = t,, r = 



cos V — cos 6' cos ^ — cos V 
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178. Confocal Limagons. 

By inversion of this system of conies with respect to 
either focus, say S, taking 2c as the constant of inversion, 

x+iy = 2c cosech^Ku + iv), 

x+2c+iy = 2c coth^^(u + iv) ; 

and rsinh% = 4c(cosht6+cos0) (vii.)> 

r sin^v = 4ic( cosr— cos0) (viii.), 

are the inverses of the ellipse (i.) and hyperbola (ii.), and 
the pedals with respect to N or S of the circles v oi% 175, 
a system of confocal lima9ons (§ 155). 

Denoting now by r' the distance from the other focus 
S\ the position of which is unaltered by inversion, then 

/ o xi> 1/ I • \ j.1. 1/ 'NO coshu+cosv 

r = 2c coth i(u+ IV) coth hiu — ^V) = 2c — , , 

^"^ ^ j\ / cosh u — cos V 

while 

r = 2c cosech Uu + i'v)cosech Uu — iv) = — r ; 

^^ ^ ^^ ^ cosh u— cos -v 

so that r'+2c= — r r = r cosh u (ix.), 

cosh 16— cos v ^ ^ 

, ^ 4c cos v , . 

r — 2c= — ,- = r cos'y (x.), 

cosh u — cos v ^ "^ 

equations of the form r' — Ir = constant, which are 
Cartesian ovals in the general case, the inverse of a 
conic with respect to any point on the transverse axis. 

The lima9on in which t;= Jtt, Z= J, is called the Trisectrix. 

For if A' OR is any angle to be trisected (fig. 47, ii.), and 
if A'R is joined cutting the curve in p; then, if A'0p = 6, 

Op = OA\2cose-l); and A'p = 20 A'siniO^A'q, 

if Op meets the arc A'R in q. 

Therefore e = qpA- OA'p^pqA^-^OA'p^OA^q-OA'p 
= 'pA'q = \qOR\ so that Op is a trisector of the angle A' OR, 
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179. Kepler* 8 Laws of Planetary Motion, 
Kepler, by long continued observations and measure- 
ment of the Sun's diameter and motion in longitude, 
noticed that 

(i.) the variation of the Sun's apparent diameter d 
could be expressed by the formula i)(l+ecos0), 
where D denotes the mean angular diameter 
(about 82') and Q the Sun's longitude from peri- 
helion, e being a small constant, about 1/60 ; 
(ii.) that the Sun's daily motion in longitude was 
proportional to the apparent area or square of the 
diameter. 
Since the apparent diameter is inversely proportional to 
the distance, he deduced the laws called Kepler's Laws — 
(i.) that the relative orbit of the Earth (or any 
planet) and of the Sun is an ellipse, with a focus 
at the Sun, if the Sun is supposed fixed, given in 
polar coordinates by Ijr = d/D = l+e cos 6. 
(ii.) that rH6/dt = h is constant, and that the Earth 
therefore sweeps out by its radius vector from the 
Sun equal areas in equal times. 
The third law of Kepler — (iii.) the squares of the 
periodic times of the planets round the Sun are 
proportional to the cubes of their mean distances — 
was easily inferred by arithmetical calculation, when once 
the distances of the planets from the Sun were measured, 
in terms of the Sun's distance from the Earth. 

Newton inferred from law (ii.) that the earth is attracted 
by the Sun ; and from law (i.) that the attraction must be 
inversely proportional to the square of the distance: 
while law (iii.) showed that the attraction of the Sun was 
of the same nature on all the planets (§ 84). 
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By Newton's Law of Universal Gravitation, employing 
CG.s. nnits (§ 142), the attraction between two spherical 
bodies, for instance, the Sun and the Earth, weighing S and 
E grammes, when their centres are a centimetres apart, 
will be given by the expression GSEa~^ {dynes)] and G 
is called the constant of gravitation^ being the attraction 
in dynes between two spheres, each weighing one gramme, 
when their centres are one centimetre apart. 

Then Kepler's Third Law, in a mathematical form, 
asserts that if T is the period in seconds of the Sun and 
planet, and if n = 27rlT denotes the mean motion, then by 
ex. 1, p. 175-, 7iV = 47r2a7r2= C(S+E). 

According to the Cavendish experiment, now being re- 
peated with improved apparatus by Mr. G. V. Boys, we 
may take (7=10-8x6-48, 1/(7= 10' X 1-54. 

(Everett, Units and Physical Constants,) 

Denoting by g the acceleration (in c.G.S. spouds) of the 
attraction of the Earth, then 

g = CE/R^, or GE==gR^; 

so that we can determine E when C is known, and vice 
versa; and till this is done, Newton s PrimcijAa is merely 
Kinematics. 

With the above value of 0, and sr = 981, i2=10»-r-|7r, 
E=gm/G= 1027 X 6-12 grammes ; 

giving a mean density /o = -£'/(^7ri2^) = 5'67. 

Denoting by 11 the Sun's parallax (8"-76), and by T 
the number of seconds in one year, then the attraction 
between the Sun and the Earth is 

F^ 4^7r^ER cosec U/T^ = 10^7 x 365 dynes, 

while S = Fa^GE = 10^3 x 1 2 grammes, 

with the above numerical values. 



CURVES IN GENERAL, 377 

180. Elliptic Planetary Motion, 

Considering an elliptic orbit AP (fig. 55) in which 
c cosh u = OA = a,c sinh u = OB = 6, c sinhHt/cosh u=l, the 
semi-latus rectum, and sech u = e ; then the sector 

ASP = sector A OP — triangle OSP = ^a6(i; — e sin t;). 

If the ellipse is described in period T round the Sun at 
S, with mean angular velocity n = 27r/T, and if the time 
to P from perihelion at A is denoted by t, then by 
Kepler's Second Law, 

^_^a6(i;— esinv) 

or nt=V'-esinv'y 

and now nt is called the mean anomaly^ being the true 
anomaly of a planet moving uniformly in a circular orbit 
in the same period T, 

Expressed in terms of the true anomaly 6 from peri- 
helion A, 

€+cos6 . ;^(1 — e^)sin0 

l + ecos0 1+6COS0 

l+ecos0 l+ecos0 

Reckoned from aphelion A with the Sun at S\ 
nf=v+esinv 

" l-ecos0' "^ l-ecos0' ' 

For completeness we may consider the hyperbolic 
branch KP, of excentricity e\ described round the Sun 
at S, and now the mean anomaly 
n't = e's\uhu''U 

_ eV(e-^- l)sin g .,, s/(e''- l)3iD 6 
~ l + e'cos0 l+e'cosd ' 

with n'a'«=G{8+E), 
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a' denoting the semi-transverse axis OK, and the true 
anomaly Q being now KSP, reckoned from perihelion K. 
With the Sun at S\ the hyperbolic branch KP will be 
described under a repulsion from S\ and 
7i,Y = e'sinhu+u 

_ gV(e^^-l)sin0\ . ^J{e2^)_^mff 
"" e'cosd'-l ■^^'''^ "^'cosd'-l • 

181. Rectilinear Motion under Varying Gravity. 

By making e = l or u = 0, the ellipse becomes a straight 
line, and we obtain the solution for rectilinear motion 
under attraction to a fixed centre, varying inversely as 
the square of the distance ; for instance, the motion of 
a body shot vertically upwards from the Earth with 
very great velocity. 

Then, if R denotes the radius of the Earth, 



d^x jB2 ,1 dx^ 

and — 



-<-rJ 



dt^ ^x^ ' 2 dt^ 

if the body reaches a distance 2a from the centre of the 
Earth ; and if E denotes the weight of the Earth, and C 
the gravitation constant, n^a^=gR^ = CE. 

Now x = a{l+cosv), where nt = v+smv, 
for an attraction ; but for a repulsion 

X = a(cosh u + 1), where nt = sinh u + u. 
If denotes the centre of the Earth and the body 
starting from A falls to P in the time ty and if PQ is 
drawn at right angles to OA to meet the circle described 
on OA as diameter in Q; then the angle J.OQ = Jv, and 
the area A OQ = ^a\v + si nv) = \a^nt ; 

so that the time from J. to P is proportional to the area 
AOQ, and the point Q moves freely on the circle -4Q, 
under an attraction to 0. 
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To shoot a body weighing W\h. vertically upwards away 
from the Earth to a distance 2a feet from the centre, say 
the distance of the Moon, requires energy, in ft. lb., 

iW/^=Tfie^Q-4)=Fii(l-|); 

and with the Moon's parallax of 57', 2a/i2 = cosec 57' = 60, 
about ; while iZ, the radius of the Earth in feet (§ 99), 

= 90x60x608p-T-Jx. 

If we could shoot a body up to a height R at the 
North Pole, then fired horizontally the body would, in 
the absence of resistance, skim the surface of the Earth 
like a grazing satellite, in a period 2Trln = 2Tr^{Rlg) 
seconds; and n^=g/R = CE/R^ = ^'7rpG. 

Assuming that gravity inside the Earth varies as the 
distance from the centre, as would be the case if the 
Earth were homogeneous, 27r,^(R/g) sec. is also the period 
of oscillation in a diametral tunnel from pole to pole. 

The period of the Moon being about 27 days, the period 
of a satellite about one-ninth the distance of the Moon 
would be, by Kepler's Third Law, about one day, and the 
satellite would appear stationary in the sky ; and this 
would make the period of a grazing satellite 27-^(60)^ or 
about one 17*^ of a day. 

In the plane of the equator such a satellite would fly 
past a point 16 or 18 times a day, every 90 or 80 minutes, 
according as it moved eastward or westward ; and the 
velocity over the ground would be 7^ or 8J kilometres 
per second. 

If the Earth was made to rotate 17 times faster, without 
change of shape, bodies at the equator would be on the 
point of flying off* into space, and elsewhere the plumb 
line would point to the Pole Star. 
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182. Kepler's Problem, 

This is the problem in Astronomy, to express the true . 
anomaly Q and the excentric anomaly v in terms of the 
mean anomaly nt ; and for this purpose we acquire the 
series of Lagrange (§ 150) and Fourier (§ 183). 

Denoting by m the mean anomaly nt, or more strictly, 
in astronomical notation, n^ + e— or, where e denotes the 
epoch, and or the longitude of perihelion, then 

v^m+e^m v, 
which, by Lagrange's Theorem (§ 150), gives in a series 
proceeding in ascending powers of e, 

, ^ e^ dP-'^ , . . 
^= "^^^ ^ di^^^'"^"^^"^ 

sm v = sm m + z,p — -=^ j — - (sm m)P^\ 

{p+l)\dmP ^ ^ 

Similarly sin 2y, sinSt;, ..., sinp;, and cosv can be 

expanded in powers of e ; and 

r ^ dm 

a dv 

a dv - ^eP dP , , . 
- = ^- = 1 + 2— , :7-^(sin m)P ; 
r dm p\ dmP ^ 

while from tan|0 = coth Ju tan Jt;, 

we deduce by logarithmic differentiation 

dQ dv dO sinhu 1 — c^ 

or 



sin sinv* dv coshu — cosv 1 — 2ccos'y+c^' 
with c= e~^ ; and, resolved into partial fractions, 

de^ 1 1 

= 1 + 22c^cos pv, 

qp 
6 = v + 22— sin pv, 
p ^ 

which can now be expressed in terms of m. 
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We can rearrange these series so as tx) proceed accord- 
. ing to cosines or sines of multiples of the mean anomaly 
TTiy the coeiBcients being functions of e, and this is more 
convenient for Kepler's Problem. 

A series of this nature is called a Fourier Series, and 
we proceed to show how the coefficients of such a series 
can be calculated, for any arbitrary single valued function 
fflj, in a series proceeding by cosines or sines of multiples 
of xjl, when I is any arbitrary quantity. 

183. Fourier' 8 Series, 

Assume that, between the limits x= ±1, the function 
tx can be expressed by a series of the form 

ix — ^A^ + 2ilj>cos(297raj/i) + 25psin(p7ra?/i), 
where 2 denotes the sum of the series obtained by giving 
p all positive integral values from 1 to oo ; it is required 
to determine the -4's and J5's. 

Divide the function ix into its odd and even part (§ 46), 
and denote them by f^aj and i^ respectively ; thus 

{^ = Wi^) + f( "■ ^)} = i^o + 2^i>cos(p7raj/Q (i.), 

fjcc = \{i{x) — f( — aj)} = 2J5j,sin(p7ra?/i) (ii.). 

To determine Ap, change x into v in (i.), multiply both 
sides by cos(p'7rv/l), and integrate with respect to v be- 
tween the limits and I ; then since 

/ cos 2j- cos ^ dv = / \ ^ cos(p - ?)^+|cos(p+g)-j- \dv 



_ rsm{(p-q)Trv/l} sin{(p + q)'7rv/l} ']^ 
L 2{p'-q)irll 2(p + q)7r/l X ' 

while ycos^r =/\h + h cos ^-^)dv = il ; 



therefore /f^v coB{p'jrv/l)dv = ^lAp ; /f^vdv = ^lA q. 
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Similarly, changing x into v in (ii.), multiplying by 
sin (pirv/l), and integrating with respect to v between 
and I, we find 



/f-^v ain(p7rv/r)dv= ^IBp. 



Therefore, between the limits x= ±1, 



=5/ 



^ ^ ^ "• cos ^-^r-dv 



fo = T /fgt^cZv + yScOS , 



y?:. 



z 



+ -tS s\n^-j—/{^v sin -^^dv, 



and this is called Fourier's Series. 

By supposing p to take all integral values, positive or 
negative, we may write Fourier's series, 

faj = 5yV /ft; cos ^ , ^c?i;. 
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At the limits x=±l, the value of the series is fgZ, 
so that in general there is discontinuity at the limits; 
and outside these limits Fourier's series represents periodic 
repetitions of the function fee between the limits. 

This is exhibited by drawing the graphs of 

(i.) ^Aq + ^ApCos{p7rx/l) and (ii.) XBpSm(p7rx/l), 
as exhibited in fig. 56, (i.) and (ii.) 

When the limits between which faj is to be expressed 
by Fourier's Series are and I, then either series (i.) or 
(ii.) may be chosen at pleasure, to represent ix, but it is 
best to choose the series which introduces the least dis- 
continuity at the limits. 

For instance, suppose fx = x; then between and I, 



or 



(ii.)x= -L slul^^^in^^ ; 



but outside these limits the series (i.) and (ii.) represent 
the dotted lines in fig. 57, (i.) and (ii.). 
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The student is recommended to draw the graphs of 
the first two or three terms in a Fourier series, to see how 
quickly the series approximates to the given function. 
For instance, draw the graphs of 

sinaj — |sin 2cc, sinaj — Jsin 2x+ \ sinSa;, 

cos aj + 1" cos 3x, cos aj + ^ cos 3cc + -jt cos 5aj, 

Fourier made great use of his series in problems on 
the Conduction of Heat, where a solution of the partial 

differential equation — = A:;^ is required, giving the 

temperature u at any time t\ for if the initial state of 
temperature TJ is expressed by the Fourier Series 

Cr= I J[^+ 2(.4pCOs j%c 4--BpSin_2:wj), 
then at any subsequent time t the temperature 

u — \A^-^ 2(^j>cos yx + JS^sin pic)e -^^**. 

The geometrical meaning of the coefiicients -4^ and J5p 
is explained by Clifford {Proc. London Math. Soc, vol. v.) 
according to the method employed in Tidal Harmonic 
Analysis (G. H. Darwin, BHt Ass, Report, 1883), where 
a pencil attached to a float registers tx, the rise and fall 
of the tide, on a cylinder turning uniformly by clockwork 
round a vertical axis. 

To analyse the tides of period P, the cylinder is made 
to revolve in a period P/p, when a closed curve will be 
traced on the cylinder after p revolutions. 

The plane on which the projection of this closed curve 
has a maximum area A, with attention to the sign of the 
area (§ 137), will cut the cylinder in an ellipse of area A/p. 

Now, if the pencil is made to follow this ellipse as the 
cylinder revolves in the period P/p, the pencil will have 
the component harmonic travel 

ApC0s(p7rx/l) + BpSm(pTrx/l). 
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184. Fourier's Series in Kepler's Problem, 

Given the relation connecting the mean and excentric 

anomalies, v = m + e sin v, 

. , dv 1 1 + e cos a - . ^ , 

then ^— = = — ~- — s — = - = 1 4-2-d.»cos pwi, 

dm 1 — ecosv 1 — e^ r f r ' 

suppose, when expressed ty a Fourier series ; and now 
^Ap = — I ^ — QO^'pmdm,=— /cosp(v — esinv)dv; 



and this definite integral is a function of jpe, called BesseVs 
function of the order p, and denoted by Jp{pe); and it 
can be verified, as in § 207, that i; = /^(gr) satisfies Bessel's 
differential equation on p. 185 (§ 88). 

Now - = -^— = 1 + 2XJJpe)cos pm ; 

r dm 

and integrating with respect to m, 

v = m + 2i:^^£^ampm, 

P 
giving the excentric anomaly v in terms of the mean 

anomaly m. Again 

and integrating with respect to m, 

cos v=G-\- 22 ^^^ ^ cos pm. 

P 
where 

TrC=/cos vdm = /cos v(l — e cos v)dv = — ^Tre, (7= — Je; 



and sm v = = 22- sm x>m : 

e pe ^ 

whence a; = a cos v, 2/ = 6 sin v, and r = a(l — e cos 'v) 

can be expressed by a Fourier series in terms of the mean 

anomaly m. . 

2b 
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ExaTTiples on Fourier Series. 

(1) Prove that, in a series of cosines, between ± ^l, 

Prove that, between and tt, 

sinaj sin 3a; sin 5a; ^sin(2jp — l)a; i_o / \ 

-p~+ 33 "♦" 53 +••• ^^ ^ — (2»-l)3 =4'^^v'?r-a;); 

and give the value of the series outside these limits. 

(2) In a series of cosines, 

cosh mo! = 2x sinh mij^^ + ^ mfcyV ^"'^}' 

(3) In a series of sines, 

sinh vix = 27r sinh ml2 \j.y 9 a P sin^-^— . 

Deduce the expansions of 6"** and e"*^. 

(4) Expand cos mo; in a series of cosines, and sin rax in a 

series of sines, of multiples of irxjl. 

(5) Show that the Fourier series which represents the 

density or temperature in an endless wire must in 
general contain both sines and cosines. 

Determine the Fourier series to represent the 
temperature of a circular wire, in which the temper- 
atures of the four quadrants in order are 1, 2, 3, 4. 

(6) Prove that the equation 

represents a staircase, of vertical and horizontal 
steps of length Z. 

(7) Prove that the equation 

I4 = 2;^^^cos|p(a; + 2/)cos Jp(a;- y) 
represents circles of radius tt ; and draw them. 



CHAPTER VII. 

INTEGRATION IN GENERAL. 

185. In the preceding chapters a sketch of Integration 
has been given, and then a number of applications, in- 
tended to show the practical use of the method ; and now 
it is proposed to resume the consideration of Integration 
from a more general and systematic standpoint, so that 
the student may more readily perceive and write down 
the required result. 

The process of Integration is necessarily of a tentative 
nature, depending on a previous knowledge of Diflferen- 
tiation ; and in general the most convenient order of the 
mental operations employed in the integration of a given 
function will be found to be : 

(i.) to guess the function required for the integration; 

(ii.) to assign the argument of this function ; 

(iii.) to write down the proper constant or numerical 
factors of the integrated function. 

Of these three operations, the first is of the most 

fundamental importance, depending as it does on the 

principles of the Calculus, but it is the second operation 

which presents the greatest practical difiiculty, while the 

third only requires verification by a mental differentiation. 

387 
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186. General Integration of Algebraical Functions, 
The most general algebraical function of x which is 

capable of integration by means of the preceding stock 
of functions can be written 

S+TJR 

u+ vjb: 

where 8, T, U, V are rational integral algebraical func- 
tions of X, and iJ is a linear or quadratic function of x, 
and therefore of the form ax^+2bx + c. 

If JR is of the third or fourth degree in x, elliptic 
functions are in general required for the integration ; 
and if R is of the fifth or higher degree, hypereUiptic 
functions are required. 

We first rationalize the denominator, when 

8+T^/R _ (S+TJR)(U- VJR) M N 1 
U+ VJR~ [72- yiR "D^D JR' 

where D= m- V^R, M=SU-TVR, N={TU--SV)R; 

and D, M, N are thus rational integral functions of x. 

To integrate the rational function M/D, this function 
is split up into its quotient and partial fractions, in the 
manner explained in treatises on Algebra ; and then the 
integration of each term is in general easily effected. 

To integrate the irrational part, N/DJR, we may 
resolve the rational function N/D into its partial fractions, 
and integrate each term by appropriate substitutions. 

187. Integration of Rational Algebraical Functions. 
To integrate any rational function MjD or tx/Fx, the 

numerator and denominator of which are rational 
integral algebraical functions of x (§ 15), and therefore 
of the form 
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Yx" Ax''+Bx''-^+Cx''-^+.., 

where m and n are positive integers, the function is first 
resolved into its partial fractions by the ordinary rules 
of Algebra (Smith, Algebra, chap, xxiii. ; Hall and 
Knight, Higher Algebra, chap, xxiii.); if the degree m of 
the numerator is equal to or greater than the degree n of 
the denominator, the quotient must be first obtained by 
division. 

Thus to integrate ^^.^ ^ f we must suppose it re- 
solved into the form 

x + ^ being the quotient and Al(x^l), B/(X''2) the 
partial fractions of the remainder. 

To determine the numerator A, multiply both sides of 
the identity by its denominator a; — 1 ; then 

and now put aj — 1 = 0, in the identity : then -4 = — 1. 

Similarly to determine B, multiply both sides by its 
denominator x-2, and then put this denominator x - 2=0; 
this gives B = 8. 

= |ic2 +3x- \og{x - 1) + 81og(a; - 2), 
and the integration is effected. 

Every rational integral algebraical function of x, such 
as FXy can be resolved into real linear factors (factors of 
the first degree in x\ or real quadratic factors (factors of 
the second degree) ; thus, for example, 
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a;2^1 = (aj-l)(x+l); 

aj2+l is not decomposable into real linear factors ; 

aj3-l = (aj-l)(a;2+a;+l); 

a?+l = {x+l){x^-x+l)\ 

a;^ - 1 = (a: - l)(a; + l)(aj2 + 1) ; 

x^ + l^ix^-J^x+l^x^ + J^x+l)] 
ic«-l = (a;-l)(aj+l)(aj2-aj+l)(aj2+aj+l); 
a^+l = (x2 + l)(cc2-^3a;+l)(a;2 + >v/3x+l) ; 
x^^\=^{x^l){x + \){x^+l)(x^- Ji>x+\){x^+ j2x+\)\ 
X^+I^{x^'-J{^ + J2)x+I}{x^+J{2 + J2)x+I}x 

and so on. 

Corresponding to a quadratic factor in the denominator 
Yx we must assume a partial fraction with a numerator 
of the form Hx + K; thus, as resolved into partial 

fractions, we must put -5 — 7 = =-+ -5-; -r. 

^ xr-^l 03—1 x^+x+1 

The numerator A is then determined as before by 
multiplying both sides by cc — 1 and then putting a;-l = 0; 
thus -4 = J : and now, by transposition, 

Hx+K ^ JL 1 J^ -x-2 

x^+x + 1 x^-l Sx-i^Six^+x+l) 
(sothatif=-J,Z=-f). 

T.^ fdx f/l 1 1 x + 2 \, 

^^^ 7x3^=7 Vs^i-s^-qi^Tl)^ 

1, . _. 1 /'2aj+l , \ r dx 



^+x+l 

= 3log(i»--l)-^log(»Ha;+l)----^tan-i-^±l 
1, (a;-l)3 1 ^ ,2a; + l 
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Examples. — Eesolve into partial fractions and integrate 

/•t\ -L) Xy oO f V& i.f Xj X f X l.f Xf X f • * t 

^ '' (a!-l)(a;-2)' (x-aXx-b)' (a;-l)(a;-2Xx-3)' 
(a;-a)(a;-6)(a;-c)' {x''^-a%x^^\P') 

/ov l, X, X , X f , •, 



^^^ (rm;+ny (s"^^*!*"*® ma + « = y ), 

-7 — ^ . .^ (substitute aaj"*4-c = l/2/), 

1 1, Xy 01? 

/jj»2 __ ^2 J^ 

188. Generally, if aj— p denotes a real linear factor of 

the denominator Fa?, so that 

Yx = {x—p)({>Xy 

. {x A R 
we put =-= , 

¥x x—p (px 

and now ^^ix-R{x-p) ^ip^ 

<j>x <pp 

on putting x — p = 0. 

But {x — p)<f>x = Fxy 

^x+(x—p)(f/x = ¥'x; 
so that <pp = F'p ; and A = ipfF'p. 

Now the integral of the corresponding partial fraction 
A/(x - p) is A log(a; —p). 
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If (flj— a)*4-i8^ denotes a quadratic factor of Faj, splitting 
up into the conjugate imaginary linear factors 

flj— a— i)8 and x—a-^-i^y 
we take two corresponding partial fractions 

x^a — ijS x—a+i0 
coalescing into the real partial fraction 

the integral of which is 

ilog{(aj-a)2+^2|_2Jkftan-i{(a;-a)/i8}. 
We generally begin by assuming 

ix Hx+K 8 



Fx (x-a)^+l3^ ' V^ 
where Fx = { (a? — a)^ + fi^yyfrx ; 

and Hx+KJ^Z:^i^f^l±^=^, 

on putting (x — a)^ + )8^ = 0. 

This gives imaginary values of x ; but we may avoid 
the use of imaginaries by continually writing 2aic— a^— )8^ 
for ic^, when finally {x and {Hx+K)\frx are each reduced 
to the form Bx + G smd Bx+CT; and thus B = R, C=C, 
whence H and K are determined. 

189. When Yx has a factor (aJ— g)*", a linear factor 
x—q repeated r times, so that Yx = {x'-qyxXy we must 
assume r corresponding partial functions of the form 
Br Br-i . ^1 ^ f a? r . 

and then to determine the B's put x — q = yy or x = q+y; 
so that multiplying by y^, 

Br+Br.iy + ...+B,yr-'' = {{{q+y)^Tyr}/x(q+y); 
Br, Br-u ••• ^1 ai^e thus the coefficients of y^, 3/\ ... j/^-^ 
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in the algebraical expansion of f(9+2/)/x(?+2^) i^ ascend- 
ing powers of y ; or, in other words, 

Br+Br^xy + .^.+B^y^-^ 
is the quotient and Ty'^ is the remainder, when i{q+y) is 
divided algebraically by x(?4-2/) in ascending powers of y. 
Now the integral of a partial fraction J5«(fl? — g)"' will 

be , — ,";^' , , ; and of -A. will be BAogix - a). 

When Fa; has the repeated factor {{x — of + ^Y, 
we can proceed in the same way, and assume r corre- 
sponding partial fractions of the form 

{H^+K,)l{{x--af+^Y) 
but in this case it is generally preferable to employ the 
conjugate imaginary factors ; or else to employ the sub- 
stitution a; — a = /8 tan 6 or jS sinh u. 

^ - , where m and n are positive 

integers, required in the vertical motion of a body when 
the resistance varies as the n^ power of the velocity. 
K a? — a is a factor of aj~— 1, then 

a = l, or cos(2r7r/n)+isin(2r7r/n) = exp(2ir7r/7i), 
where r has the values ±1, ± 2, . . . 

For then, by De Moivre's Theorem (§ 111), 
a** = cos 2r7r 4- i sin 2rTr = 1. 

Put 4-1 = ^+2—'- ; 
x^—l x—l x—a 

here fa;=a;^*-\ Fa; = aj'*— 1, F'a; = naj~-^; 

fct cb^"^ cC"^ Oj^ 
so that jly=:=-= — r-7= — ^= — (since a" = l) 

1/ 2mr7r . . . 2mr7r\ 1 ^vmrir 
= -I cos h'2' sm ) = -exp , 

by De Moivre's Theorem ; also A = l/n. 
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Therefore / — - — r- = - log(aj — 1) 

, I ^/ 2mr7r , . . 2mrir\ ( 2r7r . . 2r7r\ 
+ — 21 cos h^sm )log( 03 — cos ^sm )• 

To express this result in a real form, the partial frac- 
tions with numerator Ar and A.r, corresponding to 
conjugate imaginary values of a, must be combined into 
a real form as follows — 

x-cos(2r'7r/n)-ism{2rTr/n) x-cos{2r'7r/n)+i sin(2r7r/?i) 
1 2mr7r 2a; — 2 cos(2r7rM) 

= cos ^ ' ^ 

n n x^^2x cos(2r7r/ti) + 1 

_ 2 . 2mr7r sin(2r7r/n) . 

n n x^-'2x cos (2r7r/^) + 1 ' 
and the corresponding integrals are 

—cos logfaj^ — 203 cos h 1) 

n n ^ n 

— sm tan"M{a3— cos I /sin }-. 

n n W n // n J 

When the denominator is x^+l, the result is of the 
same form, with 2r— 1 written for 2r; and from the 
above we infer that the typical quadratic factors of 
03^ — 1 and x'^+1 are 

03^ — 203 cos(2r7r/7i) + l, and 03^ — 2a3COs{(2r— l)7r/7i} + l. 

If m is changed into n — m, as for instance by the 
substitution 03 = 1/2/, then C08{2mr7r/n) is unchanged, but 
sin(2mr7r/n) changes in sign, so that.(Euler) 



/ 



'^m-l + ^n-m-l 

-ao3 



03'*±1 

is expressed entirely, either by logarithms, or inverse 
circular functions. 

Degenerate cases occur for n = m, or 2m, 
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When n is even, =2jp, and m is odd, =2g + l, then, 
by a rearrangement of terms, 

229*^ 2j3 ic^— 1 

Taken between the limits and ex, the part depending 
on the hyperbolic functions vanishes ; and the result is 

We may put a;^^ = 'v^, and (2q+l)/2p = m/n, where m 
and n are any integers ; and now, with the restriction 

thatm<'W', / -, . 1 =- 



Thus, for instance, the area of the loop of the curve 

on putting y==xv (% 63), is given by 

h I x%lv = ia^ I ^i— — -~ = -pr- cosec . 



Putting v^-=t, then with the restriction that m<l. 



cosec — . 

71 



TT cosec mTT. 







To determine y (tan 0)^/^dft put tan0 = ajP; then the 

. 

integral becomes 

' 

(tanh uY^du, 
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Examples, — ^Resolve into partial fractions and integrate 
l,x x^a?,„ j^^ ^ = 4, 6, 8. 



(2) 



(the integrations required when we integrate 
(sec Oy, (sec Oy, ... or (cosh uy, (cosh uy, ... or 
cos^ft cos^0; ... or (sech uy, (sech uy, ... by means 
of the substitution tan ^6 or tanh ^u = x, 

191 . IntegrationofanlrrationalAlgebraicalFunction, 
In integrating the irrational part NjD^R we may- 
suppose the rational function NjD resolved into its 
quotient and partial fractions ; and now we shall con- 
sider the integration of the simplest elements, 

1 1 . Hx+K 

JB; {x-p)JR' {Ax''+2Bx+G)JR' 

corresponding to the constant term in the quotient of 
N/D, and to the partial fractions corresponding with the 
linear and quadratic factors x—p and Ax^+2Bx + C of 
the denominator D ; and afterwards investigate formulas 
of redvyction for 

^ 1 . Hx+K 

JR' {x^qYJR' {Ax' + ^Bx+GyjR' 

By formulas (i), (w), (x) (§ 39), 

,, . f dx , .X , ,/-, x^ \ 
(i-) /— 77 — 9 27= sm-^— = cos-^;^il 5); 



("•) / // 9 . — o7 = sinh-^— = cosh-^^(H — „) ' 


(iii-) / // o 9\ = cosh - ^— = sinh ■ ^^( — « — 1 ) ; 



TO 



INTEGRATION IN GENERAL. 397 

in which it will be noticed that the analogy with the 
first integral breaks down when only the logarithm 

is employed to express the second and third integral. 

The quadratic R = ax^ + 2bx + c can always be expressed 
as the sum or difference of two squares in one of the 
three forms rn? — x^, m^ + x^, or x^ — m^; and then the integ- 
ral ycfo/^^it has the corresponding form (i.), (ii.), or (iii.). 

But substituting R = y\ we find 

ax + b = j^(ay^ + b^ — ac); 



/*dx _ r 
JR J ^ 



dy 



so that . __^-^^^__^____, 

and now, as in ex. 13, p. 77, 

(i.) when a is negative, but W — ac positive, 
dx 1 .1 ^(-a)y 



y*dx 
JR 



/: 



/-. 



(ii.) when a is positive, and b^—ac positive, 

(iii.) when a is positive, and 6^ — oc negative, 
^ =1 cosh-i x/Wx/£ ^1 3inh-i-gg+L. 



We cannot have both a and h^ — ac negative, because 
R would then be negative, and JR imaginary, for all 
real values of cc. 

When a = the integral assumes an indeterminate 
form (§ 120), and reduces to an algebraical function 

JRjb, or J{2bx+c)lb, 
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When b^—ac is positive, R breaks up into real linear 
factors, say x—a and x—^; and now 

(1) when a is negative, we may without loss of gener- 
ality put a= — 1 ; and 

/' dx _ r dx . i ,J{a—x.x — fi) 

= 2 sin-\/^= 2 cos-\/^=2 tan-\/^:^. 
\a—p ya — p y a — x 

(2) when a is positive, we may put a = 1, and 

/dx _ f dx _ . , _i </(a!— a.a; — j8) 

JIt~JJ{x-a.x-^)~ ^'°'' Ka-/8) 

= 2sinh-7^=2co.h-y^^=2 tanh-7|^ 

We may write, by analogy, 

y R~J{ac-l^) ^^^ J{ac-b^)' 

thus exhibiting the result as a function of R. 

For instance 

/** do? __ 1 _^^{x^-2axQ,o^a+a^) 

Vic^ — 2aa3COSa + a^ a sin a a sin a ' 

r^ dx _ 1 , _ J (x^-2ax cosh^+a^) 

y x^-2axcosh^+a^~~ asinh^ asinh^ ' 

y^"^ dx _ 1 ^_^^{x'^-2ax^\vAiy-a?) 

-9 ^^ ; ir 9 — V cosecn 5- , 

x^-2ax sinh y-a^ acoshy a cosh y 
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-——-J-, we may substi- 
tute x — p = l/y, when we obtain an integral of the same 
form as fdjx\JR ; but it is more direct to substitute 
y = ^RI{x''p)\ and now 

2 R+\{x-'pf _ {a+\)x^+'2.Q)-p\ )x+c+\p' . 

b ^~ CLP 

a perfect square when >^= 2^2bp+c ' ^^ ^^^^ 

and taking the negative sign with the radical 

(ap+b)x+bp+c ^ _ V{(apH26p+%H6''-ac} 

dx _ ydy 

{x-pf~ »y{(ap^ + 2bp + c)y^ +b^+ac} 

J {x-f)JR ~J(x-ffy ~J J{{af+Wp+cyy^+W-ae} 

1 . _i s/(.—ap^—2bp—c)y 

J{-af-2bp-cf ■ J{b^-ac) 

1 „;,,-i x/(-V-%-c),>/.R .. 

~ Ji-ap'-ibp-c) J{b^-acXx-p) ■•••^•'' 

1 . u - i JJf^f' + 2bp + c)jR 
^"^ ~ y/{ap^+2bp+c) ^^''° V(^' - «c)(a; -p) ^"•■' ' 

^{ap^+2bp+c) '^'^'^ J(ac-¥){x-p) ••■•^"'•''' 
the real form to be chosen (ex. 18, p. 78). 

Here again, when a'{?+2bp-\-c = 0,so that x—p is a 
factor of R, the integral assumes an indeterminate form 
and the limiting value is the algebraical function 

V(fe*-ac)(a;-p)" 
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193. We may * consider the integxa,iydx/(x—p)^R as 
the degenerate form of the elliptic integral ydx/^X, 
where X is a quartic function of x, as in ex. 39, p. 82, 
when X splits up into the factors (x—pY and R. 

Similarly the canonicaZ elliptic integral J'dsj^Sy 
where, as in ex. 39, p. 82, 

8=-^^-g^8-g^= 4(s - e^s - e^s - e^\ 

when resolved into factors, becomes one of the degenerate 
forms of exs. 14, 15, p. 77, if two of the quantities 
e^i eg, ^3 are equal. 

We suppose e^>e2>e^] and that the middle quantity 
^2 is made equal either to the smallest Cg, or the greatest 
e^ ; and then with (i.) ^2 = 63, as in ex. 14, p. 77, 

y^da _ r ds 

JSJ2{x^e^)J{x-e^) 

= -y. ,sin-\/ ■=0-77 rsm-^ ^^ ^ ^ ^; 

,.., fds 

(u) e,=e„J-^ 

Similarly, as in ex. 15, p. 77, 

r da _ 1 ^^., l e^-x , 

J2{e^-x)y/{e^—x) Jie^-e^ Se^—x' 



=x7(^'^''^'VJ4('^>*»>- 



or 

. Up. .^p. \ \ yv. A^ 
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In physical problems we require the integrals 



^r_J^ and C=f- 



00 



A ' - — ' ^^ 



M /"" d\ _ 2 .1 /c«+X 



or = — TT-o iiT cosh 



-V; 



cHX 



and, differentiating with respect to a* or <?, we deduce 

il = ^*^-^ ^ h , cos '*/ , 

(a^ - c^Xa^ + X) («.2 - O^ > a^ + X 



or 



(c2 - a)(a2 + X2) (c2 - a2)t V a^ + X 

n 2 2 .1 /c2+X 

0= ; COS ^x » 

2 2 , , /c2 + X 

— -. cosh \/ 

(c2 - d')J{c^ + X) (c2 - a2)^ \ a2 + X 

2 



so that 2A'{-C= ________ — _. 

(o^HX)V(cHX) 

194. To determine the integral 

Hx + K dx 



A 



Ax^+ 2Bx + C J{ax^ + 2bx + c)' 
where we may suppose A and therefore Ax^-h^Bx + C 
positive for all real values of x, substitute 

y = (ax'- + 2bx + c)/{Ax^ + 2Bx + C). 
Now, with the notation of ex. 8, p. 146, 

dy __ 2{A h — aB){x-^ — x){ x — x^) 
dx~ {Ax^+2Bx+Cf~' ' 

^^ ^ Ax^+2Bx+C ' ^ 2/2- Ax^+2Bx+C ' 

where y^y y^ denote the maximum and minimum of y^ 

and x^y x^ the corresponding values of x. 

2c 



/ 
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We may write L(x—x^ + M{x^ — x) for Hx + K, and 
V=Ax^+2Bx+G=P(x^-'xf+Q(x-'X^)^ 
R = ax^ +ibx +c ='p(x^ — x)'^+q{x — x^^y 
as before, on p. 147 ; and now we find 

' L{x - x^ + M{x^ - x) ^^ 
VJR "^"^ 
_ r L(x - x^ + M{x^ - x) VHy 

J V^Jy 2{Ah - aB){x^ - x){x - x^) 

so that the expression for the integral is 

- LL' cos - VC2//2/1) + ^I^ cosh ~ V(2//y2)> 
or ii' sin - V(2//2/i) + ^^^ si^h - V( - y/y2l 

according as 2/2 is positive or negative ; that is, according 
asac — 6^ is positive or negative; that is, as R is always 
positive, or as i? can vanish, for real values of x. 

When V = R, these forms become illusory ; and now 

fR-Hx = {ax + 6)/(ac - }?-)JR, 
/(ax+b)RMx = R-^; 

whence y(Hx + K)R~^dx can be determined, as an alge- 
braical function of x; and similarly y(^ia::+ir)JB"^cia;, 

Again, when Ab—aB = 0, the results are again illusory; 
but in this case we can choose a new variable, by changing 
x + b/a or x+B/A into x, so that we may make 6 and B 
vanish ; and now the integral 

*Lx+M dx 



/: 



Ax^+C ^{ax'^ + c) 
consists of two paris, of which the first is a function of 
x\ of the form given in § 193 ; while, as in ex. 19, p. 79, 



or 
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/• Mdx _ M J IfC ax^+c \ 

M ■ x. X If C aa?+c \ 

JCJ{aC-Ac) ^'""^ iK—cA^D)' 

reducing, when a(7-^c = 0, to (T-77~2X"3V 
By differentiation of the integral 

y ^ Hx+K dx 
Ax^ + 2Bx+C J(ax^+2bx+c) 

with respect to A, B, or (7, we can deduce the results of 

ix dx 



Jh 



{Ax^ + 2Bx + C)^ s/{ax^+2bx+c)' 
The general linear substitution may be written 

,_ e(x^-x)+f{X''X^) 

E{x^-x)+F{x-^x^' 

and the form of the integral will then be unchanged; 
in particular V and R are interchanged if we put 
E=J(Ppl F=J(Qql e = L\f=Fx,, 

Exam'ples. 

(1) Integrate 

x — \ x+ 1 

{x^+x+\)J{x^'-x+l)' {x'-x + \)J(x^+x + \)' 

L(x-l)+M(x-2) 
(3^2 - 1 Ojb + 9)^(5a;2 - 1 6a; + 1 4/ 

L(x-l) + M(X''2) 
(5x^'-lQx + U)^(Sx^-'10x+9) 

Hx+K 

(3x2 - lOic + 9) V(a;2 - 8aj + lOy 
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(2) Prove that the maximum and minimum of 

{ax^ + "Ihxy + h'f)\{A7?' + "IHxy + By^) 
are given by the common conjugate diameters of 
ax^ + ^hxy + hy'^ = c, Ax^ + 2Hxy + By^ = C. 

(3) With the notation of ex. 8, p. 146. prove that 

((a? + 6 = — ^{x-^ —x)+ q(x — x^\ 
hx+c= 2^i(^i — oc)— Q^2(^ "" ^2) j 

Bx+C= Px^{x^'-x)'-Qx^{x~'X^\ 
and (52-^a)i?2+(^c+aa-256)iiF+(62_ac)F*^ 

= (^6 - aB)\x^ - a;)2(a; -ccg)^. 
If a;' = - (fiic + C)/(^a; + JS), or - {hx + c)/((/ir + &), 
jR R_ Ac+aG-2Bb F. V _ Ac-{-aC-2B b 
V'^V'~ AC-B^ '^^]K'*"iJ'~ ac-62 "' 

(4) Prove that if X denotes the reciprocal quartic 

ax^ + 46x^ + 6cx^ + 46aj + a, 
/'^''-l do; 1 J(-a)^X 

^U ^r JX-Ji- a) ^'° s/{-n)x' 

where aX = (ax^ +2bx+af+ Dx\ D = 6<tc - 46^ - 2a\ 



("•l/l 



ir^-l dx 1 . .i^(2a-6c)V-X' 



or 






and (Gc - 2a)Z = 4 { hx^-{a - 3c)x + 6 } « + D(x« +l)\ 
(5) Determine in a similar form 

/x^ + \ dx , fx^ + l dx 
' IT 'JX' ^'^ V ^f^\ 'JX" 

where A"'=ax*— 46a;'+6cx^+46a;+a. 
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195. Integration by Rationalization. 
Any integral of the form 

r S+T(ax±byi^ . 
J U -^V^ax-^hf"^^' 
where S, 1\ U, V are rational integral functions of x, 
can be rationalized, that is, can be made to depend on 
the integration of a rational algebraical fraction by the 
substitution ax + b = y^'j for then aclx = ry^''^dy ; and the 
integral becomes the rational algebraical integral 

rS + TyPryr-^ 

J u+v^~ir^y\ 

Any integral of the form 

f S + TJR^ ^ 

J u+vjr"^' 

where jB is the quadratic form ax^ + 2bx + c, can be 
rationalized, 

(i.) when b^ — ac is positive, and a negative, by the 
substitution 

V(-a)V.R _ 2y ax+b _ l-y^ 

J{b^-ac) ~l+i/ ^^" Jib'^-acyi+y^ 
, —adx _ ^ydy 

^"^ Jib^-acrJX+y^r 

(ii.) when b^ — ac is positive, and a positive, by 
^{a)^R _ iy ax+b 1+y^ 
J(,b^-ac) 1-y^' J{b^-ac) l-y^' 
. adx _ 4iydy 

^""^ J{b^-ac)-ly^^' 

(iii.) when b^—ac is negative, and a positive, by 
J{a)JR 1+y^ ax+b 2y 

^{ac -b'') 1 - y^' y/(ac - b^) 1-y^' 
adx . 1 ~l~ ^ 
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In (l) we may suppose 2/ = tan J0, and in (ii.) and (iii.) 
2/ = tanh Ji6; and then the integration is changed to an 
integration of a rational function of cos Q and sin Q, or 
cosh u and sinh u^ with respect to Q or u. 

Examples. — Integrate with respect to x 

x^ x^ 1 

s/(x-a)' ^(x+a)' *V(«^-«). (r+x)Jx 

1 1 X 

x*^(l+x^)' {x+l)^(x+2y {x+2)^(x+l)' 

1 1 1 

{x+l)J{2x+l)\l-x)^il-x^)\l+x^)y/(l-a^)' 

V(a"+x^) 1 . 1 

a+x ' Xy/(x^+-Sx+2)' {x + l)^{x^+x+l)' 
1 1 bx+c 

(x^+1)^' (ax^+2bx+c)^' (ax^+2bx+c)i' 
Hx+K 1 ^, 

(ax^+2bx+cf a^^{ax^+c)' ^^~^* ^'' 

Jix'^+aKx^+b^)' {x*-l)i' ^(^+^)' 
{ax^ + b)'~ir, (substitute a+bx-^ = y). 

196. Integration of Circular and Hyperbolic Func- 
tions. 

To integrate powers and products of cos a? and sin a? 
the most general plan is to convert them into cosines and 
sines of multiples of x^ which are immediately integrable 
(ex. 9, § 40). 

To integrate any odd power of cos a; or sin a?, say 
(cos x)^^'^^ or (sin 0?)^**+^, we write them in the form 

(1 — sin2ic)'*cos X, or (1 — cos2a;)'*sin x, 
and expand by the Binomial Theorem ; then each term is 
immediately integrable, since by § 40 and (a), p. 84, 
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(sin a;)"*cos a^ = 



> 



m+1 
(cos a;)"*sin fl%?a; 



_ (cos g;)^+^ 



m+1 

A similar method will serve to integrate (sina;)^(cosic)^, 
where either p or 5 is an odd integer ; also to integrate 
any powers of vers x. 

The same processes apply when cos a? and sin a? are 
replaced by the hyperbolic functions cosh x and sinh x. 

Examples. — Integrate with respect to a?, coswa?, 
sin(ma; + 7i), sin 2a; cos 3fl?, cos 3aj cos 5a;, sin 3a; sin 5a;, 
sin (ma; + 7i)cos(pc + q\ sin x sin 2x sin 3a;, sin x cos a;, 
sin2a;cosa;, (sina;)^cosa;, (sin a;)"*cos a;, sina;cos^, sinajcos^a;, 
sin-a;(cosa;)*^, sfn^, sin^a;, sin^a;, cos^a;, cos^a;, cos*a;, 
cos^a; cos 'yia;, cos^ma; cos -na;, .... also the same functions 
with cosh X for cos a;, and sinh x for sin x. 

1 97. The integration of the remaining trigonometrical 
functions is a little more complicated ; thus, by (v), p. 85, 

y/^COS X 
c>oixdx = / - — dx= log sin a; = J log vers 2a; ; 
^x sin X 

y #^sin X 
iaxLxdx = / dx = — log cos x = log sec x ; 
^/ cos X 

while ycoth xdx = log sinh x, yianh xdx = log cosh x. 

Again Jiditi^xdx=J{ sec^a; — l)da:= tana;— a;, 

fQ.o\?-xdx =y(cosec^a; — V)dx = — cot x-^x; 
with similar results for tanh^a; and coth^a; ; while other 
powers of tan x are integrated by expressing them as the 
sum of terms of the form (sec a;)^tan a; or (tan a;)^sec^a;, 
the integrals of which are (secxy/p and {ta,nx)^'^^/(q+l); 
similarly for powers of cot a;, tanh a;, or coth x. 
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To integrate sec x and cosec a?, we may rationalize by- 
means of the substitution tan ^x = y ; and then 

1-y^ . 2y 

cos aj = r--^, sin a; = =— ^, 

1+7/2 l^y2 

while x = 2 tan - ^y, dx = 2dy/(l + y^). 

Then /sec ircZa* = / ., % .. . ^9 =/^ — ^o 

y y 1-2/^ 1+2/ y 1-2/ 

= log(l+2/)-log(l-2/)=log(l+2/)/(l-2/) 
, l+tanAaj , , ,, , , ^ , /l+sin« 

= log(sec X + tan x). 
Similarly /cosec ajcZx = / ^ ^ TX2 

/^ /I "^ cos a? 

=/d3r/2/ = log 2/ = log tan h^x = log^ j^:^^ 

= log(cosec X — cot x). 
Similarly to integrate sech x or cosech x, we may 
rationalize by the substitution tanh ^x = y ', and then 

/ sech xdx = f ^ ^ = 2 tan ~^y = 2 tan ~ ^tanh ^x, 

/cosech xdx = /dy/y = log y = log tanh ^x. 

The substitution of y for tan|0 will apply for any 
function of the form f (cos 0, sin 0)/F(cos 6, sin 6) ; the 
function being thereby reduced to the form 

n(2/-6) nsinKg-^) . 

U(y-a)^ Usin^ie-a)' 
and this again by partial fractions (§ 187) to the sum of 
terms of the form A/{y — a), or 

A cos|0 cos Ja cosecj(0 — a) = ^ cos^^acot J(0 — a) — 1-4 sina. 
(Hermite, Froc. London Math. Society^ vol. IV.). 
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1 98. By the substitution of z for sec x, cosec Xy sech Xy 

cosech X, the results can be expressed more directly ; thus 

dz 
with z = sec x, dx = — jy-^ — r-, 

and /sec xdx = /— 77-2 — f\ = cosh ~ ^0, by (x.) (p. 85) 

= cosh " ^sec X = sinh " Han x = tanh " ^sin x = 2 tanh ~ Han |a;. 
Similarly, with = cosec x, and as a corrected integral, 

yri^ f dz 

cosec aj<^ = / , g — r-r = cosh ~ '^z 

= cosh " ^cosec x = sinh " ^cot x = tanh " ^cos x. 
Again, with z = sech aj, or cosech x, 

/sech cccZx = /— ^— — 2^ = cos ~ ^0 

= cos~%ech X = sin"^tanh x = tan"^sinh a? = 2 tan"^tanh Ja;. 

y cosech irdaj = / ,., . — ot = sinh " ^0 

= sinh " ^cosech x = cosh " ^coth x = log coth ^x. 

199. To integrate 1/P, where P = a+6 cosx+csinaj, 
we can proceed in a similar manner, and put 1/P = y; 

and now — a = b cos ir + c sin x, 

y 

so that — 2 -^ = /;sinir — ccosa;= a/J62 + c- — f — a\ !■ ; 

■j /* dx f_ dy^ 

y a + 6 cos a; + c sin aj"l/^{( — a^+ 6^+0^)2/^+ 2a;j/— 1}' 

which by the results of § 191 can be expressed by 
(i ) __ J cos - ,«+( -aH6Hc% 
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^^("•> ^(-aHfeHc^) x/(62+c2)P -' 

as in ex. 34, p. 81. 

Similarly to integrate 1/Q, where 

Q = a + 6cosh cc+csinhic, 
we put 1/Q = y ; and obtain the results of ex. 35, p. 81. 

Between the limits and tt, (i.) ='7rl^{a? — b^—(^) 
but (ii.) is illusory, as the function to be integrated 
becomes infinite between the limits of integration. 

Again, between the limits and oo , 

fdxjQ oy:fdxj{a + h cosh x+c sinh x) 
is given by means of ex. 35, p. 81, as 

^(.a2+62-c2)^^ a+b+c 

2 tanh-ir V(a^-6Hc^) l 

^^ V(^~'-62+c2) coth-H a+6+(; f 

200. The integration of (a +6 cos 0)-** is effected 
(i.) when hja < 1, by the substitution 

a+6cos0= T , or idi,n\6=^Jv u^^h'^y 

equivalent geometrically to a change from the true 
anomaly 6 to the excentric anomaly v (§ 177) in an 
ellipse of excentricity b/a ; and now 

= 7 /(a — 6cos'y)**~^rfv: 

(a + 6cos0f (a^--b^y-iJ 

(ii.) when b/a > 1, by the substitution 

a change from the true anomaly 6 to the hyperbolic 
excentric anomaly u in a hyperbolic orbit of excentricity 

b/a ; and now 



/ 

J a 
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= /(bcoshu — a)"-^du. 

When n = l, we obtain, as in ex. 31, p. 80, 

dO _ V _ 1 _^ a cos 6+b 

a +b^^^e ~ y/{a^ - b^) ~ ^(a^ - h^) ^^^ ' a+b cos ff 

_ u _ 1 ,_jaeos0+6 

^^ - 7(62^:^) - ^{b^^a^) ^^^^ ' a+b cos ff 

The integral for 71 = 2 is required, as in § 180, for the 
mean anomaly nt in an elliptic or hyperbolic orbit. 
Reciprocally, since cos 6 = (a cosh u — b)/(a — b cosh u), 

f- — = ^- i f{a+b cos ffT'He, 

J (bcoah u-ay (b^-^a^'^J 

including all possible cases required in the integration of 

(a+)8coshu)-^ 

The results for 7i = 1 will be found in ex. 32, p. 80. 
To integrate (a+tsinhit)"**, substitute 

a+bsinhu=^(a^ + b^)/(a — bsiiih v) ; 
and now 

/ r-rh — r = —^ — v^ — i /(^^ - ^ sinh vf ' ^dv. 

J {a+b sinh u)~ {a^ + b^ '^J^ ^ 

If we put 6 = gdu, <p=:gdv, and &/ci = tany, 

then smhu = ia,n6 = t-^ — =tan(0 + y) : 

a — otan^ ^ ^^ 

so that 0=<p + y ; and now 

and the result for n = l is given in ex. 33, p. 81. 
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With the notation of the eonfocal conies of § 177, 

tan \Q = coth \u tan \v, tan |0' = tanh ^n tan ^v ; 

and taking the logarithmic differentials 

dQ --du dv dQ' du dv 
sin 6 sinh u sin v sin 0^ sinh iv sin r ' 
while 
(cosh u — cos v)(cosh u + cos 0) 

= (cosh u + cos 'y)(cosh u — cos 0') = sinh^u, 

(cosh u — cos t;)(cos v — cos 0) 

= (cosh u + cos 'y)(cos 6' — cos v) = sin^t; ; 
so that 

y'*(sin'y)^'*~^(ii6 /7 ^, s^j/\f i. 

.^ , ^ tt: =/(C0S 0' - COS vYd& \ ctC. 
(cosh it 4- COS v)^ y ^ 

As applications the student may evaluate the expres- 
sion for the area PpP'j/ bounded by the elliptic arcs 
It, u' and the hyperbolic arcs v, v' of the confocals of § 177, 
or bounded by the circular aras n, n and v, v' of the 
system of dipolar circles of g 175, given by the integral 

II - 1 ' ^' dudv = hcyAcofih 2u - cos 2v)dudi\ 

c^dudv 
(cosh V + cos u)^ ' 

Prove also that the c.G. is at N (fig. 54, i.) of 
(i.) a circular wire in which the line density varies 
as SP"^; (ii.) a circular area in which the surface density 
varies as SP'^ ; (iii.) a spherical shell in which the surface 
density varies as SP'^ ; 



or H — ^^_^^^^^ — _ . ^^^ determine the centroid. 
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(iv.) a solid sphere in which the volume density varies 
as SP'^; the boundary in each case being defined by 
u = a positive constant, so as to include N in the boundary. 

Determine also the CG. of a quadrant of an elliptic 
wire, as in fig. 55, in which the line density varies as p, 
1/p, xp, x/p, . . ., 2^ denoting the length of the perpendicular 
from the centre on the tangent. 

201. The integrals 
(i.) /(H cos e + K)dei{A cos'^0 + 2B cos + C\ 
(ii.) yiH cosh u + K)du/(A cosh^it + 2B cosh i6 + (7), 

(iii.) y{H sinh u + K)dv/(A sinh^t + 2B sinh u + G), 
are reduced to the form of the integral of § 194, by the 
substitution x = cos 6, cosh u or sinh v ; and now jR is 
replaced by sin 6, sinh u, or cosh v ; while by writing 

sin^g = {(1 — cos a cos 6)^ — (cos 6 — cos a)^} /sinV 
F=J.cos20 + 25cos0 + (7 
= P(l — cos a cos Of + Q(cos 6 — cos a)^, 
and y = sin^a sin^Q/ F; 

l-P2/ = (P+Q)/(cos0-cosa)VF, 
l + Q2/ = (P+Q)/(l-cosacos0)VF; 
then 1/P and- l/Q are the maximum and minimum values 
of y; and the integration is expressed by inverse sines, 
circular and hyperbolic, of s/i^V) ^^^ s/iQv)' 
Similarly for cosh u and sinh u, by writing 
sinh^u = { (cosh /3 cosh u — 1 )^ — (cosh u — cosh 13 f} /sinh^^, 
cosh^i^ = {(sinh y sinh u + 1)^ + (sinh u — sinh y)^}/cosh^y. 
As numerical examples, integrate 
{X(3 - cos 0) + Jl/(1 - 3 cos 0) }/(5 cos^g - 6 cos + 5) 
{i(cosh u — 2)4-i/(2 cosh u— 1)}/(5 cosh^u — 8 cosh it +5), 
( H sinh V + K)/(sinh^v — sinh v+1), 
(H sinh i; + K)/(Q sinh^v - 4 sinh t; + 9). 
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Eaximples, — Integrate with respect to «, 
(1) tanajsec^a;, tan^a; sec^o?, tan^a; sec^aj, (tan a;)"»sec^aj, 
cot a; cosec^x, cot^a; cosec^a?, (cot a;)'^cosec^a;, tanx, 
cot a?, tan^a?, cot^a;, tan^a;, cot^aj, tan^a;, cot^a?, 
sec a? tan aj, sec*a?tana;, (secajy^^tana?, (coseca;)"*cota;, 
sec^a;, sec^aj, sec^a;, cosec^a;, sec a;, sec'a;, cosec^a?, 
sec X cosec x, vers aj, vers^aj, 1/vers x ; also the cor- 
responding hyperbolic functions of x, 

1 1 1 1 

^ ^ cosaj+sina;' bcosaj+csina;* a + 6tana;' a + 6cosaj* 

_1 1 1 1 

5 4- 4 cos aj' 4 + 5 cos x a^cos^x ± )8%in V a^cosh^a; ± jff^sinh'^x' 

jB cos a; 4- (7 sin a; J.+ ^cosa;+(7sina ; 
6 cos aj+c sin aj' a+6cosa:+csina; ' 

~„ i . .y . o ,seca;sec2a;,tana:tan2a;. 

(3) Prove that 

(i.) sin X sec(a; - a)sec(a; — b) 

= cos a cos 6 cosec(ci-6){sec a sec(aj-a)-sec 6 sec(a:-6) } ; 

,.. . sinaj _^ sin a cot(aj — a) 

^ '^ sin(a? — a)sin(a;--6)sin(a;— c)~ sin(a — 6)sin(a — c)' 

..... sin^a; _^y sin2acosec(a; — a) 

sin(a; — a)sin(a; — 6)sin(aj — c)~" sin(a — 6)sin(a — c)' 

,. . sin ma; 1 v>/ in/. • i,, ^ (wherea = A;7r/7i, 

(iv.) — =;^2(-lrsinmacotA(a;-a) , J ' 

^ sin7ia; zn ^ ^^ / sma rii<n) 

= ^j-2(-l)^*sin ma cosec(a;-a), or ^2(- l)*sin ?7iacot(a?~a), 

according as "iu + n is odd or even. 
Write down the integrals of these functions. 

(4) Integrate sec(a;-f a)^{sec(aj + 6)sec(a; + c)}, 

by means of ex. 16, p. 78, putting tan a? = 2/. 
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(5) Prove that the definite integrals 

tan xdx =/cot oiAx = log,^2 = 0-34657. 



(ii. 



(in. 



(iv. 



(V. 



(vi. 



(vii. 



(Vlll. 



(6) (i. 



J sec xdx =J cosec xdx = log(^2 + 1) = 0'88137. 

iir 

/ '"^ dx _ IT 'f >/ 

a + b COS x^ >s/(fl^^ — 6^y 



/i^ dx _ Ji _!& 
a + b cos X "" ^{a^ — 6^) a' 

P^T^a =0-23180. 
^ 0+3 COS 





/: 



h"^ dB 

^^ =0-27465. 



3 + 5 cos 





/ 



h^ d6 

""^ . = 0-03494. 



(5 + 3 cos 0) 



2 





/"IT H cos 6+K ,^ _ jff sin|a + Kcos^a 

jcos^Q - 2 sin a cos + 1 "" ^(2 cos^a) 





if cos a is taken positive. 



... r^ HcosO+K 

^ ' y a^QOS^Q — 2ab cos a cos + b^ 



de 



where c* = a* — 2a26^cos 2a + 6*. 



uv 
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202. Integration by Parts. 

Corresponding to the formula for the differentiation of 
a product, by parts (§ 12), 

duv ^du dv 
dx ~^dx dx 
is the formula, obtained by integrating both sides, 

/dii ^ , r ^'^ 1 
-7- vdx + 1 u-r-dx. 
dx J dx 

y' dv , r du 7 

u-T-dx^uv—i -r-vdx, 
dx J dx. 

or, as it may be written, 

Judv = uv "Jvdu, 

the formula of the Integral Calculus, for integration by 
parts, as it is called. 

Interpreted geometrically the formula merely asserts 
that, as in fig. 19, with u and v for x and y, the rectangle 
OMPN is the sum of the areas OMP, ONP into which 
the rectangle is divided by any curved line ; or that the 
area ONP, represented hyyiidv, is the difference between 

, the rectangle uv, and the area OMP, represented hyjvdu. 

The formula shows how the integi^al Judv may be 

made to depend upon another integral Jvdu, which is 

either integrable, or which is made to depend upon 

another integral, by another integration by parts. 

Thus (i.) to integrate xt^, take x = u,-^ = e^\ then v = €*, 

and fxe^dx = xe^ —fe^dx = x^ — ^. 

(ii.) to integrate \ogx, take tt = loga?, dv = dx, v = x; 
then yiog xdx = x log x '-fxdx\x = x log x — x, 
(iii.) sin ^x, put ^x = z, then 

ysin ^xdx =ysin z . ^zdz 
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= — 22; cos 0+ 2/cos zdz = — 2^ cos z+^ sin z 
= — 2^x cos ^x + 2 sin ^x. 

\1 (iv.) /sin"^ircZa: = a:;sin"^x— l—iT^~,r. 
J J J^-^^) 

= oj sin ~ ^aj + ^(1 — x^). 

(v.) to integrate cos wjx cosh nx \ here cos rax or 

cosh iiaj may be taken indiscriminately for u, and we 

must integrate by parts twice. 

•r, . . .« , dv sinhtia; 

J^ or instance 11 cos ma? = i6, cosh nx = -=-, t; = , 

' dx n 

M 1 771/ t 

/cosma7COshnajcia; = -cos77iajsinhna;H — /sinma;sinh7iajcte 
^ n nj 

=-cos majsinhTiajH — ssinmxcosh^ia; x /cosmxcosh nojdaj; 

71 tv' nU 

and therefore, by transposition and division, 

, r 17 msin7/ixcosh7ia?+7icos77ia;sinh7iic ,. . 
: /cos majcosh nxdx = «- — 5 (i-) 

^ 77l* + 71^ ^ ^ 



In a similar manner 



/. 

/ 
> 



. , , 7nsin7naJsinh7ia;4-7icosmaJCOsh7ia; ... , 
cos77iajsinh7iajaa5= k-, — 3 (11.) 

771^ + 71^ ^ 

, , TisinmajsinhTioj— 77icos7najcosh7ia; .... . 
sin7?ia5cosh7ia;aa; = «-; — s (111.) 

. , , 7lsin7?^a;cosh7la; — 7ncos7n^sinh7ia; /. v 
sin77WJSinh7ia;ow3 = ^- — 5 (iv.) 

By addition of (i.) and (ii.), (iii.) and (iv.), or independ- 
ently by integration by parts, 

7 ^^TicosTn/OJ+Tn/sinTnaj , v 

e^Q,o%Tfvxdx'^ t^ 9- — 5 (v.) 






7 ^^Tismmaj — TMcosTMo; , .. 

e^sin77ia;cte = e"* 5-1 — 5 (vi.) 

2d 
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With a subsidiary angle a given by tsina = m/n, 
ie^cos mxdx = e***cos(ma; — a)eos a/n, 

/e'**sin viccdx = e^ sin(ma5 — a)eos a/n. 
Again, by successive integration by parts, 



d'^'H du d^-y 



dx dx^'^ 



the integrated terms of which can again be integrated by 
parts. 

We can establish Taylor s Theorem by successive inte- 
gration by parts ; for 

f(a+A)-fa=/faj^; 

a 

= - {(a+h-x){'x} + /{a+h-x)i"xdx 

a 
(l '\a-\-h I /-a+h 

= M'a-|2,(a+A-aj)2f"a;| +j^/ia+h^xyrxdx 

a 

=hra + %ra+ ... +'^7^"^ +^^{a+h^xr£^+^xdx; 



a 



SO that, as in § 114, 



f(a+A)-fa-.Afa-5V«- - "-^f"« 
= -B=^ /(a+h'-'x)H^+^xdx. 
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It is assumed here that fa;, i'x, i"x, ... f**a;, P*+^flj are all 
finite between the limits a and a+h oi x; and now if P 
denotes a certain average value of f'^+^oj between the 
limits a and a+h, which we may denote by i^'^\a+dh), 



i 



Examples. — Integrate by parts, 

\ a;"*Iogaj'*, (loga;)^ (logo;)^, (log a;)*, (loga;)^ ... ; 
^ oj^cos X, a^sin x, . . . ; cos^/^c, sin^a?, . . . ; 

e^+*cos(7naj + n), e***+*sin(wia; + n), . . . ; 
^cos'^a;, tan'^a;, cot"^a;, sec^a;, cosec^a;, vers'^a;; 
cosh"^, sinh"^a;, tanh"^a:, ... ; 
(sin'^a;)^ (sinh"^a;)^ (cosh-^aj)^ (vers'^a;)^; 

tan'^^^a;, tanh"^4/^» ••• > 

s/ia^-x^l ^(x^-a^l ^(a^-^a?\ (a« - a;^)*, 

ajsin^^a;, a^^cos'^a;, a;*cosh"^a; ... ; 

log a; . sin ~ ^a;, log x . cosh * ^a; ; 

(sec aj)^(sec aj)^,(cosec aj)^ ... sec a;tan2a;,seca;tan*a;, . . . 

203. The integration of J{a^ - x^), J{x^ - a?\ ^{a^ + x^) 
is required in the quadrature of the circle, ellipse, and 
hyperbola (§§ 50-54); thus for example, integrating by 
parts for the quadrature of the circle, 

(and, by transposition and division) 

= J^(a2-a;«) + Ja%in-i(x/a) (§ 50). 



= j(x- 
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Similarly fsj^^ ■" ^ • ^ "" iS)^^ 

SO that taken between the limits a and ^, at which 
^{a — x.x—^) vanishes, and between which it is real, 
the integral is ^7r(a—^)^. 

More generally, denoting ax^ + 260: + c by i^, and 
integrating by parts 

y^ t-rn .ax + b ,„ ,b^^ac f dx 

and the result has one of the three forms of ex. 28, p. 80. 
When a is negative and V^ — ac positive, the value of 
x is restricted to lie between the roots, a and fi suppose, 
of the equation ii = 0, for R to be positive and ^R real ; 
and now, taken between the limits a and ^, 

Suppose for instance that the area of the ellipse is 
required, when it is given by the general equation of the 
second degree (§ 13) 

ax^ + 2 hxy + by^-\-2gx-\- 2fy + c = 0. 

Solving this equation as a quadratic in y, 

hx + by+f=^{{h^-ab)x^ + 2(fh-bg)x+P^bc}; 
so that (§ 49) 

%i - 2/2) = V{ ('^' - «^>' + 2(/A - 6^)0: +/2 - 6c} ; 
and the area of the ellipse 'i^j{yi — y^dx, taken between 
the limits for which the quantity under the radical is 
positive; and therefore the area is tt A/(a6 — A^)^, where 
A denotes the discriminant (ex. 9, p. 139). 

This result must be multiplied by sin o), if the co- 
ordinate axes are inclined at an angle cd. 



/: 



a 
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204. Foi'mulas of Reduction. 

A foi^inla of reduction in the Integral Calculus is a 
formula, obtained in general by integration by parts, by 
which one integral, say u^ is made to depend upon a 
simpler integral, say Un-i> or Un-2\ then by successive 
substitution in the formula of reduction, we finally arrive 
at an integration which can be effected. 

Suppose, for instance, Un=J{^i^ O^dd; integrating by 
parts, with u= (sin 6)^"^, do = sin QdQ, v= —cos 0, 

ttn=y(sin eTde=/{sm 0)'^-isin OdO 

= -(sin 0)~-icos e + {n-l)/(sm 0)»*-2(cos 0)^6 

= +(n-l)an-2-(n-l)Un; 

or nUn= +(n-l)Un-2f 

1 -/i-l 

Un = — (sin0)""^cos0H !(„>.>, 

n n 

a forynula of reduction. 

Similarly 

/cos OYde^^siii 0(cos 0)""' + ^^ /cos OY-MO, 

another formula of reduction. 

Taken between the limits and ^tt, 

(sin erde=-^y (sin ey-^e, 



(cos 0f rf0 = ^^^^^^y (COS er " '^^^• 



First suppose n is an even integer 2m ; then 

/ \C08 / 2?7l ^ \C08 / 





2m -3 2m -1 A^'^PeY""' " ^^g 
~2m-2 ' -im y \co. / 







\ 
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and so on, finally being 



cos ' 2.4.6 ... 2m 2' 



Secondly, suppose n an odd integer 2m + 1 ; then 
/-J^ /sin \2m+l 2.4.6*... 2m fh'^ ""^ ,^ 

y \co8 / 3.5.7 ... 2m4-l^ co» 





2.4.6... 2m 



3.5.7 ...2m+r 

These are called Wallis's Theorems, and arS of great 
practical use in the Integral Calculus. 
Both cases are included in the formula 

J \cofi / 7i(n — 2)(tl — 4)...^^ ^ 

If we put cos Q — sech u,or 6 = gd u, then cZ0 = sech udu ; 
(cos eyde =y (sech ^y^+^cZu ; so that 



y/°° 1 xo,«_Li 7 1 .3. 5 ... 2m— 1 TT 



/^" v N2,«j 2.4.6... 2m-2 



When 771 or ti is not an integer, these integrals depend 
upon the Gamma Function, defined in the next article. 

Examples. — Prove the following formulas of reduction, 
and evaluate the integrals for 7i=l, 2, 3, 4, ... . 

(1) Un =^/x^e^dx = x^^ — nfx'^ " H^dx = x^^e^ — 7iu„ _ i. 

(2) Un=y(loga;)^c?aj = aj(logfl:;)^*— TiUn-i. 

(3) itn=/x'' (log x)''dx = --7w~ — — TT- 

\ / n ^ \ G f 771 + 1 771 -hi 
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(4) u^=/x-J{a^-x-)dx L_^ + _^a%„... 

(5) Un =/c"^(2aa; — aj2)cfe 

_ a;^-H2aa;-a;^)^ , 2n+l 

(6) Un=y(a2-iz,'2)"-i= {a:(a2-rc2)~-i + (2^-l)a2wn.i}/2n. 
,^x /" cfqg 

^ 1 1 Ax+B 271-3 J.Un-1 

2n-2 AG-B^ (Ax^+2Bx+C)'''^'^ 271-2 AC-B^' 

/ox ^ f ^^ 5 . np-p-\ Un-J 

^ ^ "" J (xT^^-a'^f' (n-l)pa^(x^+a^)'»-i'*" (n-l)p a^ ' 

^^^ ''''V(x-aYJ{x--h) 

1 1 sJ{x-}>) 2n-3Ttn.i 
n — 1 6 — a (a:— a)"'^ 2ti — 2 6 — c* 

(10) Un=/x''dxlJR 

= {aj""^>v/-R— (2n — l)6un-i — (n — l)ci/n-2}/^a. 

(11) Determine X, ilf, iV in the formulas of reduction 

(13) u» = /(tan x)"£Za; = ^ ^_ Un-2. 

(14) Un=y(cota;)«da;=-^-^^j^ iin-2. 

, , „. /T ,^ , (sec a:)" - ^tan x . rt — 2 

(15) Un=/ (sec X)^dx = - :i ^i-^n-2. 

/i/>\ /" \,ij (cosecaj)''*"^cotir , 7i— 2 
(16) Un= / (cosee xy^dx =^ —^ --r h-;;^ttn_2- 
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(17) Un=^ (vers x)''dx= ^—^ — ^- — + ^ Un-i. 

/io\ f.r ^"^ ' '^^'^"^ n(n-l) 

(18) Un=/x*^cosmx = — smmaj+ — s-cosmaj — ^^ — o-^i^n-s- 

/i /^\ /'•v../ w 7 e^(cos a;)'*-Vp cos aj+ti sin a;) 

(19) Un =/eP^(coB xYdx = — i^- — J^2 ^ 2 ^ 



(20) V(,n, n) = /(sin aj)'Xcos a;)'^da5 



n(n^l) 



_(sin aj)^+Xcos x^-^ n — 1 
(sinaj)"*"Vcosa;y*+^ . m — 1 

or — ^ ^ p ; U(m-2,nh 

m+n vi+n 

(7n-l)(sinxy'''^\cosxy'-^-{n-l)(iimxy''-\cosxY-^'^ 
' (m + n)(;m + 11 — 2) 

(m -- 1 X^^ - 1 )^(m-2. n-2 ) 

(m + '7i)(7ii + ?i — 2) 

(21) /(Hhixyicosxr^+'dx - 7- ^,\t'^_^o^"\ lo ^1V 
\ /^ \ / \ / (m+l)(m + 8)...(m + 27i+l) 



(^''> "-t/ti 



cZ0 b sin 



(a + 6 cos 0)'» ('^i' - 1 ) (tt2 ~ 62)(a + b cos 0)^ - 1 

(2'yi~3)at<n-i-(->i~2)ttn- 
(7i-l)(a2-62) 

Investigate formulas of reduction for the integration of 
(sin~^.r)^ (vexs''^xy\ (cosh x)'^, (sinh xy\ (tanh a;)**, (sech xY, 
(sinh'^aj)^ (cosh-^ic)'*, sin ^u?, sinh ^o;, (a + 6 cosh u)-", 
(a+6cos0 + csin0)"'*, (a + 6coshu + csinh it)"'*. 
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205. Definite Integrals, 

Any integral taken between definite fixed lii\ 
called a definite integral (§ 45), but the name is tfsed 
especially for a large class of integrals which can only be 
evaluated between certain limits, and of which the in- 
definite integrals (§ 39) cannot be found. 

The subject of Definite Integrals is a large one, and 
would lead us beyond the scope of the present treatise ; 
we may however consider a few elementary definite 
integrals, which are of frequent occurrence, and great 
practical use. 

Consider first the definite integral 

J e-^x'^-'^dx, ox y {-log z)''-^dz=y (log I'/zy-hlz, 



if e~^ = z; this is called the Gamma function, because it 
is usually denoted by Tn. 
By integration by parts 

(\ OO /'OO 

^^-x^n^i\ +(n-\)/e--''x''-^dx = (n-l)T(n-l)] 



so that if 'Ji is a positive integer, 

r7i = (n-l)!ri = (n-l)! 

since Fl =/€ '^dx = 1, also TO = x . 



Changing x into x^ gives 

•* 00 

Tm^l/er'^ x^'^'-^dx ; 



so that we may express the product of Fm and Tn by 
F?n,Fn = 4 /e - * V'« - \lx le - y'-f' - ^dij 
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/ 












(on changing to polars, with a* = r cos 6, y = r sin 0) 
= 2r(m + n)/(sin 0)2«* - ^ (cos 0)2» - HQ. 



The definite integral 

^ ^ 

on putting x = {smQY\ and is called the First EvZeHan 

Integraly and denoted by B (m, n) ; so that 

r7/ir')i = r(m + ^)-B(m, n). 

Also, on substituting 
1 a-aj = (a-^)sin20, x- 13 = (a- P)co8% 

^ y(a-ic)'^-i(a;-i8)^-id:aj = (a-/3)^'*+~-i5(m, n). 

This integral has been evaluated for positive integral 
values of m and n (p. 424) ; and for fractional values it 
depends on the tabulated values of the Gamma Function, 
by means of the above relation ^ 

B(m, 7i) = TmTn/r{m + n). 

When m + 7i = l, then (§ 189) 

Tm r(l -m) = 5(m, 1 -m) = 2/(tan0)2'^-ic?0=7rcosecm7r. 



Putting m = i, gives rj = V7r, and Tf = jr^ = ^^^tt ; 
and now if Tm is calculated from m = 0*5 to m = l, then 
Tm is known from m = 0'5 to m = 0; and thence Tn or 
logPn can be tabulated from n = l to n = 2 (Bertrand, 
Galcul Integral, p. 285), and the graph of Tx can be drawn 
for all positive vahies of x. 



2n\ n 
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The substitution of x=^zl{\+z) or 1/(1+0) shows that 

As an application, let us determine the area A of the 
positive quadrant of the curve x^'\-y'^^a^\ then 

A =/(a«-aj»)icZaj = ^a2/Jsin 0)^l(cos &pHQ 


(on substituting x^^ = a^sin^Q) 

=«X-.-+i)=-rir(i4-i)/r(^+i) 

"ti \n n J n n \n // \n J 

^ n 

For instance, if n = 2, we obtain the area of the quad- 
rant of the circle Xird?. 

Also 5 = y = ia(r ^)y (rir-). 

By integrating, as in § 134, throughout the space in 
which X, y, z are positive, but 

(a;/a)»«+(2//6)^+(^/c)^-l 
is negative, we can show in a similar manner, that the 
volume F of this space within the surface is given by 

_^=r(i+i)r(l+i)r(i+i)/r(i+l+-Vi) 

abc \m J \n J \p /i \m n p / 
while 

agreeing with § 134 when we put m = n=p = 2. 

Again, the length of a branch of the curve r'* = a'*cos nO 

(§ 172) can be expressed by Ma"* a)' 
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206. Bemovlli's and Eulers Numbers as Definite 
Integrals, 

Making use of the relation, with mu = x, 

- 

then, with the notation of § 112, 





2prji 2^"^ /^^ 

so that Sj, = TTT— ^ = TTn — TvrT - /cosech u . u^ ~ ^dtt ; 

^ 2''— 1 (2^— i)rp y 

2n r^ 

and £n = /22n,x\ 2n /cosech u . u^** - Hu ; 

so that, on putting u = Ittv, 

cosech ^TTV . t;2^ - Hv = 22«(22« - l)5„/27i, 



the coefficient of x^^~^/(2n — l)\ in the expansion of tan x, 
and this is always an integer ; namely 1, 2, 16, 272, .... 
Similarly (§§ 112, 113) 





which, since En is an integer, shows that this definite 
integral is always an integer, namely 1, 5, 61, 1385, 

Put sech ^TTf = cos 0, \irv = log(sec + tan 0), 

and we find (S 118) 
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{log (sec ^ + tan <t>)YM<j> = {^-wf^'-^^En = 2(27i)! f^2n+i, 



J {log (sec ^+ tan 0)}2"-id^ = 2(2n-l)!f/2n. 



Thus, for example, the area between a branch of the 
curve m6 = gd {rja), or rja = log (sec m 6 + tan mO), 
and its asymptote is 

rirrim 

hay (log (sec mO + tan me)}^de 



207. Differentiation and Integration of a Definite 
Integral, and Applications, 

Denoting the indefinite integral with respect to a? of a 
function i(x, c), involving a parameter c (§ 105), by {^{x, c), 
then the definite integral 

i - liix, c)dx = f i(a, c) — f i(6, c) ; 

b 

and now, if a, 6, c are functions of some variable f, 
dJjdl^ da,'^ db dl dc 
di~da di'^db Rt'^dc dt 

b 

of which a geometrical interpretation can easily be 
constructed. 

We can also integrate / according to a similar rule. 

By means of this method of difierentiation and inte- 
gration we can show the connexion between different 
classes of integrals, as already done in §§ 193, 194, and 
also verify the solution of certain differential equations, 
when the solution is given as a definite integral, as in § 184. 
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Thus, if a and /3 are positive, and /=— 7— ^tan-^A/-, 
r\^ d6 1 .V acos^e-^sin^e y^ ^ 

P__A^ ^ 1 ^^^.i/ gcosh^^-^sinh^ N"^ ^ 

y acosh^itri-^sinh^ii. ^^Jia^) Vacosh^u+jQsinhWo"" ' 

Ind /i^ ^e _ C-1)Y 3 , 3Vr 

y(acos20+^sin20)~+i"" n! \9a"*'3/3/ ' 



A 





00 



Cil6 



^(idlYLAY/. 



fi 



(acosh2i^+^sinh2u)~+i ti! \3a 3)8/ 
We can verify that the definite integral 

cos(29V — X sin i;)rfi;, 

which we have denoted by irJp(x) in § 184, and called 

the Bessd function of x of the jpth order, is a solution of 
Bessels differential equation given on p. 185. 
For dividing the definite integral into two parts 

y = /cos pv cos(aj sin v)dv, and z = / sin pv sin(aj sin v)dv ; 



and using accents to denote differentiation, 

xh/'^ + xy' + (x^ —p^)y 

— I { — flj^in^v cos jDV cos(a; sin v)^x sin v cos py sin (aj sin v) 
+ (x^ —pr')co%pv cos(a; sin v)} dv 

= {x cos V cos pv sin(aj sin v) —p sin jpv cos(ic sin v)}^ 
=p8inp'jr; 

while ccV + xz' + (x^ -'P^)z 

= { — X COS V sin pv cos(a: sin v) +p cos pv sin(a; sin v)} 
= — fljsinpTr; 

and therefore 2/ and z each satisfy Bessel's equation if p 
is an integer ; and y or z vanishes as p is odd or even. 
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The student may verify in a similar manner that if Pnip) 
(p. 314) satisfies the difierential equation of ex. 6, p. 276, 

irPniiJi) =j{tJi - s/{p? - 1 )cos ,j>Yd4 




=/{l + s/(/^^-l)C08 ylr)-^-H-^, 





the one integral being transformed into the other by (§ 200) 

Putting yti = cos u, or cosh v, according as /i is < or > 1, 
S+ cotu||+«(«+l)P=0, 

,g+cotht;^-n(.^+l)P=0; 

J 

so that the solution may be written 

7rP„(cos u) = / (cos 16 — i sin u cos <j>Yd<l>, 



or 7rP„(cosh v) = /(cosh v — sinh v cos 0)"d0. 



Another solution Qn is obtained by putting 

Qn(M)^{M + x/(M'-l)cosh 0}-»- W, 

u 

/•coth~V 



by means of the substitution (§ 200) 

{/^+V(M'-l)cosh0} {M-V(M'-l)coshVr} = l; 
and QnC^) = ^ when yti = 0. 

Example, — Prove 



or 
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208. Mr. W. M. Hicks has shown {Phil Trans., 1881) 
that similar functions are required in the solution 
of "^^¥=0 (§ 146), in connexion with the anchor rings or 
toroidal surfaces, generated by the revolution of the 
dipolar circles of § 175. 

With the columnar coordinates y, 6, x (§ 132) round 
Ox as axis, we obtain the transformation 

^ 'dx^'^y dyV^'dyJ'^y^ d0^ ' 
or, on putting V^rfry'^, 

and, on changing from x, y to new variables u, v, con- 
jugate functions given by u+iv = i{x+iy) (§ 133) 

With the dipolars of § 175, 

x+iy = c tan \{u + iv\ 

1 /dx\^ /'^y\^ <^ 

J~~\du/ KduJ "~ (cos u + cosh v)^""^ * ' 

Then, if t/r varies as cos nu cos mO, 

— '^^V=y-hm\\h\ ^ —n^—(m^—\)\}r cosech^k 

Further, if we put \}A = ^{csinhv)P, —'^^V= 
{cos u+ cosh vy^f-^P , d^P , , , 3P , u2 3^^ ■ 1 dN 



c« 
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If we suppose P = P„cos7iucosm0 is a solution of 
y2y__o^ \^Qii p^^ as a function of v only, satisfies the 
differential equation 

d?P dP 

-^-o + coth v-j — (n^ — i)P —m^P cosech^t; = ; 

and this equation, in the particular case when m = 0, is 
the same as that satisfied by the zonal harmonic (p. 276) 
with 71 — I for 71 ; so that writing G for cosh v or jul, and 8 
for sinh-y or ^(jj.^ — l), the solution is APn+BQn, 

where ^^n^^fj _^ g^^^ ^y,^^^ ^^^y (G + ^^cosh^)-+i- 
For, as in ex. 6, v., p. 276, 

dG J (C+/Sfcos0)~+^ 

; ^^y ((7+scos0)~+^ ^ 
= (7i+i)(Pn+i-c'P„); 

while taking the form TrPn = / {C—S cos <f>)^~^d^, 

' 

'S^^=(»t-i)((7/'«-P«-x); 

and thence ^(s^^») = (.^-i)P„; 

and similarly with Q^. 

209. In the general case, when on is any integer, we 
can verify that the solution of 

is p=ap:+bq:, 

2e 
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where 

p,» _ /**• cos m^fp ^,» _ /''*' cosh m<ltd<l> 

" ~J{C+S cos ^)»+i' '^^•' 'y(C+Scoah<f,)"+i' 
For now, taking Qj,", we prove that 

(s«^-m(7)Q:+(«-m- i)<SQ:+i 

- r 'S8ipK^+l)^+<7sinhm0 \"_^ .^ . . 

~ t x(C+Scosh^)»+i 1^ -«.»*»«•<«. 

and thus, as in § 87, 

= (_l)-(n-^)(ti-»)...(n-m-i)Q;r. 
Similarly 

(-^2^+ mCf)e + (n +m- J)S« " ' 

_ / iSsinh(m— 1)^+ Osinh m^ \ *'_ » 
~\ 7r(0+Scosh^)»+i Jo"" ' 

and thence ^^Q^ 

= -mQ:-mC^,Q:-in+m-^)(^Q:-^+S^^Q:-^) 

thus verifying the differential equation; and the procedure 
is the same with P"*. 

The substitution tani0 = ^W"v- . v ^^^ will give 

^ cosh ^u smh <f) ° 

f' dd _ n jm _ . Q . 

and generally if we denote / rr? j^n by xi^n. then 

° ^ (C/— COS0)* 
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^--^-(n^\\R,.^.-Qn^) = \ (C-cosg) |o=^' 

it being an integer ; so that R^ satisfies the same equa- 
tions as Qn ; and since R^ = V^^o* therefore R^ = s/^Qn- 
By difierentiation with respect to (7, 

*.^.t...(m jy (O-cose)"*^* 



= 7rV2(^i-i)(n-f) ... (7i-m~i)e 



((7 + /Sfcosh^f+i- 

When the dipolar system of § 175 is revolved round 
Oy, we shall find in the same manner, on putting 

V= \[rx " *, and x/r = ^(c sin u)P, 



c^ 



and now if P=P^cosnvcoa7nd is a solution of ^^V=0, 
then i^, as a function of u only, satisfies the equation 

^ + cot tt^ - {n^ + i )P - m^P cosec% = 0, 

of which the solution is obtained from the preceding results 
by replacing cosh v and sinh v by cos u and sin u. 

In the particular case of m = 0, the order of the corre- 
sponding zonal harmonics would be - J + in, and therefore 
imaginary; these are called MefUer functions. 
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210. LcLgTange*8 Theorefm. 

We can establish Lagraiige!s Theorem (§ 150) by the 
rule for the diflFerentiation of a definite integral ; for, if 

so that differentiating again n— 1 times with respect toi», 
and denoting 

a 

we find that 

while RQ = iy — ia ) so that 

which is Lagrange's Theorem. 

211. Integration conMdered as Summation. 

Hitherto we have considered Integration as an opera- 
tion to be performed by the reversion of Differentiation ; 
but now, by the aid of Taylor's Theorem, we can exhibit 
it in its fundamental conception as the Summation of a 
number of infinitesimal elements, and show that the 
result is the same as before. 

For let any function ix, between the fixed values a 
and 6 of x, be considered for the consecutive infinitesim- 
ally graduated sequence of values a, x^, x^, ... Xr, ... x^, 6, 
of the variable x. 
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Then by Taylor's Theorem (§§ 114, 118), using the 
symbol .4^r to denote the difference aJr+i— aJr, 

tx^^ia -Aa fa =\^a?' £"( a+0 ^a ), 



ft — f aJn — AXni'Xn = ^ AaJn^"(aJn + 0n AflJn). 

Adding these equations, and using the symbol 2 to 
denote summation for all integral values of r from 1 to n, 
f 6 - fa - XAx4% = 2iAx^T(Xr + OAxr). 

Now if denotes the greatest of all the evanescent 
quantities, such as AXri\Xr+drAxr), then 

ib-ia-^i'XrAXr is less than ^o^Xr+i-Xr) or ^o(b-a) ; 
and this quantity ultimately vanishes ; so that replacing 
2 and A by y and d, as in § 44, 

f6-fa-yf4cfe = 0, 

a 

which establishes the fact that Integration considered as 
Summation is a process the reverse of Differentiation. 

For instance in Quadrature (§ 41), we determine the 
arithmetic mean of the ordinates of a curve, and call this 
the mean ordinate ; and now the area between any two 
ordinates is a rectangle of the same breadth, and height 
equal to the mean ordinate. 

Suppose however that M, the geometric mean of the 
equispaced ordinates of the curve y=ix between x = a 
and x = b 18 required ; then M=lt !i/{ix^ . fajg. . .fa^r- • .fa^n), 

and log JIf = It -2 log fee =lt 2 log ixAxlQ) - a) 



= /log ixdxjQ) — a), 
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so that M = exp-^ /log ixdxjQ) — a) J-. 

a 

As an application, prove that the geometric mean of 
the ordinates of the curve 2/ = sinaj, from to tt is J. 

212. Approximate Quadrature. 

Various rules are employed in practice for determining 
A the area of an irregular figure, when a Planimeter is 
not used. 

The simplest is that usually employed with an Indicator 
Diagram ; the breadth of the closed contour is measured 
at a number of equidistant ordinates which divide the 
area into & number n of strips, and the arithmetic mean of 
these is taken as the mean breadth; more strictly we take 

^=Khyo+y 1+2/2+ '"+yn-i+hyn\ 

where the y*s denote the breadth of the figure at the 
different ordinates, and I denotes the spacing of the 
ordinates: this is equivalent to replacing the area by 
trapezoids, bounded by the same ordinates. 

In Simpson's rule a higher approximation is obtained 
by taking the curve through the tops of three equidistant 
ordinates as a parabola with axis parallel to the ordinates; 
and now if ym^h Vm* 2/m+i denote three consecutive 
ordinates and I their spacing, the parabola is taken as 
given by y = a+bx+ ^cx^ ; 

and if the origin is taken at the foot of the middle 
ordinate, as in § 77, 

a = ym, ^bl = 2/m+l - ym - 1, CP = y,n+l - ^ym + ym - 1- 

The area of the two strips is now given by 

/ydx = 2al + ^cP = ^l(ym+i + ^ym + 2/«i - 1) ; 

-I 
and A = J%o+^2/i + 22/2+%3+ • • • + ^2/2« +2/2n+i), 

which is called Sivipson's Rule. 



INTEGRATION IN GENERAL. 439 

We may replace the boundary by the curve of the 
third degree y = a+hx + \cx^ + \doi?, 

passing through the tops of four equidistant ordinates 

2/«i> 2/»»+ii yw+2, ^wi+s; and now the area of the three strips 

rzi 
\jdx = U{a + m + Ul^ + f6?3) 



fi 



= %Kym+ Sym+i + ^ym+2 + ym+s)' 
It is convenient to take an odd number of ordinates 

and an even number of strips, and a general formula for 

approximate quadrature, with 271+1 ordinates and 2n 

strips, has been given by Newton. 

With 7 ordinates and 6 strips, the expression for the 

area is, very approximately, 

called Weddle's Rule. 

In these expressions for the area, the y*s may represent 
breadths as well as ordinates ; or they may represent the 
areas of parallel equidistant sections of a solid body ; for 
instance, horizontal and vertical sections of a ship. 

Simpson's Rule may also be used for finding the centroid 
and moment of inertia of the area, by putting 

Ax=/xydx =il(xQyQ+ 4<x^iyi+ 20^22/2 +...)» 

Ay^Ay'dx=U( yo^+ ^yi'+ ^yi+-')\ 

Ah^^/lyHx = il{ 2/oH ^y^'+ 2y,^+...\ 
Ahy^ =fji?ydx = \l{x^yQ + ^lyi + ^^^y^ +...)• 
As a numerical application, calculate these quantities 
for the horizontal section of a ship, given the half 
ordinates from the median line as "5, 6, 10, 12, 12*4, 12*5, 
12-5, 12-5, 12-4, 12-3, 11, 8, and -5 feet; the ordinates 
being spaced at 12 feet. 
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213. Differential Equations, 

In Chapters III. and V. the principles of Differentia- 
tion have been illustrated in the formation of certain 
important Differential EquatioTis from their primitive 
relations ; but now it is requisite to reverse this process, 
and to discover the primitive relation or solution where 
possible, when the differential equation is given. 

The complete subject of Differential Equations would 
soon lead us beyond the scope of the present work, and 
the reader is referred to the treatises of Forsyth and 
Woolsey Johnson ; at the same time, however, it is 
possible to indicate the procedure in a few simple cases, 
often required in Dynamics and Electricity. 

The simplest case of frequent occurrence is the linear 
differential equation with constant coefficients^ of the 

where the A's are constant, and Fis a given function of x. 

Using D to denote the operation d/dx, we may write 

the equation F(D)y=^ V, (1). 

where F(D) is a rational integral algebraical function of 
D ; and now the solution will consist of two parts, one 
of them being any particular solution of this equation 
(1), while the other, called the complementary function, 
is the general solution of the equation 

F(i))t/ = , (2). 

We have shown in § 68 that F(i))e^ = 6^F6, so that 
2/ = e^ will be a solution of (2) if F6 = ; and if the n 
roots &!, 62* ^8» ••• ^f ^^^^ equation F6 = are all real and 
distinct, the complementary function will be 

y = B^e^^^ + B^e^^ + . . . + 5^6**^, 
where B^, B^, .,,, Bn are n arbitrary constants. 
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If m roots are equal to 6, the corresponding terms in 
the complementary function are, as in ex. 3, p. 143, 

((70 + C^aj + . . . + (7^_ia;"»-i)e^. 
If a pair of roots are conjugate imaginary, a±ip, the 
corresponding terms, by De Moivre's theorem, are written 

e^{A cos^x + B sin fix) ; 
and if these imaginary roots are repeated m times, the 
corresponding terms are 

€^^a^(AtOOS fix + J5^in fix). 

In determining the particular solution of (1), we write 
the equation symbolically 

where (FD)"^ can be resolved into a series of partial 
fractions of the form A{D-b)-\ or B{D-by\ 

Consider the separate terms of V, of the form 
05**, e^, c^, sin mx^ .... 

(i.) (FDy^x^ is foimd by expanding (FD)""^ in ascend- 
ing powers of Z) as far as D^, and performing the 
differentiations; terms beyond D** may be neglected, 
since i)~+iaj'*=0, .... 

(ii.) (FD) - ^e^ = e«*(Fa) - \ 

and (FD) - ^ef^x = e<^{F{D + a) } - Hx, (§ 70). 

(iii.) Thus,if F6 = 0, 

(Z)-6)-'•e^=e^i>-^l = e^(Co + (7la; + ... + a.la;»"l + ^j), 

since i)"*" represents the operation of integration, re- 
peated r times. 

(iv.) {D - b) -'•cos mx=(D'{- by{D^ - 6^^ -''cos mx 

= (D + 6)'-( - m2 - 62) - r^os ma;, 
(§ 68); and so on for (i)2 + m*)-''cos mx, ..., by employing 
the exponential values of cos ma?, ... (§ 111). 
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Consider for instance the differential equation 

d'^tC doc 

--rr^+ 2n cos P'-Ti +n^x=Ecos pt, 

for the forced vibrations, due to E cos pt, of a system of 
which the natural motion is given by the complementary 

function a: = ae " **' ^^ cos (oit sin 13+e) 

(ex. 6, p. 139), representing S.H.M. (§ 103) of period 
2Tr/{n sin )8), with a modulus of decay n cos fi. 
The particular solution is 
__ E cos pt 

^ ■" i)2 :^2nc6s~0TD + n^ 

"" (i)2 + ^2)2 - 4n2c082^ . D^ ^ ^ 

_ (n^ -'P^)E cos pt + 2np cos ^ . Eain pt 

"~ (7l2 — p2j2 ^ 47i2p2Qos2)8 

_ Jgcos[j9^— tan~^{2tij}C0SjQ/(n2— p2)|j 

>v/(^' + 27l22)2cos 2/3 + p*) 
representing a S.H.M. of amplitude 

E/J(n^ + 2nY(^os 2^ + p% 
and change of phase tan " ^ { 2np cos ^l{n^ — P^) } • 
Examples, — Solve the differential equations 



(i- 
(ii. 

(ill. 

(iv. 

(v. 

(vi. 

(vii. 
(viu. 

(ix. 
(x. 

(xi. 



y" — 02/' + 62/ = a:;^ + e^ + sin 4aj + xe^ + 052^**, 

y" + 42/ = sin X + sin 2aj, 

y" — 62/' + 92/ = cosh ^x + x cos 2aj, 

2/" + 2y' + oy = e^sin 2a; + cosh x cos 2aj, 

2/'" + y"'-'y'—y = a^ + cosh a?, 

2/'" — Sy" + 162/' = sinh 2a; + cosh 4a;, 

2/"' + 2/' + IO2/ = cosh 2a; + sinh x sin 2a;, 

2/'" - 262/' - 6O2/ = sinh 6a; + (il + i^ cosh 3a;)cos x, 

y"" + 2y" + a; = sin 2a; + sin x, 

y"- _ 22/'" + 32/'' - 22/' + 2/ = a^ + € - **(il + JS cosy 3a'), 

y"" + 322/' + 482/ = a; sinh 2a- + cosh 2a; cos (2^2a;). 
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214. With variable coefficients P, Q, iJ, ..., given 
functions of x, the general linear differential equation of 
the first order is written y'+Py^Q. 

To solve this, we can proceed as in § 88, and put y = ux\ 
when u'v + u{v + Pv) = Q ; 

and now if v is chosen so that 

i/+Pf=0, or v=exp(-/Pflx), 

then u = Q exj>y^Pdx, 

and u =y^Q(expy^Pdx)(Ix + C, 

y = ia'=exp(-/Pdx){/Q(exf/Pdx)(Jx+C }- 
Thus, for instance, the solution of the equation 

L ,^ + Ri = Ecospt, 
at 

connecting It the resistance in ohnis, i the current in 
amperes, E cos 2>t the variable electromotive force of a 
dynamo in volts, and Ldildt the counter electromotive 
force of induction, so that L is the inductance in secohms 
(Fleming, The Alter^iate Current Transfirrmer) is 

i = e-^^'^i/Ecospt . e^'hlt+C) 
__RE cos pt+pLE sin pt.p _j„,j^ 

_ ^ cos {y>^ - tan - HpLjR) } ri^^EHL 

But considered as a differential equation with constant 
coefficients, Ce~^^^ is the complementary function, while 
the particular solution is, as before, 

. Ecospt R—LD jp . 
_R Ecospt+pLEampt _ E<iOB{pt—i&n-\pLIR)} 
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215. The general linear differential equation of the 
second order (§ 88) y" + 2Py +Qy = R, 
on putting y = uv, and choosing v so that 

becomes u'^+In=S, 

where I^Q^-F'^-P', S=R/v; 

and now, if we can determine the complementary function 

Aui+Bu^y then, as in § 88, 

and the particular solution is 

Cu = Ui y^Su^dx — u^ySu^dx. 

Thus, as on pp. 183-185, Bessel functions are required 
in the solution of Biccati's equation, in which /=Aja;**; 
and the student can easily verify that 

= Axf ^^j^(maj«sin 0)(cos 0)*cW 



+ J5y^"'^^(ma;'*sinh 0)(cosh 0)'^cZ0 



is the solution when /= ±m^'n.^a;^'*~^ with certain restric- 
tions on the value of n. 

Again, zonal hurmonics (pp. 276, 314, 431) give the 
solution when 

Iz=(n+^)^+l cosec%, or — (ti + J)^ + 1 cosech^i; ; 
and the toroidal and Mehler functions of §§ 208, 209, when 

1= -~n^- (rn? - J)cosech2y, or - ti^ _ (m^ - J)cosec^iC'. 

But no general method has yet been discovered for the 
solution of the linear differential equation of the second 
order in the canonical form 

u''+lu = 0, 
where J is an arbitrary function of x ; nor a fortiori 
for differential equations of a higher order. 



u 
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216. Sirnultaneov^ differential equations vfii^ 
coeflBicients, as those given in exs. 3 and 4, p. 219, may 
be solved by assuming x = Pe^\ y = Qe^\ oxx=P sin(7i^ + e),. 
y = Q sia(nt + c'), and eliminating P/Q, and c, e', when an 
equation for the determination of 6 or 71 will be obtained. 

Thus for instance, in the theory of the Induction Coil, 
the differential equations 

L-Ti + M-T^ + Ri-^=E cos pty 

for the primary current i^, and the secondary current ig, 
aie solved by assuming 

i, = /iCOS (pt+e^), i^ = I^eos(pt+€2h 
and determining Z^, /g, e^, e^. 

The formation of partial differential equations has 
been illustrated in Chap. V., from which the solution of 
simple partial differential equations of the first order 
may be inferred ; but for further developments the 
reader must consult the Treatises on Differential Equa- 
tions mentioned above. 



